


UNIVERSITY 
OF MICHIGAN 


FEB 21 i395] 


ENGINEERING 
LIBRARY 











FEBRUARY, 1951 





CONTENTS 


The Minimum Drag of Thin Wings in Frictionless Flow ROBERT T. JONES 


Design and Performance of a Simple Interferometer for Wind-Tunnel Measurements.............. 
Harry I. ASHKENAS AND ARTHUR E. BRYSON 


Remarks on Low-Aspect-Ratio Configurations in Supersonic Flow 
P. A. LAGERSTROM AND M. E. GRAHAM 


Wind-Tunnel Tests for Temperature Recovery Factors at Supersonic Velocities W. F. Hrtton 


Structural Load Measurements on Complex Aircraft Components Using Strain-Gage Summation 


Circuits Maurice A. WALTER 


Conical Flow Methods Applied to Uniformly Loaded Wings in Subsonic Flow. ...SzyMouR LAMPERT 


Calculation of Thermal Stresses in a Wedge-Shaped Wing 
GEorRGE H. EISENHARDT AND WARREN M. ROHSENOW 


Buckling Load of a Stepped Column FREDERIC M. Hos.it 


On Some Special Problems in Linearized Axially Symmetric Flow 
BERTRAND DES CLERS AND CHIEH-CHIEN CHANG 


Readers’ Forum: ‘Generalized Complex Forces in Flutter Calculations,” by J. E. V. Nilsson and 
N. B. Langefors; ‘‘Vibrations of Rectangular Plates,” by Harold Lurie; “On Supersonic Wing- 
Body Center of Pressure,’’ by G. K. Morikawa and T. F. Coleman; ‘Airplane Pitching Moment 
of Inertia,’’ by Innes Bouton; ‘‘Base Pressures at Supersonic Speeds,’”’ by J. Leith Potter; “A 
Simple Graphical Solution of the Duhamel Integral,” by E. V. Laitone and J. T. Ahlin; ‘‘Aero- 
dynamics of Blasts,” by L. Ting and H. F. Ludloff 


Published by the 
INSTITUTE OF THE AERONAUTICAL SCIENCES, INC. 














JOURNAL OF THE 
AERONAUTICAL SCIENCES 





Editor: Hucu L. Drypen, Director, N.A.C.A. 
Associate Editor: Ropert R. DEXTER Managing Editor: Berneice H. Jarcx 


The Institute of the Aeronautical Sciences publishes monthly the JourNAL OF THE 
AERONAUTICAL SCIENCES and the AERONAUTICAL ENGINEERING REVIEW. 


The JouRNAL OF THE AERONAUTICAL SCIENCES prints papers of scientific and technical 
interest to Institute members and subscribers. In general, it provides an opportunity 
for specialists to publish articles disclosing new knowledge and new applications in the 
field of the aeronautical sciences. The JourNAL restricts itself to the aeronautical as- 
pects of advances in science and engineering, and papers printed have either been sub- 
mitted directly to the Journat or have been previously presented at a meeting of the 
‘Institute. 


The ArERonauTIcAL ENGINEERING REviEwW keeps members and subscribers informed 
regarding the contents of new books and current periodicals in the field of aeronautics, 
It reviews new developments in the aircraft industry and reports on progress from house 
organs and catalogs published by aircraft and engine manufacturers and other com- 
panies engaged in aeronautical activities. It publishes papers on applied engineering 
and aircraft production. News of Institute activities is published in the Review. 


Statements and opinions expressed in the JourNAL are to be understood as individual 
expressions of the authors and not those of the Institute. 





SUBSCRIPTION RATES 


Journal of the Aeronautical Sciences 


Published Monthly 


Foreign Countries Including Canada 
United States and Possessions (American Currency Rate) 


Single Copies ‘ Single Copies 
Aeronautical Engineering Review 


Published Monthly 


a Countries Including Canada 


merican Currency Rates) 


United States and Possessions 


Single Copies Single Copies 


Notices of change of address should be sent to the Institute at least 30 days prior to 
actual change of address. 
Manuscripts for publication, proofs, and all correspondence should be addressed to the 
Editorial Office of the Journat. 
Correspondence regarding subscriptions may be addressed to Publication Office, 20th 
and Northampton Sts., Easton, Pa., or to the Editorial Office of the JourNAL OF THE 
AERONAUTICAL SCIENCES, 2 East 64th Street, New York 21, New York. 

Copyright, 1951, by the Institute of the Aeronautical Sciences, Inc. 


Entered as Second Class Matter at the Post Office, Easton, Pa., May 1,1937. Acceptance for mailing at a special rate of 
Postage provided for in the Act of August 24,1912. Authorized April 29, 1937. 

















JOURNAL OF THE 
AERONAUTICAL SCIENCES 





VOLUME 18 


FEBRUARY, 1951 


NUMBER 2 





The Minimum Drag of Thin Wings 


in Frictionless Flow 


ROBERT T. JONES* 
Ames Aeronautical Laboratory, N.A.C.A. 


SUMMARY 


The assumptions of the thin airfoil theory are found to pro- 
vide certain necessary conditions for the minimum drag of airfoils 
having a given total lift, a given maximum thickness, or a given 
volume. The conditions are applicable to steady or unsteady 
motions and to subsonic or supersonic speeds without restriction 
on the plan form. The computation of drag and the statement of 
the conditions for minimum drag depend on the consideration of 
a “combined flow field,’’ which is obtained by superimposing the 
disturbance velocities in forward and reversed motions. 

If the plan form of the airfoil and its total lift are given, it is 
found that, for minimum drag, the lift must be distributed in such 
a way that the downwash in the combined field is constant over 
the entire plan form. If the plan form is given and the thickness 
of the airfoil is required to contain a specified volume, then the 
thickness must be distributed over the plan form in such a way 
that the pressure gradient of the combined field in the direction 
of flight is constant at all points of the wing. A specification of 
the thickness along some line drawn on the plan form is found to 
lead to the requirement that the gradient of the pressure vanishes 
on either side of this line. For the drag to be a minimum with 
respect to small changes in the plan form, the foregoing conditions 
must extend continuously for a small distance beyond the edge of 
the plan form. 


INTRODUCTION 


he DRAG OF A BODY moving through a fluid is at- 
tributable partly to the action of tangential, or 
friction, forces and partly to the action of normal pres- 
Sures. With well-streamlined bodies, the friction forces 
are ordinarily confined to a relatively thin boundary 
layer adjacent to the surface, and in such cases that 
part of the drag arising from the normal pressures can 
be determined on the assumption of a frictionless po- 
tential motion in the region outside the boundary layer. 

A practical way of minimizing the pressure compo- 
Ment of the drag is to make the body slender and its 
angle of attack small. Many of the results of the mod- 
: Presented at the Aerodynamics Session, Annual Summer Meet- 
ing, I.A.S., Los Angeles, July 12-14, 1950. 

* Aeronautical Engineer. 


ern airfoil theory, such as the theory of the induced 
drag of airfoils, are based on the assumption that the 
velocities imparted to the air are small in relation to the 
velocity of flight and, hence, are limited to cases of thin 
airfoils and small lift coefficients. 

In the present paper, certain necessary conditions 
for the minimum pressure drag of thin airfoils are de- 
rived. The analysis makes use of, and extends, certain 
“reversed flow’’ theorems originally derived by Hayes! 
and von Karman? and is based on the idea of super- 
imposing the disturbances in forward and reversed 
motions originally advanced by Munk.* Although 
the analysis is restricted to motions involving small dis- 
turbances, the conditions for minimum drag are found 
to be applicable in a wide variety of circumstances of 
such motions. Thus, the motion may be steady or un- 
steady and the velocity of flight may be greater or less 
than the velocity of sound. 


MINIMUM DRAG DUE TO THICKNESS 


In the linearized theory the effects of camber, or lift, 
and thickness can be treated independently and later 
superimposed in any desired linear combination. Wing 
problems can therefore be divided into two classes— 
i.e., those involving thickness but no lift and those in- 
volving a distribution of lift over a surface having zero 
thickness. 

Consider first the case of a thin nonlifting body or 
wing, such as shown in Fig. 1. The surface is assumed 
to be everywhere nearly parallel to a horizontal plane 
in the direction of motion. A distribution of thickness, 
symmetrical above and below this plane, is supposed 
to be given over the plan form. Computation of the 
drag of such a body will involve, first, the determina- 
tion of the pressure distribution over its surface. The 
drag may then be obtained by one of two methods. In 
the first method the body is divided into elements of 
area, and the element of drag is obtained from the in- 
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Fic. 1. Thin airfoil with elements representing slope of surface. 
clination of the local pressure force. In the second 
method, which is related to the first through an integra- 
tion by parts, the body is divided into small elements of 
volume, and the drag is computed from the force of 
buoyancy on each element arising from the local gradi- 
ent of the pressure distribution. 


As is well known, the pressure distribution over such 
a thin, nearly flat body can be calculated by consider- 
ing only the inclinations or slopes of the surface ele- 
ments—without regard for the actual displacement of 
the surface away from the mean plane. The body can 
then be replaced by a large number of small elementary 
areas having slopes equal to the local slope of the sur- 
face. Each such element can be conceived to create a 
disturbance field of its own, and the geometric charac- 
ter of this field will depend on the properties of the 
medium (its elasticity, density, etc.) and on the cir- 
cumstances of the motion. The disturbance fields of 
all elements will be alike except for a factor of strength, 
which is measured by the product of the area of the 
element and its slope. The total disturbance field of 
the body may be obtained by superimposing the in- 
dividual disturbance fields of its elements. 


With the aid of this concept of superposition it may 
be shown that the drag of such a body is unchanged 
by a reversal of the direction of motion. The pressure 
at each point of the body is obtained by summing up 
the influence at this point of every element of the sur- 
face. The drag is then obtained by a second summation 
involving the total pressure force on each element and 
the slope of its surface in the dragwise direction. Thus, 
two integrations over the surface are required— or four 
integrals in all. A typical single element of this sum 
(see Fig. 1) will be represented by the pressure force 
on element 2 caused by the pressure field of element 1 
and multiplied by the slope of element 2. Now the 
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pressure field of any element, regardless of the distriby. 
tion or zone of action of the pressure, will be propor. 
tional in magnitude to the slope and area of the ee. 
ment. The magnitude of the small contribution to the 
pressure on element 2 caused: by element 1 will there. 
fore be given by the product m,Px, where my, is the 
“strength”’ of element 1 as determined from the prodyet 
of its slope by its elementary area (i.e., mm is equal to 
the frontal area of element 1) and where Pz is the pres. 
sure arising at element 2 from a unit disturbance at 
element 1. The increment of drag produced will be 
equal to the product of the pressure ,P2 by the are, 
and the slope of element 2—i.e., P22. In reversej 
motion, element | will lie in the same relation to element 
2 formerly occupied by 2 in relation to 1. Hence, P,, 
is equal to Py, provided the fluid itself is homogeneous, 
The signs of m, and m2 are both changed, but their 
product retains the same sign. The corresponding ele. 
ment of drag in reversed motion is, therefore, m2P,ym, 
which is exactly equal to the element of drag in forward 
motion. Such an equality can be assigned to every pair 
of elements and therefore must apply to the total drags 
on forward and reversed motion. 

It should be noted that the reversal principle involves 
hardly any restrictive assumptions aside from those in- 
volved in the linearization. It applies, therefore, to 
steady or accelerated motions and to subsonic or super- 
sonic speeds. 

The fact that the drag of the body is the same for 
either direction of motion leads to the consideration of 
a method of computing the drag wherein the forward 
and reversed disturbance fields are considered simul- 
taneously. It is found that the sum of the two drags, 
or twice the actual drag, can be obtained by super- 
imposing the two fields of disturbance velocities ani 
considering the pressures that would arise from motion 
of the combined field in the flight direction. The 
superposition of the two disturbance fields produces a 
fore-and-aft symmetry in the disturbance field of each 
element, with the result that the combined flow field is 
generally simpler in structure than the actual physical 
field and the calculation of drag is simplified. The 
pressure in the combined field is simply the difference 
between the pressures in forward and reversed mo 
tions. 

It can be shown that in the combined flow field the 
mutual interference drags of two distributions of thick 
ness are equal. Consider two distributions A and B as 
illustrated in Fig. 2. The two distributions can be di 
vided into an equal number of elements, and elements 
from the two distributions may then be paired in ai 
arbitrary way. Considering two elements such as ! 
and 2 shown in Fig. 2, it is evident from the symmetry 
of the combined field of each element that the drag d 
element 1 caused by element 2 is equal to the drag d 
element 2 caused by element 1. Since this equality 
holds for every pair of elements, it must hold for the 
total distributions A and B. 








AIF 


We 
for m 
Ness ¢ 
Let t 
that 
certai 
ness 
drag? 
specif 
drag | 
a sm 
produ 
tion i 
hence 
equal 
Sa fe 
ment 





distriby. 
€ propor. 
f the ele. 
On to the 
vill there. 
M, 1S the 
2 product 
equal to 
the pres- 
bance at 
1 will be 
the are, 
reversed 
’ element 
ence, Py 
geneous, 
ut their 
ding ele. 
m2P ym, 
forward 
very pair 
tal drags 


involves 
those in- 
efore, to 
or super- 


same for 
ration of 
forward 
d simul- 
‘oO drags, 
y super: 
ties and 
1 motion 
1. The 
yduces a 
of each 
v field is 
physical 
1. Th 
fference 
ed mo- 


ield the 
f thick- 
nd B as 
1 be di- 
lements 
1 in al 
‘h as | 
nmetry 
drag a 
drag ol 
quality 
for the 


MINIMUM DRAG 








es % 
/ ‘ 
fic. 2. Elements showing the mutual interference of two wings 
in combined flow. 




















AIRFOIL SHOWING ADDITION OF EQUAL OPPOSITE 
ELEMENTS OF VOLUME. 
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PRESSURE DISTRIBUTION FOR WHICH DRAG 
OF ADDED ELEMENTS IS ZERO. 


Fic. 3a (top). Fic. 3b (bottom). 


We are now in a position to determine the conditions 
for minimum drag with various specifications of thick- 
ness or volume. Consider, first, the following problem: 
let the plan form of the wing be given, and suppose 
that the thickness distribution is required to contain a 
certain specified volume. What distribution of thick- 
ness over the plan form will result in the minimum 
drag? First, let a distribution of thickness having the 
specified volume and designed to achieve the minimum 
drag be given. If the drag is actually a minimum, then 
asmall variation in the shape of the distribution will 
produce no first-order change in the drag. The varia- 
tion in shape must not change the original volume and, 
hence, may be divided into pairs of elements having 
equal and opposite volumes. The drag of each element 
isa force of buoyancy equal to the volume of the ele- 
ment multiplied by the local gradient of the pressure. 





Ni 


“Ni 


OF THIN WINGS 
Consider two such elements of volume added to the 
original distribution as shown in Fig. 3. Following the 
principle of superposition, the drag added by the ele- 
ments will be composed of three parts: (1) the drag of 
the elements in their own pressure fields; (2) the drag 
added to the original distribution by the pressure fields 
and (3) the drag of the elements in 
Since 


of the elements; 
the pressure field of the original distribution. 
the variation in thickness distribution is small compared 
to the original thickness, item (1) will be of smaller 
order than (2) or (3).* Furthermore, if use is made of 
the combined flow field, items (2) and (3) will be equal 
because of the mutual drag theorem previously demon- 
strated. The added drag is then equal to twice the 
buoyancy of the volume elements in the pressure field 
of the original distribution. Denoting the pressure in 
the combined field by P, and the volumes of the ele- 
ments by » and m2, as in Fig. 3, we have the following 
two equations: 


Vv; oe ve = 0 


(1) For no change in total volume: 
=) 4 (=) 7 
Ox 1 ™ Ox 9 - 


The two equations will be consistent if 


(2) For no change in drag: n( 


(OP, Ox); = (OP, Ox) 


Since such equations must hold for all positions of 
the elements within the specified plan form, we con- 
clude that the pressure gradient 0P,./0x must have the 
same value at all points of the wing. For minimum 

* The fact that the second-order variation in drag is always 
positive ensures that a stationary value will be a minimum and 
not a maximum. 

















PLAN 
VIEW 
A— 
Tl A-A 
SHOWING PRESSURE 
DISTRIBUTION 
Va <—— | es 


7Fc_ _consTANT 
dx 








Fic. 4. Condition for minimum drag with a given total volume 
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drag, therefore, the thickness must be distributed in 
such a way that the drag per unit volume is constant 
over the entire wing in the combined flow field. The 
significance of this condition in terms of the pressure 
distributions on an oblique airfoil in supersonic motion 
is illustrated in Fig. 4. 

Examples of bodies satisfying this condition are not 
difficult to find. One example is that of a thin, flat- 
tened ellipsoid accelerating in an incompressible fluid. 
The drag force in this case is in opposition to the accel- 
eration of the motion and, hence, is attributed to a 
“virtual additional mass’’ of the body. The surface 
pressures due to the acceleration are proportional to 
Oy/Ot, where ¢ is the surface potential. As is well 
known, the surface potential of any ellipsoid moving 
in a direction x parallel to a principal axis is of the form 
gy = KVx, where A is a constant depending on the 
proportions of the ellipsoid and V is the velocity. In 
accelerated motion, O0g/0Of and, hence, the pressure are 
also proportional to x, and this is true for motion in 
either direction except that the pressures are changed 
in sign. The pressure distribution in the combined 
flow, therefore, has a constant gradient in the direction 
of x, and the ellipsoid satisfies the necessary condition 
for minimum drag with a given volume. In this case 
the condition is perhaps better stated as the condition 
for minimum virtual volume with a given actual vol- 
ume. 

Another example is that of a biconvex airfoil of in- 
finite aspect ratio in steady motion at supersonic speeds. 
If the upper and lower surfaces are parabolic arcs, the 
Ackeret theory gives a straight-line distribution of pres- 
sure from nose to trailing edge for motion in either 
direction. 

The foregoing method may be readily extended to 
other problems of minimum drag involving different 
specifications on the thickness distribution. Perhaps 
the simplest of these is the case in which the thickness 
of the wing is specified along some line such as A—B 
in Fig. 5. Over the remainder of the wing plan form 
the thickness is to be distributed so as to achieve the 
minimum drag under the specified conditions. 

In this example we wish to consider variations in the 
thickness distribution which do not alter the shape of 
the cross section A-B. Any small element of slope 
having the strength m, as in Fig. 5, must therefore be 
followed by an equal opposite element m2 at some down- 
stream position ahead of the line A—B so that the thick- 
ness along the line A—B is not changed. The two equa- 
tions determining the necessary condition on the pres- 
sure distribution are: 


(1) For no change in thickness along A—B: 

m, + m. = 0 
(2) For no change in drag: mPa + mP2 = 0 
and these lead to the relation 


Pa = Pa 
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thickness. 


Hence, the pressure must not vary in the streamwis 
direction at points ahead of the section A—B. If the 
airfoil is further required to close along the trailing 
edge, a similar condition will be imposed on the pres 
sures in the region behind the line A—B (see Fig. 6). 
It will be noted that the condition for minimum drag 
does not impose a restriction on the spanwise variati0l 
of the pressure and that this variation will presumabl} 
be determined by the given shape of the section A-B 
If the shape of the section A—B is not given exactly 
but is merely required to have a certain frontal ares 
then it is found that the pressure in the combined flor 
field must have the same constant value at all points 


ahead of this section. The existence of a constant pres 
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MINIMUM DRAG 


sure means that the drag per unit of frontal area is con- 
stant. 

It is interesting to note that the condition for mini- 
mum drag with a given frontal area requires that the 
drag per unit frontal area be constant over the entire 
wing, while the condition for minimum drag with a 
given volume requires a constant drag per unit volume. 
If the maximum thickness of each streamwise section 
is specified, then the drag per unit thickness must be 
constant over each section. 


MINIMUM DRAG FOR A GIVEN LIFT 


The foregoing principles of calculation can be carried 
over without essential modification to the problem of 
The drag arising from the lift can be 
hence, the 


lift distribution. 
calculated independently of the thickness; 
wing can be considered a “‘lifting surface’’ slightly cam- 
bered or twisted in such a way as to support the de- 
sired distribution of lift. A given distribution of lift 
will induce a certain downwash velocity, and the result- 
ant inclination of the stream lines must agree with the 
slope of the lifting surface at all points. The drag of 
the lifting surface will be computed by integrating the 
product of the local lift by the local inclination of the 
lifting surface. 

The fact that a given distribution of lift has the same 
drag for either direction of motion can be seen by con- 
sidering the mutual drags of a pair of lifting elements. 
Each element of lift produces a distribution of down- 
wash over the plane. The actual configuration of the 
downwash field does not need to be specified, but its 
value at all points is assumed to be proportional to the 
lift of the element. Hence the downwash of an element 
|,at the position of another element /, may be expressed 
in the form w2, = 1, W2,, where W2; is an influence func- 
tion that depends on the geometric position of element 
2 in relation to element 1. The drag of element 2 
caused by element 1 is, therefore, /; W2/2. 

In reversed motion the roles of the elements are re- 
versed and the drag of element 1 caused by element 2 is 
now 1, Wil; and Wi, = We. Similarly, the drag of ele- 
ment 2 caused by element | in reversed motion is equal 
to the drag of /; caused by /, in forward motion. These 
statements do not exclude the possibility that each ele- 
ment lies wholly or partially outside the zone of in- 
fluence of the other element, and, hence, the rever- 
sal principle holds for unsteady or supersonic mo- 
tions. 

Suppose a plan form and a distribution of lift are 
given, and consider again the field of perturbation veloc- 
ities obtained by superimposing the fields in forward 
motion and in reversed motion. The horizontal pertur- 
bation velocities that give rise to the lift will be ex- 
actly canceled at each point of the surface. However, 
the downwash velocities over the surface will not, in 
general, be canceled, and the drag of the given distri- 
bution of lift for either direction of motion can be com- 
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puted by considering this lift to act in-the combined 
downwash field. 

If consideration is restricted to the method of calcu- 
lating drag by combining the flow fields in forward and 
reversed motion, it is easily seen that the mutual inter- 
ference drags of any two lifting areas are equal. The 
drag of a lifting surface A caused by the downwash of 
lifting surface B is exactly equal to the drag of B caused 
by the downwash field of A. This theorem follows 
from the symmetry of the combined flow field of any 
element and the equality of interference drags of two 
elements exactly as was previously demonstrated in 
the case of thickness. 

Consider now a lifting surface supporting a given lift 
L, and suppose that a distribution of lift over the plan 
form designed to achieve the minimum drag under the 
stated conditions is given. Any small variation in the 
shape of the lift distribution which meets the condition 
of constant total lift can be divided into pairs of ele- 
ments of equal, opposite lifts. If the drag added by 
every such pair of elements is zero, then the drag 
added by any small continuous distribution of the re- 
quired type will also be zero. 

Considering two such elements of a small additional 
distribution of lift, it is seen that the drag added to the 
original distribution will again consist of three parts: 
(1) the drag of the elements alone; (2) the drag added 
to the original distribution by the downwash fields of 
the elements; and (3) the drag of the elements in the 
downwash field of the original distribution. The drag 
of the elements alone will be of second order but will 
always be positive so that the second-order variation 
in drag is positive. Items (2) and (3) will be equal by 
virtue of the equality of mutual interference drags. 
For the first-order variation in drag, therefore, we need 
only consider the elements of lift acting in the down- 
wash of the original distribution. Since the elements of 
lift are equal and opposite, their drags will be equal and 
opposite, and the added drag will be zero if the down- 
wash is the same at each element. For all such pairs 
to contribute zero drag, the downwash must be constant 
over the whole wing. Hence, for minimum drag the 
lift must be distributed over the wing in such a way 
that the downwash in the combined flow field is con- 
stant at all points of the plan form. The condition is 
illustrated in Fig. 7. 


EXAMPLES OF THE MINIMUM DRAG OF LIFTING 
SURFACES 


The stated condition for minimum drag with a given 
total lift obviously includes the well-known condition 
for minimum induced drag in steady flight at subsonic 
speeds.‘ In that case the combined flow field is nothing 
more or less than the two-dimensional field of motion 
induced by the vortex wake that is normally trailing 
but which now extends to infinity both ahead of and 
behind the wing. The elliptic span loading, which 
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Fic. 8. Examples of minimum drag for a given lift. 


yields the minimum drag in this case, produces a uni- 
form downwash not only over the plan form of the 
wing but over the whole infinite vortex ribbon as well. 

In addition to providing the conditions for minimum 
drag, the combined flow field also provides a simplifica- 
tion of the actual computation of the drag in specific 
cases. Thus, in the case of steady motion at subsonic 
speeds, the superposition of the forward and reversed 
disturbance fields, which are three-dimensional, results 
in a two-dimensional field of motion. In supersonic 
motion, the combined field retains its three-dimensional 
character, but the integration for the downwash is 
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nevertheless considerably simplified when the forwarg 
and reversed fields are considered simultaneously, 

In the supersonic case, the computation of downwash 
over the surface involves the superposition of lifting 
elements whose individual disturbance fields are repre. 
sented by the potential of an elementary horseshoe 
vortex that has the form? ° 


I dy XZ 


9 » 9 (2 9 9 
TT (V+ 2)Vx2 —-— vy -— 2? 


g= 


for an element of lift [pV dy located at the origin. The 
downwash at the position x2, ye, 2 near the plane z = () 
arising from an element of lift at the position x, y 
is obtained by differentiating g with respect to z near 
z = 0 and introducing the displaced origin 


lr; dy, 
Wa = 
2a 


ixe = x1) 
[(ye = 41)” + 22°|V —~ = %5)* — (Ve baad 1)? = 2 


Normally, the computation of downwash over the 
lifting surface requires that the range of integration be 
limited to a portion of the plan form so as to exclude 
the forward branches of the Mach cones, where there is 
no real disturbance. In the combined flow field, how- 
ever, both branches of the cones are to be included, and 
the integration for the downwash at each point can be 
extended over the entire plan form provided the sign 
of the radical in the equation is chosen so that the down- 
wash of the lifting element has the same sign in both 
branches of the cone. 

One example that lends itself to calculation is that of 
an extremely long narrow wing or “‘lifting line’’ at an 
angle of yaw such that the line lies inside the Mach cone 
originating from its forward tip (see Fig. 8). Here the 
integration for the downwash in the combined flow 
field extends from one end of the line to the other, as in 
the subsonic case, and it is found that an elliptic dis 
tribution of lift results in a uniform downwash, yielding 
the value 


D = € i m2) (L? 1qb*) 


for the minimum drag. Here, the quantity m is the 
ratio of the slope of the lifting line to the slope of a 
Mach line and is less than 1.0 when the crosswise com: 
ponent of velocity of the lifting line is subsonic. As 
the angle of yaw is decreased so that the lifting line ap- 
proaches the Mach cone, the drag of this distribution of 
lift approaches infinity, a fact that must be taken to 
mean that the concentration of a finite amount of lift 
within an extremely narrow chordwise dimension 1s 
not permissible when the crosswise velocity component 
is near sonic velocity or is supersonic. 

For a lifting surface of narrow proportions lying neaf 
the center of the Mach cone, it is found that the expres 
sion for the drag can be separated into two components, 
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one of which depends on the lengthwise distribution of 
lift while the other depends on the spanwise loading. 
The minimum value in this case is found to be 

Pe. 


mgc* 


L? 
D=- 


] 
M? — 1 
xqb t 2 | ) 


where } is the overall span of the wing and c is the over- 
all length. The minimum value of the drag is achieved 
when both the spanwise and the lengthwise loadings are 


elliptical. 
DETERMINATION OF OPTIMUM PLAN FORM 


In practice, the determination of the best distribution 
of lift or thickness for a given plan form may not be so 
important as the problem of determining the best shape 
for the plan form. Unfortunately, it is difficult to 
formulate the latter problem in such a way as to lead to 
a single definite solution that will be generally useful. 
The foregoing analysis does, however, provide one im- 
portant principle concerning wings of minimum drag. 
Such wings will be characterized by the fact that the 
drag is stationary not only with respect to variations 
in the lift or thickness distribution but also with re- 
spect to variations in the shape of the outline. Sup- 
pose the plan form of such a wing, together with the 
distribution of lift or thickness satisfying the desired 
condition for minimum drag, is given. If the drag is to 
remain unchanged when the outline of the wing is dis- 
placed through a small distance, then the distribution 
of pressure or downwash corresponding to a minimum 
drag not only must appear over the plan form itself 
but must also extend for a small distance, without first- 
order variation, away from the edge of the plan form. 

That the foregoing condition leads to a stationary 
value of the drag can be seen by considering the effect 
of a small displacement of the outline of a lifting wing, 
as illustrated in Fig. 9. Supposing the total lifts of the 
original and the distorted wings to be the same, lift 
must be removed from the original area and placed on 
the added area. The first-order variation in drag can 
again be computed by considering only the effect of the 
original downwash field on the added distribution of 
lift. If this downwash, which was required to have a 
constant value at all points of the original plan form, 
remains constant in the regions of the added area, then 
a redistribution of lift from the original to the added 
area will not affect the total drag. 

Curves of variation of the drag with plan-form coef- 
ficients seldom show local maxima or minima but show, 


more often, a monotonic character. As an example 
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Fic. 9. Condition for drag to be a minimum with respect to 


variations in plan form. 


we may consider the induced drag of wings at subsonic 
speed. Here, the combined flow field, as previously 
noted, is simply the two-dimensional motion induced 
by the vortex wake, and, in the case of minimum drag, 
the downwash is constant over the whole wake ex- 
tending ahead of and behind the wing. Since the condi- 
tion for minimum drag persists along the whole wake, 
the drag is unaltered by any displacement of the wing 
boundary in the chordwise direction parallel to the 
wake. If one proceeds beyond the tip of the wing, 
however, the downwash changes sign abruptly, and 
large values of upwash are encountered. The drag is 
thus not a minimum with respect to changes in the span 
of the wing but can be continually diminished by re- 
moving area from the center portions of the wing and 
placing it at the tips. At supersonic speeds, the mini- 
mum drag for a given area and lift occurs when the 
wing surface is disposed along narrow lines lying near 
the center of the Mach cone and having the greatest 
possible length and span. 
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Design and Performance of a Simple 
Interferometer for Wind-Tunnel 
Measurements 


HARRY I. ASHKENAS* ano ARTHUR E. BRYSON?t 
California Institute of Technology 


ABSTRACT 

The design of an interferometer with a 3'/.-in. diameter light 
field and its application to a continuously operating 4- by 10-in. 
supersonic wind tunnel are discussed. The main features of this 
interferometer are (1) its extremely low cost compared to previous 
designs; (2) its light weight and flexibility; (3) its simplicity, 
permitting the instrument to be used also as a schlieren system; 
and (4) its unique feature of passing both light beams through the 
test section to minimize the effect of the side-wall boundary layers 
in the interferograms. 


INTRODUCTION 


Db pee THE PAST FEW YEARS, the Mach-Zehnder in- 
terferometer has come into increased use through- 
out this country for investigating high-speed aerody- 
namic phenomena. Ladenburg and his coworkers at 
Princeton University have been particularly active in 
this field.'~’ The main attraction of the instrument for 
the aerodynamicist is its ability to provide quantitative 
data on the density distribution for the complete flow 
field under investigation without the necessity of in- 
serting probes into the flow (which introduce changes in 
the flow and require many measurements to describe 
the complete flow field). The bibliography at the end 
of this paper lists references giving ample information 
on the optical theory of the instrument and describing 
its application to the study of two-dimensional and 
axially symmetric flow fields. 

In the course of work on boundary-layer-shock-wave 
interaction at GALCIT,f it became apparent that an 
interferometer could be 
However, the location and physical setup of the GAL- 
CIT 4- by 10-in. Transonic Wind Tunnel are highly un- 
favorable for the use of an interferometer according to 


used to great advantage. 


standard concepts: The tunnel is continuously operat- 
ing, the test section being directly above two 300-hp. 
centrifugal blowers; the power plant of the GALCIT 
10-ft. Wind Tunnel is close to the test section room; 
and the floor of the test section room is of temporary 
Thus, the vibration and noise 
The 


wood construction. 
levels near the test section are unusually high. 
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main problem, therefore, in the design of the interferom- 
eter was to obtain proper vibration isolation. Fur. 
thermore, it was required that the interferometer be 
mounted and removed easily to make room for other 
setups (such as schlieren) and also to make possible easy 
(For details of the wind 
The tunnel has a 


changes to the test section. 
tunnel see references 43 and 44.) 
flexible nozzle and variable second throat permitting 
operation from MM =0.7 to 1.5. 

The instrument developed is comparatively simple 
and has proved to be adequate. It is believed that a de- 
scription of the design and installation of the interferom- 
eter will be useful, since similar problems will probably 
arise with other wind tunnels. 

The instrument was designed and built during the 
course of an investigation that was begun under Air 
Force sponsorship but which, for the past year, has been 
carried out under the sponsorship of the National Ad- 
visory Committee for Aeronautics. 

The authors wish to thank Dr. Hans-Wolfgang Liep- 
mann, who started them on the design and whose idea 
it was to put both light beams through the test section. 
They also wish to thank Satish Dhawan, who sug- 
gested and designed the schlieren attachment to the in- 
strument. 


DIFFERENCES FROM PREVIOUS INTERFEROMETERS 


The principal innovations of this interferometer are 
believed to be (1) its low cost, (2) its light weight and 
simplicity of construction, (3) its simple attachment for 
taking schlieren pictures, and (4) its feature of passing 
both light beams through the test section. 

Table 1 contains an itemized cost estimate for the in- 
terferometer. It should be noted that the cost of re- 
producing the instrument now would be considerably 
less, mainly because of the elimination of the engineer- 
ing time. 

A rather thick boundary layer exists on the test- 
section windows so that, if only one beam traversed 
the test section, a large error would be introduced be- 
cause of the fact that the interferometer integrates the 
value of density through the flow. Most of the previous 
interferometers put one beam through the flow being 
investigated and the other beam through a glass- 
paneled compensating chamber in which the air density 
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can be adjusted to approximately the free-stream den- 
sity of the flow. Inthisinterferometer, both beams trav- 
erse the test section, one passing through the flow being 
investigated and the other passing through the test sec- 
tion 7*/; in. ahead of that beam. The boundary layer 
does not change its thickness appreciably in the distance 
between the two beams, so that both beams traverse 
effectively the same boundary layers and the latters’ 
effects are thus canceled. At the same time, the flow in 
the front beam corresponds to an ‘“‘automatic com- 
pensating chamber’’ providing free-stream density as 
the basis for the fringe shifts. This is particularly ap- 
plicable to supersonic flow where the flow ahead of the 
disturbance is uniform, but it can also be used for sub- 
sonic flows where the effect felt ahead of the model is 


small. 
OVERALL DESCRIPTION OF THE INTERFEROMETER 


(1) Optical Path 


A schematic diagram of the optical path through the 
instrument is shown in Fig. 1. Light from a mercury 
vapor lamp is passed through a filter and focused by a 
lens on a small mirror that is at 45° to the incident 
beam; this small mirror is also at the focus of a para- 
bolic mirror causing the reflected light to be parallel. 
This parallel beam is passed back over the small mirror 
(which shows up as a slight shadow on the interfero- 
grams but does not interfere at all in their evaluations) 
and is then split at the half-silvered mirror No. 1 into 
two beams; one beam passes through mirror No. 1, re- 
flects off the fully silvered mirror No. 3, then passes 
through the test section ahead of the model; the other 
beam reflects from mirror No. 1, passes through the test 
section over the model, then reflects from the fully 
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TABLE | 
Cost of the GALCIT Interferometer 


Piane-parallel plates $2,000.00 


Mirror frames 75.00 
Turning motors 50.00 
Light source 50.00 
Filter 30.00 
Lens 18.00 
Spherical mirror 27.50 
Paraboloidal mirror 100.00 
Shutter 40.00 
Camera back 18.00 
Suspension frame 34.00 
Optical bench 25.00 
Two test-section windows 13.00 


Traversing table, switch box, mirror cells, 


rails, suspension clamps, etc.—shop labor 1,500.00 


Temporary setup 40.00 
Engineering time 2,000 .CO 
Total $6,020. 50 


silvered mirror No. 2. The two beams are recombined 
at the half-silvered mirror No. 4 and are brought to a 
focus by a spherical mirror on a ground-glass screen. 
The plane mirrors are at the corners of a 7*/4- by 35%/4- 
in. rectangle, the plane of the rectangle being horizontal. 
The light beams are circular and can be varied in diam- 
eter from 2 in. to 3'/2 in. by varying the position of the 
lens. Both beams pass through the 12-in. diameter 
test section windows that are made of good quality 
plate glass of 1-in. thickness. (The use of plate glass in- 
stead of high-quality plane-parallel optical glass pro- 
vides a substantial reduction in cost; it was found that 
there was no need for the latter.) The instrument can 
be traversed forward and backward, as well as up and 
down, so that the small field of view is not a disad- 
vantage. Actually, most of the phenomena studied are 
local flow phenomena that are adequately covered by 
the 3'/s-in. diameter light field. Other test-section- 
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Schematic diagram of optical path. 








BA JOURNAL OF THE 


window arrangements for studying various types of 
flows are being designed. 


(2) Light Source 


The interferometer has been designed with a mini- 
mum of air-glass interfaces to prevent loss of light, there 
being only four such surfaces—two at the lens, and two 
at the filter. These surfaces are coated to prevent fur- 
ther loss of light. Asa result of this design, a 100-watt 
mercury light can be used, and, in conjunction with the 
filter, it gives sufficient monochromatic light for pictures 
to be taken at 1/800 sec. (Figs. 5 to 8 were taken with 
this light source at exposures of 1/400 or 1/800 sec.) 


(3) Mounting 


The four plane mirrors, the light source, the para- 
boloidal mirror, and the spherical mirror are mounted 
on a wooden box that is suspended from a welded-pipe 
framework by four sets of eight rubber bands each. 
The box with all the equipment mounted on it weighs 
only 250 Ibs., which is in contrast to all other interferom- 
eters that have come to the authors’ attention where 
the mirror base plates have been massive to provide 
rigidity. Here, the box is allowed to deflect in any man- 
ner it chooses, and, once it has assumed its deflected 
position and the mirrors have been adjusted, the shock 
mounting of the system effectively prevents any fur- 
ther misalignment from taking place. The action of the 
rubber bands rubbing over each other provides suf- 
ficient damping for the system. The natural frequency 
of the box is about | cycle per sec. This was found to 
provide not quite enough vibration isolation, and, 
therefore, the frame has also been isolated by placing 12 
in. of foam rubber under each of the four corners. This 
has provided sufficient vibration isolation to permit one 
to see the fringes while running the wind tunnel (two 
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Fic. 2. North half of the interferometer mounted on the 
GALCIT Transonic Wind Tunnel showing parabolic mirror, 
plane-parallel mirrors Nos. 1 and 3, and the vibration mounting. 
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300-hp. centrifugal compressors power the tunnel) and 
also while the GALCIT 10-ft. Wind Tunnel is running 
(its 750-hp. motor and fan are near the interferometer), 
The rubber bands are attached to jackscrews, permit. 
ting the box to be traversed in a vertical direction and 
allowing a fine vertical adjustment to be made in getting 
the light beam parallel to the surfaces of two-dimen. 
sional models. The whole interferometer can be tray. 
ersed along the test section on four casters; it can algo 
be rotated through small angles on these casters, also 
allowing a fine horizontal adjustment to be made in 
getting the light beam parallel to the surfaces of two. 
dimensional models. The interferometer can also be 
completely traversed away from the test section and the 
regular schlieren system traversed into its place. (Al 
though schlieren pictures can be taken with the inter. 
ferometer as described below, these pictures are not of 
such high quality as those taken with the regular 
schlieren system.) 

The camera box is mounted on the pipe frame so that 
putting film holders in the camera does not disturb the 
instrument’s adjustment. 


(4) Adjustments 


The various component parts of the interferometer 
are fitted with a large number of adjustments. This 
practice, while giving rise to a large number of possible 
misalignments, also results in considerable ease in set- 
ting up the instrument. Locking devices are present to 
prevent any change in the final setup. 


DETAILED DESCRIPTION OF COMPONENT PARTS 


(1) Monochromatic Light Source 


The monochromatic light source is a mercury vapor 
lamp, General Electric Type AH-4, housed in a Cenco 
No. 87269 lamp housing. The lamp house is attached 
to a brass tube that contains a Bausch & Lomb achro- 
matic coated lens of 142.9-mm. focal length and 46.0- 
mm. diameter. The brass tube also contains a recep- 
tacle for a Wratten No. 77A filter, which filters out the 
green line at 5,461A. in the mercury spectrum. At- 
tached to the top of the brass tube is a U-frame that 
supports a rectangular frame on which a small diameter 
wire is mounted carrying the 2- by 4-mm. first-surface 
plane mirror that acts as a secondary light source since 
it is at the focus of a 5'/-in. diameter, 15-in. focal length 
paraboloidal The 
mounted in an aluminum mirror cell of the type de- 
scribed by Strong.*! A central flexure post and four ad- 
justing screws provide the fine adjustment of position, 


mirror. paraboloidal mirror is 


while the course adjustment of position is provided by 
the sliding tubular mounting that is attached to the in- 
terferometer box. 


(2) White Light Source 


A 2- by 2-in. first-surface plane mirror has been 
mounted on a hinge in the mercury lamp housing and a 
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South half of the interferometer showing the plane- 
parallel mirrors Nos. 2 and 4, the spherical mirror, the shutter, 
the micrometer knife edge for schlieren pictures, and the camera. 


Fic. 3. 


hole cut out of the housing so that by turning a screw 
the mirror is brought to a position that is at 45° to the 
light from a 500-watt motion-picture projector lamp 
that is mounted adjacent to the lamp housing. This 
light source is the same distance from the mirror as the 
mercury light so that both sources are optically at the 
same point with respect to the lens (see Fig. 1). 


(3) Plane-Parallel Mirrors 


The plane-parallel mirrors are carried in a frame 
shown in Fig. 4. Each mirror itself is held in a split 
brass ring by means of a soft rubber, full-circumference 
support. Two gimbals allow rotation of the mirror 
about a horizontal and vertical axis. Rotation of the 
gimbals is obtained by means of remotely controlled, 
reversible electric clock motors. The output shaft of 
these motors (Haydon No. 3200) makes one revolution 
per 8 days, resulting in an extremely fine control of the 
angular displacements of the mirrors. The gearbox on 
the motor is equipped with stamped metal gears, and 
a not inconsiderable amount of backlash is present; 
hence, a simple torsion spring is placed at each bearing 
to supply a backlash-eliminating force. 

The frames of all four mirrors are identical. 
min castings are used for the outer frame and the inner 
gimbals. There is no provision for coarse angular ad- 
justment of the mirrors, which makes the initial adjust- 
ment somewhat tedious but thereafter no coarse ad- 
justment is needed. Mirror Nos. 1, 2, and 3 (see Figs. 2 
and 3) are bolted to the top of the plywood box. Mirror 
No. 4 is mounted on a traversing table to provide the 
adjustment for equalizing the optical path lengths of 
the two arms of the interferometer. The traversing 
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mechanism features a Le Blond lathe bed type way; a 
traversing screw of l-mm. pitch, driven by a hand- 
operated 50:1 worm-gear reduction unit; a simple nut; 
and a pin-driven gravity-emplaced traversing table. 


(4) Camera and Schlieren System 


The photographic components consist of a 3'/2-in. 
diameter, 24-in. focal length spherical mirror, Eastman 
Synchro-Rapid 800 shutter (speeds up to 1/800 sec.), 
and a 4- by 5-in. Graphic camera back with a ground- 
glass viewing screen. Because of the shutter’s small 
diameter, it has been placed at the end of a truncated 
conical piece extending from the camera box (see Fig. 
3). The focal point of the light beam reflected from the 
spherical mirror is just inside the shutter in the conical 
part of the camera, and a micrometer with a razor blade 
mounted on its end is placed there as the knife edge of a 
schlieren system. The knife edge can be traversed along 
the light beam to obtain best sensitivity, and it can be 
rotated to give ‘‘knife-edye horizontal’? and ‘‘knife- 
edge vertical’’ pictures. To take schlieren pictures, the 
forward light beam of the interferometer is blocked 
out, the filter is removed, and the white light source is 
used. 

The film used is Ansco Triple-S-Pan, which is suf- 
ficiently fast to take good pictures at 1/400 sec. Slightly 
better pictures are obtained using Eastman 103-G 
spectroscopic plates, which are especially sensitive in the 
green range containing the wave length of the mono- 
chromatic light used (5461A.). These plates are fast 
enough for pictures at 1/800 sec. and have the added ad- 
vantage that there is no film shrinkage that could cause 


errors in fringe shift measurement. (Actually, by com- 





Close-up of plane-parallel mirror frame showing gimbals 
and rotating motors. 
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paring the film and the glass plates, no significant er- 
rors in fringe shift measurement have been found.) 


PERFORMANCE OF THE INTERFEROMETER 


(1) Initial Adjustment 


The initial adjustment of the interferometer was 
made using the pentaprism method as described in 
reference 30. This adjustment required several days to 
perform because of the slow mirror movements, but, 
once completed, it does not have to be done again. 


(2) Daily Adjustment 


Because of slight deflections of the box caused by tem- 
perature and humidity changes, the interferometer gets 
slightly out of adjustment from day to day. How- 
ever, with an experienced operator, fringes of the de- 
sired spacing and orientation can be focused on the 
ground-glass screen in less than 10 min. using the 
“near-and-far cross-hair’’ method of adjustment de- 
scribed in reference 30. For this, an auxiliary plane 
mirror is placed at 45° to the light beam just ahead of 
the spherical mirror, and an autocollimating telescope 
is placed in the reflected beam from this mirror. If the 
telescope has been previously focused on a distant ob- 
ject, it will now focus on the two images of the secondary 
light source (the 2- by 4-mm. mirror) if the light is 
parallel. These images serve as the ‘far cross-hairs’’ 
being optically at infinity; the ‘‘near cross-hairs’’ are 
the shadows of this 2- by 4-mm. mirror, which is only a 
few inches from the half-silvered mirror No. 1. Now, 
focusing on the ‘‘far cross-hairs,’’ mirror No. | is rotated 
about its vertical and horizontal axes until the two 
images coincide; then, focusing on the “near cross- 
hairs,’’ mirror No. 4 is rotated about its horizontal and 
vertical axes until these two images coincide. This 
latter adjustment slightly disturbs the previous adjust- 
ment of mirror No. 1, which is then done over again, 
and, similarly then, mirror No. 4 is readjusted. This 
process is rapidly convergent, usually requiring only the 
first adjustment of each mirror. 


(3) Horizontal or Vertical Fringes 


The geometry of this particular interferometer is such 
that by now rotating both mirror No. 1 and mirror No. 4 
in the same direction about horizontal (or vertical) axes 
in a ratio of 7 to 5, horizontal (or vertical) fringes ap- 
pear very nearly in the plane of the ground-glass screen. 
The spherical mirror is rotated slightly, placing the re- 
flected beam outside the camera so that a white card 
can be placed in the beam to determine the plane in 
which the fringes are focused. By reference to the 
“linear interferometer representation’ described in ref- 
erence 34, the operator can then tell which way to 
rotate mirror Nos. | and 4 to focus the fringes precisely 
on the ground-glass screen with the desired spacing and 
orientation. Previously, the spherical mirror has been 
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located so as to focus a plane that is one-third of the 
way from the exit window to the entrance window on the 
ground-glass screen. This minimizes refraction errors 
in the interferogram.’ Thus, both model and fringes 
are in focus on the ground-glass screen. 


(4) Equalizing Optical Path Lengths 


The next adjustment is to equalize the optical path 
lengths in the two legs of the interferometer. This js 
done by placing an auxiliary mirror in front of the para- 
boloidal mirror so that it cuts out half of the monochro- 
matic light beam. An auxiliary white light source is 
then reflected onto this auxiliary mirror so that its light 
travels the same circuit as the monochromatic light. 
The image formed near the ground-glass screen is then 
half monochromatic light fringes and half white light 
field. Mirror No. 4 is then traversed until the white 
light fringes appear in the white light field; usually 
only a few turns of the traversing wheel are required, 
and often the white light fringes have not ‘‘strayed” 
since the last use of the interferometer. The correct 
direction to traverse is easily checked by the fact that 
monochromatic fringes should get more ‘‘contrasty” 
if the traversing is done in the correct direction. 


(5) Eliminating ‘‘Twist”’ 


The next adjustment is to eliminate ‘‘twist’’ from the 
beam—i.e., to place the two images directly above each 
other (for the case of horizontal fringes). This gives the 
maximum contrast in the fringes. This can be done by 
rotating mirror Nos. | and 4 about their vertical axes 
slightly (for the case of horizontal fringes) until the 
monochromatic fringes are parallel to the white light 
fringes and they are both horizontal. This adjustment 
works only if divergent white light is used in the auxil- 
iary beam. Getting rid of ‘‘twist’’ can also be done by 
trial and error without the white light by simply rotat- 
ing about the vertical axes (as above) until the fringes 
appear to have maximum contrast. 


(6) Infinite Fringe Adjustment 


The infinite fringe adjustment is obtained by first 
making the adjustment above for either vertical or 
horizontal fringes and then slowly increasing the fringe 
spacing keeping the fringes focused on the ground-glass 
screen (this is easily done by rotating mirror Nos. | and 
4 in the same direction in a ratio of 7 to 5). 


(7) Aligning Light Beam to Models 


The final adjustment is to align the light beam to the 
surface of two-dimensional models or, in the case of 
axially symmetric models, to align the light beam per- 
pendicular to the axis of symmetry. 


(8) Procedure During a Run 


Before the tunnel is turned on, an undisturbed fringe 
picture is taken. Then, after the flow has been turned 
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Undisturbed horizontal fringes. 
positive of Fig. 5 and a negative of Fig. 6. 
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on and has reached an equilibrium flow condition, the 
monochromatic fringe pictures are taken. While still 
running, the filter is removed, the mirror in the lamp 
housing is rotated into place, and the white light then 
becomes the light source. The white light fringe picture 
is then taken. Next, the front beam of the interferom- 
eter is blocked out and the schlieren system is adjusted; 
then, a schlieren picture is taken. A shadowgram can 
also be taken at this point, if desired, simply by moving 
the knife edge back out of the beam. All of these opera- 
tions can be performed in a few minutes so that a mini- 
mum of tunnel time is used for each run. 


(9) Evaluation of Interferograms 


Knowing the stagnation temperature and pressure 
in the settling chamber and the Mach Number in the 
test section, the density of the air in the front beam is 
readily calculated by the isentropic flow relations. If, 
in addition, we know the span of the model, the wave 
length of the light used, and the Gladstone-Dale con- 
stant for this wave length, the interferogram can then 
be completely evaluated to give the density distribution 
throughout the entire flow field about the model. 

For two-dimensional flows, the evaluation is par- 
ticularly simple, since the density change is directly 
proportional to the distance the fringes shift from their 
undisturbed positions. By comparing the fringe pic- 
tures taken with and without flow, this fringe shift can 
be determined. A particularly simple way of doing 
this, mentioned in reference -18, is to superimpose the 
negatives of these two pictures; faint lines then appear 
which are lines of constant fringe shift. Since the nega- 
tives are small, this cannot be done too accurately if the 
fringe spacing is narrow. To overcome this difficulty, 
enlarged transpareticies of the undisturbed fringe pic- 
tures are made. Placing this transparency on the en- 
larger table, the other picture is focused until the two 
pictures ‘‘line up’’—i.e., the fringes ahead of the shock 
wave fall on top of each other. Then, the printing paper 
is put under the transparency and a picture is printed, 
giving the lines of constant fringe shift as shown in Fig. 
7. To determine the proper number of fringe shifts as- 
sociated with each line, the white light picture must be 
used to trace the central fringe through the shock wave, 
or the fringe shift at a pressure tap on the model is cal- 
culated. 

The infinite fringe picture serves as a check on the 
above method, but its accuracy is dependent on the 
accuracy with which the mirrors were made; _there- 
fore, it is usually only accurate (for this interferometer) 
within +'/, fringe shifts. The method of finite fringes 
above “‘calibrates out’’ the inaccuracies of the optical 
system. 

The evaluation for axially symmetric flows is de- 
scribed in reference 1; it is considerably more difficult 
and tedious than the evaluation for two-dimensional 
flows, since the fringe shift is related to the density 
change by an Abel integral equation. 





Fic. 9. White light picture showing how fringes may be traced 


through a shock wave. 


(10), Obtaining Two-Dimensional Flows 


One difficulty in the use of the interferometer for two- 
dimensional flows is that two-dimensional flows are 
difficult, if not impossible, to obtain, and, consequently, 
errors are introduced into the interferogram. As pre- 
viously described, passing both beams through the test 
section is a help in overcoming this difficulty by mini- 
mizing boundary-layer effects. Permitting models to 
span the tunnel only up to the beginning of the side-wall 
boundary layers seems to minimize the three-dimen- 
sional effects, although this has not been thoroughly 
checked as yet. The three-dimensional or ‘“‘tip’’ effects 
probably vary considerably with the model being used 
and with the flow velocity. Also, the interaction of the 
shock waves with the side-wall boundary layers alters 
the tip effects and the boundary layers. The tip effects 
and the shock-wave-boundary-layer interaction effects 
cannot be separated and must be considered together in 
determining the accuracy of the interferograms. Pre- 
liminary tests on detached shock waves indicate that 
both boundary-layer and three-dimensional effects are 
practically eliminated by careful choice of the model 
span. 


CONCLUSION 


Figs. 5 to 10 are typical interferograms taken with 
this instrument. Fig. 5 is an interferogram showing 
horizontal undisturbed fringes taken through the test- 
section windows with the wind tunnel turned off. Fig. 
6 shows the flow field about a 10° semiangle wedge with 
a detached shock wave (the pegs on the flat portion fol- 
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Schlieren picture made with knife edge mounted on 
the interferometer. 


Fie. 10. 


lowing the wedge are for aid in alignment of the light 
beam to the leading edge of the wedge.) An investiga- 
tion of the flow field behind such detached shock waves 
isat present being carried on using the interferometer. 
Fig. 7 is an interferogram made by superimposing the 
negatives of an interferogram with flow and one without 
flow; the shadowy lines are the lines of m + '/2 fringe 
shifts (m an integer) and are lines of constant density. 
Fig. 8 is an interferogram taken with the ‘‘infinite 
fringe,’ and the dark lines are approximately lines of 
n+ 1/; fringe shifts (m an integer) and are lines of con- 
stant density. Fig. 9 is an interferogram taken with 
white light instead of monochromatic light and shows 
the manner in which the fringes can be traced through 
ashock wave. Fig. 10 is a schlieren picture taken with 
the special attachment to the interferometer described 
above. 

The performance of this interferometer on the 
GALCIT Transonic Wind Tunnel has proved satis- 
factory to date, and it is hoped that many interesting 
aerodynamic phenomena may be studied with it in the 
future. 
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First U. S. National Congress of Applied Mechanics | 


Plans have been formulated to hold the First U.S. National Congress of Applied Mechanics on 
June 11-16, 1951, at Chicago, Ill. The Illinois Institute of Technology will bz host to the Congress. 

Members of the Institute, one of the sponsoring societies, are cordially invited to submit papers, 
which should constitute original research in applied mechanics. (Their length must not exceed 5,000 | 
words or the equivalent in equations, tables, and diagrams.) Abstracts must be submitted to the Chair- | 
man of the Editorial Committee before April 14, 1951; complete papers must be submitted prior to June 
11, 1951. | 

Additional information about the Congress may be obtained by writing to N. M. Newmark, 
Secretary of the Congress, University of Illinois, Urbana, III. All members are invited to attend. 
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Remarks on Low-Aspect-Ratio 
Configurations in Supersonic Flow’ 


P. A. LAGERSTROM? anpb M. E. GRAHAM} 
Douglas Atrcratt Company, Ince. 


SUMMARY 


Simple low-aspect-ratio configurations are examined to ob- 
tain a qualitative understanding of wing-fuselage interference 
and “afterbody”’ effect of a missile. Slender-body theory gives 
for a wing-body combination a value of lift which does not de- 
pend upon afterbody length. This value of lift is compared with 
that given by the more exact linearized three-dimensional theory. 
Nonlinear and viscous effects are not considered. 

A theorem is proved concerning the lift of certain wings (or 
wing-fuselage combinations) with a long afterbody of span 
equal to the maximum span of the wing. It states that at any 
supersonic Mach Number the value of lift according to linear 
theory tends to that predicted by slender-body theory as the 
length of the afterbody tends to infinity. The effect of a finite 
afterbody can also be determined by linear theory, and a com- 
parison can be made with slender-body theory for certain planar 
systems. 

These results may be applied qualitatively to wing-body com- 
binations. They indicate that, for no afterbody, slender-body 
theory overestimates the lift of the combination and that, for a 
short afterbody, it may equal or even underestimate the lift. 
For an infinitely long afterbody, it agrees exactly with linear 
theory if the afterbody span equals the maximum wing span; if 
the afterbody is narrower, it overestimates the lift 


NOTATION 

a = radius or half-width of body 

A = aspect ratio 

A, = BA, reduced aspect ratio 

c = chord 

C, = pressure coefficient 

E = FE V1 — 9), complete elliptic integral of second 
kind 

k = a/s, ratio of body radius (or half width) to maxi- 
mum semispan of wing-body combination 

l = lift per unit span of wing 

L = lift 

M = free-stream Mach Number 

q = (1/.)pU? 

s = maximum semispan of wing-body combination 

od = s + x; tan w, local semispan of wing-body com- 
bination 

U2 = component of flow velocity in x2 direction 

U = free-stream velocity 

1, Xo, Xz; = Cartesian coordinates, directed in spanwise, verti- 
cal, and streamwise directions, respectively 
Origin is at intersection of body axis and wing 
trailing edge 

a = angle of attack (rad.) 


Received May 25, 1950. 

* This paper is nearly identical in content with the first part 
of reference 1. 
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Qa = apparent angle of attack (geometrical plus that 
induced by upwash field) 


= VM? -1 


B 
w = half apex angle of delta wing 
Q = Btanw 
Subscripts 
A = afterbody 
B = total body 
W = wing 


INTRODUCTION 


Ph pqeny A TYPICAL MISSILE such as shown in Fig. 1. 

It will be assumed that the fuselage has a circular 
cross section and is cylindrical except near the nose. 
The lifting surfaces are flat plates of some simple plan 
form like a delta wing; the wings are set at zero inci- 
dence to the body. A solution for this missile as a 
whole may be obtained formally by slender-body 
theory.*~* From this solution it would appear that 
the tail is completely ineffective; its geometrical angle 
of attack (assuming zero incidence to body) is com- 
pletely canceled by the upwash field induced by the 
wing and the fuselage. However, this is a highly un- 
realistic result obtained by using slender-body theory 
beyond its intended application. 

To analyze this situation, the missile will first be 
considered without the tail. Even in this case there is 
an unrealistic consequence of slender-body theory- 
namely, that the afterbody carries no lift and, hence, 
that its length is unimportant. This is connected 
with the fact that slender-body theory shows no 
Mach Number effect. Therefore, a critical discussion 
of the relation between slender-body theory and the 
more exact linearized theory is given. 

The question of the lift of the tail requires a study 
of the development of the vortex sheet. This problem 
will not be discussed here. However, once the upwash 
field in the tail region has been found by some method, 
the problem of wing-tail-fuselage interference reduces to 
a study of tail-fuselage interference. The tail is a 
lifting surface like the wing, although its apparent 
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angle of attack may vary spanwise in some complicated 
manner. In general, the tail-fuselage combination has 
no afterbody. Because of this, a different estimate 
may have to be used for tail-fuselage interference than 
for wing-fuselage interference. 


WING-FUSELAGE INTERFERENCE, AFTERBODY EFFECTS 


In treating the lift of a wing by linear theory, it is 
always possible to use a superposition principle. First, 
the upwash field of the fuselage alone is computed. 
Let the apparent angle of attack of the wing be defined 
as the geometrical angle of attack plus the angle of 
attack induced by this upwash field. The problem 
then reduces to the study of a wing with this apparent 
angle of attack in the presence of a fuselage at zero 
angle of attack. This reformulation of the problem 
is convenient for wings of high aspect ratio and with a 
high ratio of span to body diameter. In this case, the 
fuselage at zero angle of attack influences only a small 
part of the wing, and its influence may be understood 
and estimated with the aid of the idea of preservation 
of lift.° 

However, a different approach to the problem is 
advantageous for the cases considered in this paper, 
especially when the leading edges are subsonic. Since 
slender-body theory* ‘ gives directly a solution for the 
original combination of the wing and fuselage, it has 
been selected as the basis of the investigations. Cor- 
rections and comparisons to this theory are presented 
in the discussions. Slender-body theory does not show 
that the lift has a dependence on Mach Number or the 
length of the afterbody. A linearized theory would 
show these effects, but a complete treatment with the 
linearized theory would be cumbersome. Therefore, 
these effects will first be studied for planar systems 
(wings) for which the two methods may be directly 


compared. 
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(1) Low-Aspect-Ratio Rectangular Wing 


Linearized Theory.—Equations for the lift and center 
of pressure of a rectangular wing for a reduced aspect 
ratio varying between 1 and are easily computed— 
e.g., by the method of conical flow.® For lower aspect 
ratios, the method of cancellation of lift® 7 or Evvard’s 
Near the leading edge the pres. 
sure distribution while further 
downstream there is a complicated system of expansion 
and compression waves which eventually decreases the 
pressure distribution to a small fraction of its two-dimen. 
sional value. (This is discussed in detail in reference 7. 

The lift of a rectangular wing for which !/. < A, < o 
is shown in Fig. 2. It is interesting to note that, if the 
span is held constant, the lift has its maximum value 
for a reduced aspect ratio of one; for lower reduced 
aspect ratios the lift decreases. This means that, if an 
additional strip (ABCD of Fig. 3) is added to a wing 
of reduced aspect ratio smaller than 1, at a positive 
angle of attack, the lift carried by this strip must be 
negative. The addition of this strip thus has the ef. 
fect of decreasing the lift and forcing the center of 
pressure forward. 

Slender-Body Theory.- 
rectangular wing as a limiting case of a delta wing with 
a flat-plate afterbody (Fig. 4). (By ‘“‘afterbody”’ is 
méant the flat plate downstream of AB in Fig. 4.) Ac- 
cording to slender-body theory, all the lift is carried 
upstream of AB; the afterbody has zero perturbation 
pressure if it is at the same angle of attack as the 
forward delta wing. Furthermore, the lift is inde- 
pendent of the opening angle w; it depends only on the 
maximum span according to the formula 


method* may be used. 
is two-dimensional, 


It is convenient to treat the 


L = (m/2)qao(2s)? (1) 


The rectangular wing is then obtained as a limiting 
case by letting the angle w tend to 7/2. Since Eq. (1) 
is independent of a, it is still valid for the rectangular 
wing. Thus, for the rectangular wing all the lift would 
be concentrated at the leading edge, the rest of the 
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wing would have zero perturbation pressure, and the 
lift would be independent of the aspect ratio. This 
result is plotted in Fig. 2. 

Comparison of the Two Theories.—Fig. 2 indicates 
that, as the aspect ratio approaches zero, the lift, ac- 
cording to linearized theory, tends to the value pre- 
dicted by slender-body theory. It appears that 
slender-body theory anticipates the value of lift which 
would be obtained if the wing were of a sufficiently 
small aspect ratio. Note that, as the Mach Number 
M approaches 1, more and more of the reflections, in 
the system of compression and expansion waves which 
occurs according to linearized theory, are compressed 
toward the leading edge. According to slender-body 
theory (which corresponds to the limiting case M = 1), 
this whole wave system is, so to speak, concentrated at 
the leading edge. In the following section it will be 
shown that the asymptotic result indicated by Fig. 2 
is actually correct not only for the rectangular wing but 
for a large class of plan forms. 


(2) General Theorem About Lift of Zero-Aspect-Ratio 

Wings* 

Consider the wing shown in Fig. 5. The following 
assumptions are made: 

Plan Form.—No part of an edge is trailing. 
stream of a certain chordwise station (AB in Fig. 5), 
the side edges become parallel to the free stream. The 
wing extends infinitely far downstream; its aspect ratio 
is thus zero. 

Angle of Attack a.—a may be arbitrary everywhere ex- 
cept sufficiently far downstream where (for any fixed 


Down- 





* A special case of this theorem is stated (without proof) in 


reference 10. 


spanwise location x,) a becomes a constant, a(x,).t 
Theorem: For the wing described above at super- 
sonic speeds, linear theory gives the following results: 


(a) The lift is independent of Mach Number. ) 
(b) The angle-of-attack distribution on the 
forward part of the wing contributes 
nothing to the lift. The lift is deter- | 
mined exclusively by the values of ” 
ao(%). ( 2) 
(c) The value of the lift is identical with | 
the one that would be obtained by | 
slender-body theory (which corresponds | 
to M = 1). ) 


It can easily be seen that the proof given below 
will also be valid for appropriate nonplanar systems. 
Use of this will be made in Section (4). 

Proof of (a).—Consider a rectangular wing (Fig. 6a) 
of finite aspect ratio for which the angle of attack is 
zero, except along a streamwise strip of differential 
width where it is a. Let the Mach Number M equal 
Mo. With the aid of the Prandtl-Glauert transforma- 
tion, it is shown that wing I at 17 = My, has the same 
lift as the equivalent wing II at M = 7/2. (The 
equivalent wing has the dimensions and angle-of- 
attack distribution shown in Fig. 6b.) Now let the 
chord c grow to infinite length (Fig. 7); then wings I 
and II become identical. Hence, the lift of the wing 
of Fig. 7 is the same for M = Myas for M = +/2—i.e., 
it is the same for all supersonic Mach Numbers. Thus, 
(a) of the theorem is proved for a special type of zero- 
aspect-ratio wing. 

Next, consider a wing, like that in Fig. 8, which has 
a deflected square of infinitesimal area but which is 
otherwise at zero angle of attack. For this wing the 
lift is zero and, in particular, is independent of Mach 
Number. (This result follows if one wing of the type 
of Fig. 7 is subtracted from another such wing, since a 
supersonic wing has no upstream influence.) By 


t It would be sufficient to require that a(x, 0, x3) converges 
to a finite value ao(x;) as x3; approaches infinity and that the side 
edges become asymptotically parallel to the free stream. We 
shall not be concerned with such mathematical refinements. 
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application of the superposition principle to the above 
results, it follows immediately that parts (a) and (b) of 
the theorem are true for any wing of rectangular plan 
form which satisfies the hypotheses of the theorem. 

Finally, any wing of the type depicted in Fig. 5 may 
be considered as a rectangular wing by completing its 
plan form as in Fig. 9, giving to the added parts an 
angle-of-attack distribution that corresponds to the 
upwash of the original wing. 

This completes the proof of parts (a) and (b). Part 
(c) follows from the fact that, since the lift is the same 
for all Mach Numbers, its value may as well be com- 
puted for the particular Mach Number of 1. The 
value is then the one that slender-body theory gives. 


(3) Narrow Delta Wing with Flat Afterbody 


Full-Span Afterbody.—Consider first the wing and 
afterbody of Fig. 10a. As was shown in the previous 
section, when the afterbody becomes extremely long, 
the lift as given by linear theory is very nearly that pre- 
dicted by slender-body theory. For this particular 
wing this result will be studied in more detail. 

It is of interest to check the consistency of theorem 
(2) by the following argument: The lift of wing (a) 
(Fig. 10) equals that of (b) plus that of the afterbody 
(c) at its apparent angle of attack. The lift of (b) is 


Lw = rqao(2s)?/2E(VW1 — 92) (3) 
For an afterbody of infinite chord, only the apparent 
angle of attack at the trailing edge is significant. The 


induced part of the apparent angle of attack is directly 
proportional to the upwash in the Trefftz-plane. 
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a, — a = Ue U= — Ay E(V1 iad 0?) (4) 


a 


The lift of the afterbody is then 
i. = bal 9-\2 — bed 9s)2 (J = 
va 9a@q (2s)" = 5 qa (2s) — E (5) 


(For simplicity, the modulus of E, a complete elliptic 
integral of the second kind, is omitted.) Thus, adding 
L, to Ly gives the slender-body result 


Lwa = Lw + La = (/2)qao(2s)? (6) 


Since theorem (2) actually was used in the above argu- 
ment, it cannot be said that the theorem has been 
checked by this argument even for this particular wing, 
However, at least the consistency of various methods 
has been verified. 

Furthermore, for this plan form, as well as for the 
rectangular wing, it is possible to obtain linear theory 
results with the afterbody of finite length. Hence, 
for this case also, the validity of applying theorem 
(2) (i.e., slender-body theory) to a wing of low but non- 
zero aspect ratio may be checked. The linear theory 
results have been obtained from the charts of reference 
9. The lift is plotted in Fig. 2 as a function of the re- 
ciprocal aspect ratio of the afterbody for three wingswith 
different apex angles. These values are compared with 
the slender-body results. For no afterbody, slender- 
body theory predicts too much lift; as the afterbody 
first grows from zero to a finite length, it recovers more 
lift than slender-body theory predicts; finally, the lift 
appears to drop and, according to the previous section, 
would approach asymptotically the slender-body value 
if the curves were continued far enough. These curves 
are similar to that for the limiting case of a wing with 
a 90° half apex angle, the rectangular wing, which is 
shown on the same graph and was discussed previ- 
ously. 

Note that as w — 0 the pressure distribution as given 
by linearized theory approaches that given by slender- 
body theory. This is one aspect of slender-body 
theory. The aspect of this theory emphasized in this 
paper is that for any value of w (even for supersonic 
leading edges) the total lift, but not the pressure dis- 
tribution, approaches that predicted by slender-body 
theory as the length of the afterbody increases. In the 
first case the relevant parameter is actually 8 tan w; in 
the second case, c,/8(2s). 

Narrow-Span Afterbody.—For the applications it is 
desirable to know what happens when the span of the 
afterbody is smaller than the maximum span of the 
delta wing (Fig. 11). Theorem (2) does not apply 
directly to this configuration, but it may be applied 
separately to the afterbody as was done at the begin- 
ning of this section. Slender-body results are not real- 
ized as the chord of the afterbody tends to infinity. 
As it is to be expected, the lift is smaller than for the 
combination with a full-span afterbody (Fig. 10). It 
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infinitely long afterbody. 


may be expressed by the equations for that case [Eqs. 
(3), (5)] by introducing a factor dependent upon the 
ratio of body span to wing span. 


mgay(2s)? -rgay(2a)? l 5 
Lwa - oR + 2 L.— E -_ Lw + k*L4 
(7) 
where k = a/s and Ly, L, refer to configuration of 
Fig. 10. 


This plan form displays more of the characteristics 
of an actual wing-body combination than did the pre- 
vious examples, but only the lift for the limiting case 
of an infinitely long afterbody may readily be found.* 
This is plotted in Fig. 12 as a function of k. 


(4) Delta-Wing Fuselage Combination with Cylindrical 
Afterbody 


The wing-body combinations considered in the pre- 
vious section had flat bodies. These were discussed 


because they do contribute to the understanding of 


* Linearized calculations for the case of a short afterbody may 
be found by the method of cancellation of lift.7_ These cumber- 
some calculations have not been carried out. However, from 
the behavior of the corresponding curves for k = 1, it is to be 
expected that, for k values that are not too small, the full slender- 
body value of the lift may be realized for certain moderate values 
of the length of the afterbody. 


the wing-body problem and because solutions for such 
cases are more easily obtained than for a combination 
with a more realistic fuselage shape. Spreiter* and 
Ward‘ have both obtained the slender-body results for 
a combination with a circular-cylinder body (Fig. 13). 
Their results for the total lift of such a configuration are 
given in Eq. (8) and Fig. 12. 


Lwp = (x/2)qao(2s)? (1 — k? + k'4) (8) 


For this more realistic configuration, various results 
for the pressure distribution on the wing will also be 
examined. Spreiter* has found the theoretical pres- 
sure distribution for various slender wing-body com- 
binations. The pressure distribution on the wings of 


the combination of Fig. 13 is 


mune] 
W+O))-@Ob+O) 


The spanwise 


(Cp)w = 


s’ = s + x; tan w is the local semispan. 
lift is then 


0 
lw = —2q¢ f Cp dx3 = 4qaya X 
(x1 — s)/tanw 
(1 “s =\(5 a “) (10) 
s? a? 3° 
and the total lift carried on the exposed wing area is 
, ,|f/s , a’ 
[ l dx, = 4qa,a’ ( + ) < 
Se a ss 
AY Ss a , 
tan“! --—-7-— _ (11) 
a a s 


It is interesting to compare this loading with that on the 
corresponding portion of a slender delta wing (Fig. 11). 


Lw = 2 


The pressure is 
(Cy)w = —2ay tan w V1 — (x,/s’)? (12) 
the spanwise lift is 
ly = AgayV s? — x" (13) 


and the total lift carried on the region ABC is 


s\? a ay 
Lw = 2qaa? ( ) cos~!—- — ¥( ) ~~ | (14) 
a Ss a 
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In Fig. 14 the spanwise lift as given by Eq. (13) is com- 
pared with that given by Eq. (10). For a small 
value of a/s there is little difference between the two 
curves. Linear theory values corresponding to Eq. 
(13) are also shown. It is possible that, for configura- 
tions in which the Mach Number effect cannot be 
neglected, there is similarly not too much difference 
between the lift distribution on the outer region 
of a delta wing and the lift distribution on the 
wing of a combination with a circular-cylinder fuse- 
lage. 

Thus, for zero afterbody, the linearized theory will 
undoubtedly give less lift than will the slender-body 
theory. A similar difference between linear theory 
and slender-body theory may be illustrated by a wing- 
body combination for which the wing area has shrunk 
nearly to zero (s = a) though the wing apex angle w 
remains fixed (Fig. 15). For this tiny wing the body 
acts as a vertical wall. Then, linear theory predicts 
that the lift carried by the wing-body combination is 
just twice the lift that the wing alone would carry. 
(The angle of attack induced by the body upwash 
equals the geometrical angle of attack.) None of the 
wing lift is carried on the body. Slender-body theory, 
on the other hand, predicts that, for a wing of vanishing 
size, the lift of the wing-body combination is four times 
the lift that the wing would carry if it were alone in the 
free stream. This is shown by Ward in refer- 
ence 4. 

According to both theories, the lift carried on the wing 
is doubled by the presence of the body. However, 
according to slender-body theory, an equal amount of 
lift is also carried on the body. If, on the other hand, 
one decreased the wing span by letting w — 0 while 


1951 


keeping the wing area fixed, the theories would give the 
same result, since the combination is approaching more 
nearly a “‘slender”’ configuration. 

Theorem (2) may be proved just as easily for non. 
planar as for planar systems; hence, for the configura. 
tion in Fig. 16 slender-body values will be realized as 
the length becomes infinite. The examples of the tri- 
angular wings suggest that, for a finite length, even 
more lift than is predicted by slender-body theory may 
be realized. Finally, the case for which the afterbody 
is finite and has only a cylindrical part (Fig. 1 without 
tail surfaces) may be compared with the corresponding 
flat-plate wing discussed in Section (3). It is to be ex. 
pected that, for zero afterbody, the lift is lower than 
that predicted by slender-body theory. Then, with 
increasing afterbody length the lift will increase and 
possibly exceed the slender-body value. Beyond a 
certain optimum length, added afterbody will be 
detrimental to the lift. Although quantitative esti- 
mates are not obtained this way, wind-tunnel meas. 
urements may be studied in the light of these qualita- 
tive arguments. Of course, these arguments neglect 
nonlinear effects and viscosity. 
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Wind-Tunnel Tests for Temperature 
Recovery Factors at Supersonic Velocities 


W. F. HILTON* 
The Johns Hopkins Unwersity 


(1) SUMMARY 


Temperatures were measured on the surface of a tangential 
ogive attached to a 4-in. diameter cylinder. The tests were 
made at Mach Numbers of 1.73 and 2.00 in the Ordnance Aero- 
physics Laboratory wind tunnel in Texas. The model was 
made of Lucite, a thermal insulator, and the tests showed that 
the surface temperature experienced 88.0 per cent of the adiabatic 
stagnation temperature rise of the airflow just clear of the bound- 
ary layer. The experimental scatter amounted to about !/2 per 
cent. These results compare with Squire’s theory, which gives 
89.42 per cent for a turbulent, and 84.56 per cent for a laminar, 
boundary layer. 

The observed temperatures gave no indication of the position of 
transition, which is thus at variance with the work of Hilton and 
Wingham at high subsonic speeds. 


(II) INTRODUCTION 


|" ATTEMPTING TO FLY at progressively higher speeds, 
we have so far been limited mainly by the aerody- 
namic drag exceeding the thrust of the propulsion unit. 
However, it would appear that we shall be able to 
increase the speed of aircraft and missiles by attention 
to detail, and the question arises as to what will ulti- 
mately determine the greatest speed at which we can 
move through the air. Undoubtedly, the answer to 
this question lies in the aerodynamic heating due to skin 
friction and adiabatic compression that occurs on the 
surface of all vehicles moving at high speeds through 
air. 

This problem breaks down into two main divisions: 
(i) The temperature attained by an unheated body, 
and (2) the temperature attained by a heated body and 
the rate of heat transfer to the airstream. Only the 
first of these is treated in the present note. 


(III) PuystcaL NATURE OF AERODYNAMIC HEATING 


The air in contact with a moving body will experience 
a temperature rise even if brought to rest isentropically. 
In addition, the boundary layer will generate extra fric- 
tional heat by the rubbing of one layer upon another. 
This is offset by the direct conduction of heat from the 
boundary layer to the free stream across the narrow 
boundary layer. In the steady state there will be zero 
heat flow into the surface, which will take up an energy 
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level higher than the free stream and thus be hotter 
than the free stream. The nondimensional parameter 
that expresses the balance between a frictional heating 
and the conductive cooling is called the Prandtl Num- 
ber (c) and is given by 


(1) 


where uw and K are the coefficients of viscosity and 
thermal conductivity and c, is the specific heat. 

If the Prandtl Number is unity, the frictional effect 
just balances the heat-conduction effect, so that the 
heating of the surface would correspond exactly to the 
isentropic stagnation heating obtained when air is re- 
tarded from a velocity V to zero velocity on the surface. 
This isentropic temperature rise is given by the equation 

Ty) — T = (1/2) (V?/ep) (2) 

In the case of air, cp has the value 6,002 ft.” sec.~* per 
°F. and the Prandtl Number = 0.715, indicating that 
the thermal conductivity of air is of greater relative im- 
portance than the frictional effects. Thus, we should 
expect to find that a body placed in an air stream would 
not experience the full isentropic rise of temperature be- 
cause of this extra thermal conductivity. 

For large temperature differences, the effect of heat 
radiation from the hot surface should also be con- 
sidered. No attempt has been made to allow for this in 
the analysis given below. 

The rise of temperature 67 through a laminar bound- 
ary layer has been shown by theory’ to be given by 


(1/2) (y — 1) M*To’” (3) 


o = uc,/K 


5Tiem. = (1/2) V2o "*/cy = 


or 


25T iam lp/V? = o @ = 0.8456 for air (4) 


where 67 is the difference of temperature between the 
surface and the airflow just clear of the boundary 
layer and M is the Mach Number just clear of the sur- 
face at a point where the absolute temperature is T. 
For a turbulent boundary layer, the temperature rise 


is shown by 


iT = = ee i ; (y — 1)M°To'”* 
(5) 
or 
25T rurh. Cp/V? = o “ = 0.8942 (6) 
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Thus, we expect to recover at the surface 89.4 per cent 
(or 84.6 per cent if laminar) of the temperature rise 
that would be experienced by the air just clear of the 
boundary layer if brought to rest adiabatically. 


(IV) EXPERIMENTAL WORK IN THE O.A.L. WIND 
TUNNEL 


An accurate determination of the temperature taken 
up by a solid surface under conditions of supersonic 
airflow will be essential if we are to fly vehicles at or 
near the point where the skin of the vehicle will be on 
the point of thermal failure. Accordingly, an unheated 
model was made and tested in the O.A.L. Daingerfield 
wind tunnel. The model was constructed from a single 
piece of Lucite (part of a sheet 4 in. thick), and its exact 
shape is given in Fig. 1. Lucite (Perspex) was used in 
the construction of the model to reduce the effects of the 
thermal conductivity through the model itself. The 
temperature was measured at 23 stations along the sur- 
face of the model by means of small silver studs to which 
thermocouples were attached at the rear. Tests were 
conducted at a Mach Number of 2 and angles of yaw 
of —2°, —1°, 0°, 1°, 2°, and 4°. Angles of yaw of 4°, 
2°,0°, —2°, and —4° were tested at M = 1.73. The 
thermocouples were connected in series so that the dif- 
ference between one thermocouple and the next was 
measured directly on the instrument. This reduced 
fluctuations due to slight variations of conditions from 
time totimeinthetunnel. — 

It was found that painting a bright silver stud with 
dull black paint did not produce a measurable change 
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Comparison of experimental temperatures with theory; 
and M = 2.0. 


a=0 
in reading. Hence, the effects of radiation were neg- 
ligibly small. 


(V) REDUCTION OF EXPERIMENTAL RESULTS 


The flow over the unyawed model is axially symmetric 
and was calculated by the method of characteristics as 
applied to axially symmetric flow.* Accordingly, the 
readings on the unyawed model were reduced in much 
greater detail than was possible when the model was 
yawed. 

Two methods of presenting the observations have 
been adopted. In the lower part of Fig. 1, the raw 
data are plotted to show measured temperature dil- 
ference against axial position of thermocouple button. 
The theoretical curves for both laminar and turbulent 
flow, calculated from Squire's theory as given in Eqs. 
(3) and (5), are also shown in Fig. 1. 

Put in the more useful forms 


= si - M?T) (1 — o ”) 
AT seen. = T» — 7; = 9 
M? + 2/(y — 1) 
and 
M?T) A — o “) 


AT wv, = Io — T. = = 
M? + 2/(y — 1) 

In Fig. 2, the experiment and theory have been re- 
duced to the form of a nondimensional temperature 
recovery factor Cy, which is defined below in Eq. (7). 
It was defined in this way because it reduces theories 
of the type of Squire’s to the form C, = constant and 
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is thus convenient to use and to remember. We de- 
fine 

C, = (7; — T) (To — T) (7) 
where 


T) = stilling chamber temperature or ‘‘total’’ adi- 
abatic compression temperature (invariant 
even through a shock wave) 


T = ‘‘static’’ temperature just clear of the boundary 
layer, as shown by thermometer moving with 
the air 

T, = surface temperature as measured 

This definition of temperature recovery factor is 


different from that used by Eber,‘ who defined a co- 
efficient C given by 


ty = CF Dae T;) (To = T1) 


where 7 is the ‘‘static’’ temperature of the free stream 

ahead of the body. 
Continuing with the definition given in Eq. (7), we 
know from fundamental thermodynamic relations that 
a? = yRT = (y¥ — 1) GT 


and we may thus rewrite Eq. (2) to read 


T) — T = (1/2) V2/e> = (1/2) Vy — 1) T/a? = 
(1/2) (y — 1) M?T (8) 


Since the free-stream temperature 7 cannot be meas- 
ured by any known device, we eliminate it between 
Eqs. (7) and (8) and rearrange to find 
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Fic. 2. Experiment and theory reduced to nondimensional 
temperature recovery factor. 
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Fic. 3. Direct shadow spark photograph of flow pattern at 
a = 0°, M = 2. Note apparent transition about halfway along 
ogive. 

C , ( 2 A= T; (1 + (1 2) (¥ — 1) M?] 

nl y-1 T) M? 
(To — T;) [M? + 2/(y¥ — 1) 
— 0 ) | Y (9) 


TM? 


Eq. (9) expresses the temperature coefficient in terms 
of M (which may be calculated), the observed tempera- 
ture 7», and the observed temperature difference (7) — 
T;). This equation was used in reducing the observa- 
tions to the The 
Mach Number distribution around the model was ob- 
tained by the method of axially symmetric character- 
mentioned above. The reduced readings are 
shown in Fig. 2 in the form of temperature recovery 
factor against axial distance along the model. 

In order to reproduce Squire’s theory on Fig. 2, we 
take Eq. (3) for a laminar boundary layer and write it 


nondimensional coefficient form. 


in the form 
67 = 7,-T= (1 
But from Eq. (7), defining the temperature recovery 
factor, 
T, —T (1/2) V20'*/c, : 


A. — oe ~~ — 
lam. To - T 


~ > (10) 
(1 2) “/ Cp 


Similarly, the theoretical expression for a turbulent 


boundary layer may be written 


=o" (11) 


%turb. 


(VI) Discussion oF RESULTS 

It is at once apparent from a study of Fig. 1 that the 
readings are consistent to about 0.2°F. and, hence, from 
Fig. 2, that the uncertainty in temperature recovery 
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factor C, is only of the order of '!/. percent. This is re- 
garded as extremely satisfactory. 

The value of the temperature recovery factor is 88.0 
per cent for these experiments. This is a little cooler 
than Squire’s theoretical value of 89.42 per cent for a 
turbulent boundary layer but much warmer than the 
theoretical figure of 84.56 per cent for a laminar bound- 
ary layer. 

Fig. 3 shows a photograph by the direct shadow 
method of the flow at J = 2 over the cone at a = 0° 
and taken at the same time as the temperature meas- 
urements on which Fig. 2 is based. It was discovered 
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independently by H. H. Pearcey at N.P.L., England, 
by workers in Germany during World War II, and by 
Moffett Field, N.A.C.A., that it is usually possible tp 
determine where transition occurs from laminar to 
turbulent boundary layer by inspection of flow photo. 
graphs by the direct shadow method. There is a white 
line parallel to the surface of the body which is charae. 
teristic of a laminar boundary layer, and this white line 
disappears or becomes confused in nature in the case of 
a turbulent boundary layer. A study of the photo. 
graph shown in Fig. 3 indicates that such a transition 
appears to occur on the ogive about 4 in. back from the 
pointed nose. 


This is not consistent with the temperature ob. 
servations, which do not exhibit a fluctuation at this 
position, but is mentioned for the sake of completeness, 


(VII) ConcLusION 


It would appear that on an ogive of initial half angle 
24.91° at a Mach Number of 1.73 or 2, we can expect 
a temperature recovery factor (temperature coefficient) 
of 88 per cent of the adiabatic temperature rise calcu- 
lated on the basis of the velocity just clear of the bound- 
ary layer at the point in question. This nondimen- 
sional recovery factor defined in Eq. (7) appears to be 
a useful concept. 
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Structural Load Measurements on Complex 
Aircraft Components Using Strain-Gage 
Summation Circuits 


MAURICE A. WALTER* 
Naval Air Experimental Station, Philadelphia 


SUMMARY 


A method for the measurement of incremental structural 
loads on complex aircraft surfaces without previous knowledge of 
load distribution is presented. The judicious location of electric 
strain gages and a suitable electric circuit with a single channel 
output to a recording mechanism are the means of effecting the 
measurement. Since the output of the strain-gage circuit is 
generally fed into an amplifier and then to a recorder, equations 
for a bridge looking into a high impedance are developed for 
measuring loads in a complex surface, and application of the 
Both design and experimental data are 
shown. An ap- 
proach similar to the one presented can be used for other types of 


equations is illustrated. 
The electrical equipment used is described. 


strain-gage circuits. 

Static test calibrations using the method have indicated a 
maximum error of 3.7 per cent.in a complex surface, such as a 
horizontal stabilizer of a fighter airplane. 


NOMENCLATURE 


g = gain factor 
K = electric strain-gage factor 
K = (with subscript) a structural-electrical constant to re- 
late mechanical load input with voltage output 
L = lower strain gage 
M = bending moment 
P = shear load 
T = torque 
U = upper strain gage 
V = output voltage of summation circuit 
Vas = bridge input voltage 
Veo = bridge output voltage 
X = torque arm 
x = interaction factor 
Y = moment arm 
Subscripts 
L = left side of horizontal stabilizer 
R = right side of horizontal stabilizer 
Note: Other symbols and subscripts are defined in the text. 


INTRODUCTION 


4 THE TEXT THAT FOLLOWS, it is shown how electric 
strain gages are used to enable simple measurement 
of net loads on an airfoil surface of complex structure 
Without previous knowledge of the load distribution. 
The development of the technique was carried out at 
the Cornell Aeronautical Laboratory, Inc., as part of a 
dynamic stability flight-test program on B-25J and 


Received December 16, 1949. Revised and received June 9, 
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F-80A airplanes under contract to the Aircraft Labora- 
tory, Air Materiel Command, U.S.A.F.* 4 

Measurements of incremental vertical tail, horizon- 
tal tail, rudder, and elevator air loads were desired in 
connection with the above programs. Space limitations, 
particularly in the F-80A, ruled out a pressure-survey 
installation. In addition, some method that would 
eliminate the laborious data-reduction procedure of a 
pressure survey was desired. 

The problems were solved by means of a summation 
electric strain-gage circuit, which enabled the measure- 
ment of the net incremental load on an airfoil surface 
through a single electrical output signal. Air load was 
then obtained by correcting the net load for inertia 
effects. 

The technique stems from a method developed by 
I. G. Rosst for determining shear from a single elec- 
trical output of a strain-gage rosette in a bridge cir- 
cuit.' Further developments enabled determination 
of shear and bending moments with a single measure- 
ment on a structure consisting of more than one load- 
carrying component. Properly located strain gages 
were electrically combined into a single bridge circuit, 
using fixed resistors as weighting factors to make the 
signal output independent of load distribution.® The 
technique was further improved when I. G. Ross and 
W. Hirtreiter{ developed an electronic device to replace 
the fixed resistors.” 

The following text carries the development of the 
technique to measurement of loads in complex surfaces. 


SHEAR-LOAD MEASUREMENT IN A COMPLEX SURFACE 


If, for example, it is desired to measure the incremen- 
tal shear load on the horizontal stabilizer and it is only 
possible to mount gages in the surface itself, the most 
logical points to mount the strain gages would be the 
webs of the beams. However, the strain in the web is a 
function not only of shear load but of torsion, relief due 
to bending of a tapered beam, and secondary effects 
such as interaction between right- and left-hand sta- 


bilizer. 


+ Chief of Wind Tunnel Department, Cornell Aeronautical 
Laboratory, Inc. 

t Chief of Instrumentation Section, Flight Research Depart- 
ment, Cornell Aeronautical Laboratory, Inc. 
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No general equations can be written which will sat- 
isfy all cases for a complex structure. Each installa- 
tion must be considered as an individual problem. The 
factors to be considered in writing the equations should 
be determined from static tests of the proposed installa- 
tion. However, the method of attack will be similar 
for all cases. Therefore, an example will suffice to illus- 
trate the method. 

A measurement of the incremental shear load on the 
horizontal tail of a F-SOA was desired. Static tests of 
strain-gage installations in the aft fuselage section and 
horizontal stabilizer indicated the most suitable loca- 
tions for measuring the desired shear load were the two 
beam webs of the horizontal stabilizer. Since it is 
known from stress-analysis theory that the strain in 
the web is a function of torsion and bending in addition 
to shear, strain gages were installed on the flanges to 
measure bending in the same area as the web gages. 
Unfortunately, it was not possible to install gages 
that would respond to torsion only. 

The characteristics of each gage pair are then de- 
termined by selecting several points on each side of the 
stabilizer and applying a series of loads at each point, 
first each side separately and then both sides symmetri- 
cally. A minimum of two gages are generally mounted 
on the structure forming two arms of a bridge to com- 
pensate for temperature effects. The response of each 
gage pair is measured separately, but all pairs are 
read in sequence under a steady load with a multiple- 
selector strain reader. From these data, plots of 
strain vs. load are drawn and the slopes determined. 
Other tests are made to determine the degree of super- 
position when two loads are applied at points on the 
same side of the stabilizer, the effect of fin loads on 
horizontal stabilizer strain gages, and the scatter of 
data. For the F-SOA, tests showed that superposition 
occurred within +2 per cent; fin loads had no appreci- 
able effects on the stabilizer gages except the front 
beam flanges, and data scatter was small. 

The results of the main test showed that: 

(1) All strain-gage signal responses were directly 
proportional to load with the exception of loads applied 
at the trailing edge of an inboard station. 

(2) A shear and moment applied on one side of the 
stabilizer produced linear signal responses in web and 
flange gages, respectively, on both sides of the stabilizer. 

(3) Though gages were symmetrically located on 
both sides as nearly as possible, sensitivities of corre- 
sponding gage pairs were not identical. 

(4) Flange and web gages responded appreciably 
to torsion. 


Equation of Bridge Output Voltage 


The active strain gages generally form two arms of 
a Wheatstone bridge. In flight-test work, the flexi- 
bility of an amplifier system with an injected carrier 


voltage is of great advantage. The other two arms are 
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usually the source of power for the bridge, and bridge 
output is fed to an amplifier. Consequently, the 
bridge looks into a high impedance. For small changes 
in bridge resistance, the case in strain-gage work, it can 
be shown that 


Veo = |. (load) (1) 


where the load can be shear, moment, or torque. 


Method of Measuring Shear Load 


The method for measuring the incremental load js to 
cancel out all components of the strain signals in the 
web gages not due to shear load. Flange gages are 
used for canceling out effects of bending. Since no gages 
could be mounted to measure torsion directly, an as- 
sumption must be made—namely, it is assumed the 
shear loads carried in the front and rear spars are ina 
constant proportion except for the shear due to torsion, 
This is based on the theory that no twist takes place 
when load is applied on the surface; the beams carry 
all the bending loads, and the shear boxes carry all the 
torsion. In an actual structure, twist does take place, 
but the assumption is close enough for our purposes. 
It is now possible to set up the equations of the load- 
measuring circuit, called a summation circuit. The 
simplifying assumptions used in the equations yield 
approximate results, which are further refined by ex- 
techniques. Therefore, these equations 
will serve to indicate the 
mental techniques of setting up a suitable summation 


perimental 


direction for  experi- 


circuit. 


Setting Up of Equations and Solutions for Gain Factors 


The signal outputs of the strain gages mounted on 
the beams without amplification are expressed as fol- 
lows: 

Right front beam web— 


Veo = KirgMpe + Kor(Pr + Pixir) + KarTr (2) 
Right rear beam web— 
Vener = KigMpt+ Kor(Pr + Prixor) + KeorTr (3) 
Right rear beam flange— 

Vevsr = Kir( Mr + XsrMz) + KerTr (4) 


No equation is written for the front beam flanges, 
since static tests of the F-SOA horizontal tail indicated 
that these gages would not be reliable; this may not 
be true for other cases. 

Equations for the left side are identical in form. The 
constant coefficients (K’s) are functions of the electri- 
cal circuit and structural properties of the member. 
These are determined from static tests by a method 
shown later in the text. 

The term Kgg7 of Eq. (4) presents a problem, since 
no provisions are available to correct for this term. A 
preliminary study of the static-test data indicated that 
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q maximum error of 3 per cent would result from 
neglecting this term and Ks,7,. Consequently, errors 
in experimental results may be partially explained by 
this neglect. 

If the six pairs of gages are connected in the circuit 
of Fig. 1 and are modified to equalize signals on both 
sides of the stabilizer, the resulting voltage output will 
be 


V = (Veoirgir + Veverger + Vevar¥sr)Zr + 
(Venugiz + Veverger + Vensr¥sr)g2 (5) 


The gain factors, gir, Zor, Zar, ANd Biz, Lor, ¥sz are to be 
related so that the voltage output of each side is pro- 
portional to the shear load on each side. The gain fac- 
tors, Ze and g,, equalize the signal output of each side to 
a given load, thereby enabling direct addition of the sig- 
nals electrically. 

The right-hand terms of Eqs. (2), (3), and (4) for 
left and right sides of the airplane, neglecting Ks7 
terms, are substituted. In the right side of Eq. (5), all 
terms are collected as factors of Pr, Pz, Mr, Mz, Tr, 
and 7;. The coefficients of the Mp, Mz, Tr, and T, 
must equal zero in order that circuit signal output be 
proportional to shear load (P = Pz + P,). It can 
then be shown that, in order to satisfy the above con- 
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ditions, 
or = —(Kor/Keor)£ir (6) Fic. 1. Basic electrical circuit for summation. 
Similarly, 
a ( K3, Ker) gun (7) 
£rXar( Kirin + Karger) — 2r(Kirgir + Karger) 
83, = ; (8) 
Kirgi(1 — XsrXsz) 
Kirgirgre + Karger8r + KirXsr 83181 ; 
= 2 (9) 
Kirgr 
— Korfir + Korgor — KorXingir — KsrX2rZ2r g (10) 
c=} : : : R 
Ko 8in + K 51821 = KopXu8in = K5rXer8or 
For simplification, let Determination of K,, Kz, K;, Ks, Ks, Ke L & R 
: Eqs. (2) and (3) are utilized in determining A, to 
Ziv = Zn = £r = | (11) eo , 7 Z ae 
K, for the left and right sides. The method for solving 


All the required gain factors have now been obtained 
and are expressed in terms of structural characteristics 
and electrical circuit constants, the K and x values, of 
the horizontal stabilizer. If it is desired to measure 
the shear loads on each side separately, Eq. (5) is 
written as two separate equations, one for the right side 
and one for the left. In this case, g, and gp are 
dropped. Eqs. (2) through (4) are then substituted in 
the equation for each side, and the g values are deter- 
mined as for the case illustrated. 


In the following paragraphs, methods for determin- 
ing the K and x factors are presented. 


Kir, Kop, and Az only is developed here; the other 
values are determined in identical fashion. 

Let us apply loads at m locations, separately, on the 
right horizontal stabilizer. From the results of these 
tests the slopes (Vepir/Pe)1, (Venr/Pr)e, - . 
Pr), are determined. The subscripts outside the 
parentheses represent the different locations. First, 
Eq. (2) is modified so that the output per unit shear 
load is expressed. P, = 0, and, dividing by Pp, the 
following is obtained: 


i} com) - ( Mr) n > > 
= K Ke K; 
( Pr n 7" (Pr)n + ™ ° ” 


.» (Veoir+ 


(T'x)n 


(12 
(Pr)n ) 
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But (Mp)n = (Pr)nYn, where Y, is the moment arm 
at station m and (Tr), = (Pr)nXn, where X, is the 
torqué arm at station m. Therefore, Eq. (12) becomes 


(Venir/Pr)n = Kir(Vr)n + Kor + Koar(Xpr)n (18) 


With m equations of the form of Eq. (13), values of 
Kir, Ker, and K3z can be determined using the method 
of least squares. These A values are used for comput- 
ing the gain factors. 

In order to know the X distances, the elastic axis 
must be determined, either experimentally or analyti- 
cally. As used here, the elastic axis is defined as that 
line on which a load may be applied without causing 
twist of the stabilizer at the station of load application. 
The X distance is measured parallel to the airplane’s 
centerline from the point of load application to a line 
passing through the elastic axis at the gage installation 
station and perpendicular to the airplane’s centerline. 

Similar equations can be written to solve K4z, Ag5p, 
Ker,-Kiz, Kor, Kar, Kaz, Koz, and Kez. 

Kiz and Kz, are determined in like manner from Eq. 
(4). 


Determination of x Factors 


Because of redundancy of fuselage and tail structure, 
a load applied on the right side of the horizontal tail 
produces a response on both sides of the tail, including 
that part of the left side of the horizontal tail outboard 
of the fuselage attachment fittings. 

The interaction factor, x, is defined as the fraction of 
the signal sensitivity of a gage pair (Ke, Ks, Az) due to 


a load placed on the opposite side of the tail. By this 
definition, 

Xin = Korr/ Kor (14) 
and 

Kort = Veoir/Pr (15) 


where the RL subscript signifies a response on the right 
side of the stabilizer due to a load on the left side. 

Similar equations are written for x22, X3z, and x fac- 
tors for the left side. 
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Fic. 3. 
and elastic axis. 


gage installation location, 

Data for calculating the x factors are obtained by 
applying loads to the stabilizer on one side of the air- 
plane and recording the responses of the gages on the 
other side. 

The development of the circuit equations and the 
methods for determining the constants of these equa- 
tions enable the actual design of a summation circuit 
for measuring horizontal tail load. However, before 
taking up the actual design of such a circuit, the equip- 
ment used will be described. 


DESCRIPTION OF EQUIPMENT 


The basic circuit of the instrument enabling measure- 
ment of load on a complicated structure with a single 
channel output is shown in Fig. 1. The circuit diagram 
for a simple electronic summation unit constructed by 
the Cornell Aeronautical Laboratory, Inc., is shown in 
Fig. 2.2 As shown, the strain gages are connected at 
points A, B, and C, where they compose two arms of 
a Wheatstone bridge. The other two arms, consisting 
of an adjustable resistor, are used for initial balancing 
of the bridge. The output signal of each bridge is fed 
to an attenuator that is connected to ground. The at- 
tenuator will pass any desired fraction of the output 
signal to the grid of the isolating amplifier. Each am- 
plifier tube plate feeds the common load resistance, R; 
therefore, the voltage drop across R will be propor- 
tional to the algebraic sum of the signals across the 
grids of the six tubes. The voltage across R is then fed 
to the amplifier immediately preceding the recording 
device. 

The settings of the attenuators are governed by the 
gain factors calculated from the preceding equations. 
In setting attenuators, the highest gain factor is arbi- 
trarily given the The other gain factors 
are then given proportional values. These decimal 
values are then set into the attenuators of the proper 
channels. 


value of 1. 


In practice, it is common to supply the bridges with 
an injected carrier voltage of 1,000 cycles per sec. 
With such a system, compensating condensers are re- 
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STRUCTURAL 


quired to enable bridge balance, even if all wiring is 
shielded. To minimize electrical “‘hash,’’ tubes should 
be shock-mounted. 


DESIGN AND EXPERIMENTAL DATA 


Static tests were carried out on the horizontal sta- 
bilizer in order to calculate the summation circuit gain 
factors defined by Eqs. (6) to (10). During the tests, 
the output of each strain-gage circuit was recorded in- 
dividually for a given load. Plots of electrical output 
ys. load input were made, and the slopes (Vcp/P) were 
The load application test points on the 
The elas- 


determined. 
surface of the stabilizer are shown in Fig. 3. 
tic axis was determined experimentally and is also shown 
in Fig. 3, together with the station location of the elec- 
tric strain gages. 

An analysis indicated that it would be sufficiently 
accurate to measure torque and moment arms parallel 
and perpendicular to the plane of symmetry of the air- 


plane. 


Computation of Summation Circuit Output 


The gain factors were computed according to the 
equations already developed and static tests. The re- 


sulting gain factors are listed in Table 1. 


TABLE 1 


Comparison of Experimentally Determined and 
Calculated Gain Factors 


Gain Determined 

Factor from Test Calculated 
gir’ 0.73 0.74 
gor’ 0.34 0.40 
gsr’ 0.17 0.10 
gin’ 1.00 1.00 
gor’ 0.34 0.42 
21" 0.22 0.16 


The data of Table 1 can now be used to predict the 
output of the summation circuit for loads applied at 
each point. Only the symmetrical case is shown here. 
Eq. (5) can be rewritten in a more convenient form by 
dividing through by P, multiplying all right-hand gain 
factors by ge, and similarly for g,. The equation then 


becomes 
V Vepir ’ Vcper , Vcpsr , 
P = P £Zir + P Lor + P £3r + 
= —_ Vom 
* CDIL CD2L CD3L 
Zu," = +- a (16 
P gn + P £21 P 83x ) 


Substituting the values of the calculated gain factors 
in Table 1 and the values of Vepnr/P and Vepn,/P de- 
termined from static tests into Eq. (16), the resulting 
signal outputs of the circuit are plotted in Fig. 4. The 
points are plotted on a plan-view diagram of the hori- 
zontail tail at the points corresponding to the load ap- 
plication points. Lines of constant sensitivity are then 
drawn—therefore, the name ‘‘contour plots.” 
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Fic. 4. Horizontal tail contour plot of sensitivities. 
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STA. 25 


Contour plot of sensitivities from test data for hori- 
zontal tail using calculated gain factors. 


Fic. 5. 


An inspection of the figure reveals that a 10 per cent 
variation in sensitivity is evident with sensitivity de- 
creasing as the load moves outboard. This characteris- 
tic indicates insufficient correction for the effects of 
bending moment. The gains of the strain gages were 
adjusted experimentally to give improved results as 
explained below. 


Experimental Check and Adjustment of Summation 

Circuits 

First, a preliminary analysis of the effect of varying 
each gain setting on the entire summation circuit was 
made in order to minimize the number of trials and indi- 
cate the most logical way of adjusting the gain settings. 
It can be shown that a change in sensitivity, (V/P), 
of the summation circuit for changes in the gain factors 
can be expressed as follows: 


CD3R 


V\ Veo ear 

d (5) = - dgip’ + ; oa d gor + P dg3p + 
Tl ys Visti Wel ame 
dgi,’ + f : dg," + ~~ dgsz, (17) 


Eq. (17) shows the effect of varying the gain settings 
can easily be investigated using the unamplified chan- 
nel sensitivities (V¢p/P) determined in the static tests 
for each loading point. Subsequent tests justified the 
results of the investigation and materially cut down the 
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guesswork in adjusting the gain setting of the circuit. 
The final gain factors giving the best results are listed 
in Table 1, together with the calculated gain factors. 
This table shows that calculations yield a reasonable 
approximation of the final values. 


Comparison of Calculations with Experimental Results 


A contour plot of the measured sensitivities using 
the experimentally adjusted gain factors is shown in 
Fig. 6, while a plot of the measured sensitivities using 
the calculated gain factors is shown in Fig. 5. A dis- 
tinct improvement in uniformity of sensitivity is evi- 
dent in Fig. 6 outboard of station 40. Although the 
sensitivities at station 25 depart considerably from 
those outboard of station 40, the average over the 
chord at station 25 is in line with results outboard of 
station 40. 


Effects of Load Distribution on Accuracy 


Fig. 6 shows that the summation circuit is somewhat 
sensitive to load position. The effect of this characteris- 
tic on accuracy of measurement was investigated for 
two extreme air-load distributions and was computed 
at +2.4 per cent. This error is probably due to the 
neglect of torsion in the flange gages, as was brought 
out earlier in this paper. The only other major source 
of error detected was imperfect superposition, measured 
as 1.9 percent. Probable error was estimated to be 
+3.1 per cent. Actual static tests indicated a maxi- 
mum error of +3.7 per cent, reasonably substantiating 
the estimated error. 

Fortunately, in dynamic test work, incremental 
loads of short time duration are generally of interest; 
consequently, drift of electronic equipment does not af- 
fect data accuracy. However, where loads above a 
fixed datum are desired, drift can introduce large errors 
into the data. Therefore, equipment must be designed 
to maintain drift at a minimum. 


CONCLUSIONS 


A method has been developed for measuring struc- 
tural loads on a complex airfoil surface without pre- 
vious knowledge of load distribution through one chan- 
nel of intelligence, which is particularly useful in the 
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zontal tail using experimentally adjusted gain factors. 


measurement of net loads in flight under dynamic con- 
ditions. In the specific case illustrated, a 
mum error of +3.7 per cent, determined in static tests, 


maxi- 


was incurred in the measurement of loads using this 
method. 

The method eliminates laborious data reduction, 
and the equipment necessary to effect the summation 
of all electrical outputs into a single channel of informa- 
tion, requires little space. 
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Conical Flow Methods 


Applied to Uniformly 


Loaded Wings in Subsonic Flow’ 


SEYMOUR LAMPERT?+ 
Ames Aeronautical Laboratory, N.A.C.A. 


SUMMARY 


Conical flow solutions for Laplace’s equation are employed 
which yield expressions for the contour of uniformly loaded lift- 
ing surfaces at subsonic speeds. By properly superposing these 
conical flow fields, it is possible to determine the twist and cam- 
ber for uniformly loaded wings having various aspect ratios, 
taper ratio, and sweepback. The solutions derived may also be 
used to investigate wing contours necessary for various types of 


spanwise and chordwise load distributions. 


INTRODUCTION 


I A PREVIOUS PAPER,' Jones has demonstrated the 
possibility of extending the methods of conical flow 
analysis to wing problems at subsonic speeds. Jones 
derived expressions that give the subsonic pressure dis- 
tribution over airfoils having various thickness dis- 
tribution at zero lift. The method used is applied in 
the present paper to the calculation of the twist and 
camber of lifting surfaces of zero thickness. The 
usefulness of these conical flow solutions lies in the fact 
that they are well suited for treating problems involv- 
ing sweptback wings. 


SYMBOLS 
s, 7, 2 = Cartesian coordinates in the stream direction, 


across the stream, and in the vertical direction, 


respectively 


x’, y’,s’ = transformed coordinates correspond to x+ 
Vv 1 — M?, y, s, respectively 

u,v, W = perturbation velocities in the stream direction, 
across the stream, and in the vertical direction, 
respectively 

V = free-stream velocity 

M = free-stream Mach Number 

A.R. = aspect ratio (b?/S) 

b = wing span 

Co = root chord 

d = taper ratio 

A = angle of sweep of the leading edge 

m = slope of leading edge (cot A) 

mt = slope of the trailing edge 

Ms = slope of side edge 

t = slope of any ray from the origin (y/x) 

En» Vn = coordinates of the apex of a uniformly loaded 
sector 

ty = slope of a ray through the apex of a uniformly 
loaded sector, (y — yn)/(X — Xn) 


Received January 10, 1950. 

* The problem was suggested to the author by R. T. Jones, 
who also supplied helpful criticism. 

+t Aeronautical Research Scientist, 6- by 6-Ft. Supersonic 
Wind Tunnel, High Speed Research Division. 


the = slope of a ray through the apex of the uniformly 
loaded sector B, (y — ¥nn)/(X — Xnp) 

tne = slope of a ray through the apex of the uniformly 
loaded sector C, (vy — ync)/(* — Xne) 

Uo = value of streamwise perturbation velocity cor- 


responding to the design lift coefficient, 
4uy V= Cro 


alue of the streamwise perturbation velocity 


ll 
a 


for a uniformly loaded sector of semi-infinite 
extent 
value of the streamwise perturbation velocity for 
the uniformly loaded triangular sector ex- 


ll 
< 


tending to infinity 

value of the streamwise perturbation velocity of 
the uniformly loaded sector for 0< ¢ < m 

value of the streamwise perturbation velocity of 


ll 
e 


ll 
< 


the uniformly loaded sector for ms <¢ < m; 


u = value of the streamwise perturbation velocity of 
the uniformly loaded sector for m <t < 1 

us’ = image of u; reflected in the xz plane 

Wj = the vertical perturbation velocity associated with 


uj, where 7 = 1, 2,3,4,5 

wr = vertical perturbation velocity for a sweptback 
tapered lifting surface 

Wp, We = Vertical perturbation velocities for sectors B and 
C of the tapered wing 

4 = contour associated with the vertical perturbation 


velocity w; 


or = contour associated with the vertical perturbation 
velocit ywr 
° / 9 9 9 . . 
€ = (y + iz)/(x + Vx? + y? + 2?) conical variable 
= . / 9 9 9 : 
‘ = (y — is)/(x + Vx? + y? + 2?) complex conju- 
gate of e 
€0 = value of e which corresponds to m in the plane 


of the wing 
és = value of e which corresponds to ms in the plane 
of the wing 





Fs = a variable used in the analysis, 2e/(1 — e?) 
On(x) = Legendre function of the second kind, 
ve 1 ' x+1 
Co(x) = 2 og =a 
x x+1 
),(x) = -lo — 1 
Cr 2 . x—- 1 
Ao = m1lt+ V1 + m?)-! 
As = m (1+ V1 + m;)7! 
ANALYSIS 


If it is assumed that the stream lines deviate only 
slightly from rectilinear flow and that the rate of change 
of the area of the stream tube in the flow direction is 
small, problems of steady-state incompressible flow of a 
perfect fluid resolve themselves into finding solutions 
that satisfy Laplace’s equation 


Sun + Sy + Sq * 0 (1) 
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Fic. 1. Variation of « with x and y in the chord plane of the wing 


(z — 0). 


where SS is an analytic function of x, y, z, which describes 
any physical characteristic of a perfect fluid. It has 
been shown by Donkin? that, if S is homogeneous of 
degree zero—namely, constant along rays emanating 
from a fixed point in space—it may be expressed, in 
general, as 


y + 12 
s=s(—3 ) 
et Vx? + y? + 2? 


e( J-* ) 2) 
x + V x? + y? + 2? 


where f and g are each analytic functions of a single 
complex variable. 

For the purposes of this paper, S has been chosen as 
the perturbation velocity in the x direction and may be 
expressed as 


S =u = fle) + gle) (3) 


where ¢ and ¢ are the arguments of f and g in Eq. (2). 
Since f and g are both linearly independent solutions of 





1951 


Eq. (1), it is necessary only to work with f(e) in this 
paper. 


The Variable « 


Before entering the problem of finding the functions 
of e which define a uniformly loaded lifting surface, jt 
will be useful to investigate first the properties of the 
variable ¢ itself. We are concerned primarily with the 
variable near the chord plane of the wing (z > 0). At 
z = 0, the value of eis given by the relation 


e= y/(x + V x? + y?) (4) 
where ¢ is seen to be a pure real number constant along 
rays from the origin. In Fig. 1, the real value of ¢ as 
given in Eq. (4) is plotted vertically above the xy plane 
for all values of y/x. 

An analytic transformation that will be useful 
in many cases for prescribing boundary conditions 
is 

1/é = (1 — e&)/2e (5) 


Since, in the plane of the wing where ¢« — €¢ as 
s-+ G, 


1 l — € < 
—> —- (6) 
g 2e y 
for all z 
(1 — e€)/2e = x/(y + 22) (7) 


It is of interest to note here that the graphical con- 
struction relating « and x, y, z as given in Fig. 2 is the 
subsonic analog of the Tschaplygin transformation 
given by Busemann’? for the supersonic case. In Fig. 2, 
the center of the circle is chosen at the origin of the 
coordinate system, where e.and 1/e are projected onto 
the xy plane in the same manner as shown by Buse- 
mann. 

In the present paper, the variables ¢ or e will be em- 
ployed interchangeably, depending on the problem. 
The variable — may be more convenient to use in most 
cases since the real part of & is just y/x = ¢, making it 
easier to prescribe boundary conditions in the plane of 
the wing. 
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Graphical construction of e« and 1/e. 
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Derivation of Uniformly Loaded Sectors 


From the foregoing analysis, it is possible to define 
a conical flow field in the plane of the wing so that the 
joading prescribed will be uniform. 

The boundary conditions that must be satisfied for a 
uniformly loaded surface extending over all values of y 
and positive x are: 

(1) u = 0 for x < 0, which corresponds to —1 > 
¢>1(z— 0). 

(2) u = wu for x = 0 which corresponds to —1 < 
e<1(z—0). 

A function satisfying these conditions takes the 


form 


uo 


u, (€) = R.P. 2 log fare 
€ 


= oR L O €) 
= == §£ d P. 1 
1 of y 


. 


where Qo(e) is the Legendre function of the second kind. 
The region over which the uniform load extends is 
shown in Fig. 3a for z — 0, while the loading plotted 
against € is given in Fig. 3b. 


If the condition of irrotationality is utilized—namely, 


ow/Ox = Ou/dz (9) 
the vertical velocity component that determines the 
shape of the lifting surface may be expressed in terms 
of €as 


dw = —(i/2)[(@ + 1)/e du (e) (10) 


or 


1 2+ 1 du 
w= f: de + constant (11) 


2 e de 


Substituting Eq. (8) into Eq. (11), the vertical ve- 
locity distribution may be calculated as 


Uo l-—é 
w(e) = —— log + constant (12) 
Tv € 
or 
wi(£) = (uo/m)log — + constant (13) 


The real part of Eq. (13) evaluated in the plane of the 


wing becomes simply 
w(t) = (uo/m) log |t| + constant (14) 


The constant as given in Eq. (14) is allowed to remain 
arbitrary, since the boundary conditions will be satis- 
fied by the superposition of solutions of the type given 
in Eq. (8). 

The vertical velocity as given in Eq. (14) is seen to be 
infinite along both the x and y axis. The infinity along 
the y axis (x = 0) is of the type normally associated 
with the two-dimensional camber shape of uniformly 
loaded wings. The infinity along the positive x axis, 
however, is not characteristic of cylindrical flow. 
From the very nature of the solution, though, there 
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Fic. 3a (top). Uniformly loaded surface of semi-infinite extent 
—x2x <t<o. Fic. 3b (bottom). Plot of semi-infinite uniform 
load as a function of e (—1 < «€ < 1). 


must exist a singularity of this form along the x axis, 
which leads one to believe that the solution for the shape 
of a uniformly loaded lifting surface covering half the 
xy plane is not unique. However, it can be shown that, 
by the proper superposition of the solutions given by 
Eq. (8), the two-dimensional or cylindrical vertical 
velocity field corresponding to a wing of finite chord 
may be realized. If one lifting sector of the type given 
in Eq. (8) is placed along the line x = —1 and another 
of equal and opposite strength is placed at x = +1, 
the result is a uniformly loaded infinite wing of chord 
2. The distance and position of the lifting elements 
along the x axis is arbitrary, and the values of x were 
taken at +1 for convenience in evaluating the vertical 
velocity distribution. The resultant distribution for 
the vertical component becomes 

w= a log ed. = —2 = Qo(x) (15) 

T s=~- i T 

which is the well-known two-dimensional solution for 
a uniformly loaded thin lifting ‘wing. It might be 
pointed out that, by keeping the leading edge fixed at 
x = —1 and allowing the trailing edge to move off a 
large distance downstream, the nature of the solution 
remains cylindrical. The contour resulting from Eq. 
(15) may be obtained by the following relation: 


, Uo 


1z 
v(' ) = [w(x)] ,+0 = —2— Q(x) (16) 
dx gr 


which requires that the flow be parallel to the surface 
at the surface. Integrating Eq. (16) with respect to 
x gives the contour as 
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Fic. 4. Camber and slope of a uniformly loaded two-dimensional 
lifting surface for 2u9/V = 0.04r. 
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Fic. 6a (top). Uniformly loaded surface for —m < t < m. Fic. 
6b (bottom). Plot of uniform load in the e plane for —eo < € < &. 
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Uo © 
a(x) = — p2Qi(») + log x? — 1] (17) 
Tv 


which is shown in Fig. 4a plotted against x. It is of 
interest to compare the cylindrical nature of the con. 
tour as given in Eq. (17) for the finite chord wing with 
the conical-type contour 


Upx 
ae oe = — 


l 
(tos ob constant (18) 


7 


resulting from integrating the vertical velocity distri- 
bution given in Eq. (14) for a uniformly loaded surface 
covering half the xy plane. The contour as given in 
Eq. (18) is shown in Fig. 5 and is seen to be a linearly 
increasing function of x along conical rays emanating 
from the origin. 

It will be noted that Eq. (15) gives an infinite ver- 
tical velocity at the leading and trailing edges as a 
consequence of the fact that a finite load is carried at 
the airfoil leading and trailing edges. The flow experi- 
ences theoretical infinite accelerations at the lead- 
ing and trailing edges which are, of course, impos- 
sible in a real fluid. It is to be expected, there- 
fore, that the theoretical loading will not be realized 
in the immediate vicinity of the leading and trailing 
edges. 


Uniformly Loaded Sweptback Wing 


Infinite Span.—In a manner similar to that em- 
ployed for finding the shape of the two-dimensional 
wing, it is possible to determine the proper shape of a 
sweptback wing of infinite span supporting a uniform 
load. 

The boundary conditions for a uniformly loaded in- 
finite triangular sector are 

(1) u = Ofor —m >t >™m, where m is the slope of 
the leading edge. In the e plane, this boundary condi- 
tion is —€ >€ > € (2— 0). 

(2) uw = uw for —m < t < m, which corresponds to 
—@ See (z— 0). 


The function satisfying these boundary conditions 


is 
Uo - e¢~— @& 
Uo(e, €&) = —— R.P. 12 log 
T € €0 
uo € 
= -2"RP.iQ ( ) (19) 
T €0 


The regions over which this function is uniform are 
shown in Fig. 6a, and the distribution as a function of ¢ 
is plotted in Fig. 6b. 

It can be seen that Eq. (8) is merely a special case 
of Eq. (19), where the sweepback angle A approaches 
zero. The vertical velocity distribution resulting from 
Eq. (19) obtained by integrating Eq. (11) is 
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Tv 2€ € 


Uo} 1 + €° 1 
We(€, €&) = — ‘| 7 log (€? — e?) —— log (| + constant (20) 


The real part of Eq. (20) evaiuated in the plane of the wing becomes 


Ug V1 + m? m?* t? 1+ V1 +m? . ] 
w(t,m) = — log 7 — + sinh—! ~ | + constant 
7 m Q+Vit¢m)? (1+ V1+ 8) m t 
(21) 
where m = cot A is the slope of the leading edge and If the apices of two lifting surfaces of equal loading 
t= y/xisany ray from the apex. The vertical velocity but of opposite sign are placed atx = +1, y = O, the 


distribution has a logarithmic infinity along the lead- result will be a sweptback wing of infinite span with 
ing edge and along the x axis as was anticipated. These chord 2 supporting a uniform load. The solution may 
conditions are similar to those obtained previously for be characterized as 


a straight leading edge. However, in this case, the 

infinity along the x axis does not disappear under super- w= —[we(h, m) — walle, m)| (22) 
0sition of solutions. * a ; — 

The result conforms with previous results obtained Fhe values fer t; and &; to Se eanetiaies o Ey, C28) 
at supersonic speeds for a uniformly loaded sweptback — 

wing with subsonic leading edges.‘ It is possible, = y/(x+1); tb = y/(x — 1) (23) 
however, to eliminate the impractically high twist near 

the root chord if the loading is allowed to reduce slightly For a tapered lifting surface Eq. (22) may be rewritten 
from the uniform character in the region of the root. as 

Baldwin’ has demonstrated how this can be done for Wr = —[we(t, m) — We(te, m,)] (24) 


triangular wings with subsonic leaditig edges in super- 


sonic flow. where m, is the slope of the trailing edge. 


The contour corresponding to Eq. (21) may be given as 
| } g 


/ ‘ / ‘ ‘ 
x] I Vi m? l Vi m? i? 
a(t, m, x) = —— | s == (sint 1--—~—V1+ r) + 2 ; log |A? — ] - 
; m (1 + V1 + #*)? 





Vr m 
t? 2t Ao F =a | ” 
- + — Qo (1+ V1+#)1-+ constant] (25) 
1+Vi+e2 ™ t f 
where Ay = m(1 + V1 + m?*)—! and the contour of the Zz 
sweptback wing of constant chord is given by y 





z= —[z2(t:, m) — Zo(te, m)] (26) _ ill 


and the contour of the tapered wing is given by 






Zr = —[22(t, m) — Zo(te, m,)] (27) Rear 3/4 vew 
Eq. (26) is plotted in Fig. 7 for a 45° sweptback wing of 
constant chord. 

Finite Plan Forms.—The method outlined above 
may be used to solve for the shape of tapered wings of 
finite span. The treatment consists of superimposing 
uniformly loaded surfaces so as to cancel all lift outside 
the confines of the wing plan form so that the correct 
boundary conditions are met. For example, a finite He. * Sam ew 
sweptback wing of given span may be obtained from Fic. 7. Camber and twist of a uniformly loaded 45° sweptback 
the uniformly loaded wing of infinite span of the pre- wing for 2uo/V = 0.04. 
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V vious section by canceling the loading beyond the wing 
| tips as shown in Fig. 8. 

The vertical velocity component for region A in Fig, g 

corresponds to the value given in Eq. (24). It remains 

y now to find the vertical velocity components that result 

from superimposing constant load sectors of equal and 


A 
“ opposite sign, respectively, at the leading and trailing 
edges of the tip station. 


f The solution for these sectors, which satisfies the 

< following boundary conditions (1) « = 0 for 0 >t >m 
or 0 > € > & (2 > 0) and (2) uw = wm for0 <t<mor 
0 <€ < «@(z— 0), takes the form 


A=cot™'m 














Fic. 8. Superposition for obtaining a sweptback wing of A.R. Uo ; €— © 
= (b/2)[1 + (1/2)(A — 1)]-! and A = 1 — (b/4)[(1/m) — Uz(€, €&) = — R.P. 2 log (28) 
(1/m:) J. T € 


The region over which this function is a constant is shown in Fig. 9. The corresponding vertical velocity com- 
ponent is given as 












































ae ae Uy {2 V1 + m? t 1+ V1 + m? t 
w3(t,m) = —- log) Ap — 7 — - log 7 
2a m | l+V1l+? m l+V1+2 
l+V1i+? | 
é r -}]-+ constant (29) 
and the contour resulting from integrating Eq. (29) with respect to x is given by 
Ux | 2 Vi m? t t l t V1 + m?\. 
23(t, m, x) = —— : = (m — t) log} Ao — : + + + sinh— X 
27 I m* | 14 W714 2 Ago m Ao 
s J, ; 
l l V1 t?)? V1i-e+P tf Vi1+ m? 
+ “+ eri — of , + constant (30) 
t 4t Ay 1471447 
In sectors B and C, respectively, the vertical velocity y 
components are given in the form | 
We = —[wW(tiz, mM) — wWs(tep, m,)| (31) 
Fe a 
where m I 
— y — (6/2) — y — (6/2) [ \ “y 
1B at sas (b/2m) + he 2B i mes (d, ‘2m,) — i \ ! } 
(32) I A 
“Seal 
and 
U=Uo 
We = —[Ws(tic, —m) — Ws(tec, — m,)| (33) | 
where x m=c™a 
y b/2 y b/2 
in om 2 + ( ) a -_ a + ( ) > ju le) 
x — (b/2m) + 1 x — (b/2m,) — 1 
(34) 
The contours corresponding to the vertical velocities 77, 
in Eqs. (31) and (33) are aa _ - 
7 ~~ AI x2 
Zp = —[2s(tiz, m) — 23(tep, m,)] (35) V// Vy, e 
and “y O & if 











a Fic. 9a (top). Uniformly loaded sector for 0 < ¢ < m. Fic. 
Zo = —[22(tic, —m) — 25(tec, —m,)| (36) 9b (bottom). Plot of uniform load in the « plane for 0 < « < « 
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CONICAL FLO 
The slope of the surface for the complete wing may 
therefore be written as 


V(dz/dx) = wr — (Wg — We) 


(37) 
where the contour will be 


(38) 


z= 27 — (Zp — 2c) 

It should be pointed out here that in actual design 
practice it will not be possible to maintain a uniform 
loading all the way to a streamwise tip, and it is there- 
fore desirable to allow the loading to be rounded off at 
the tip in such a way that impractical wing twist near 
the tip may be avoided. 

This treatment of uniformly loaded wings may be 
used to obtain the twist and camber for a variety of 
wings of different aspect ratios, taper ratios, and sweep- 
back. The rectangular wing, for example, is the limit- 
ing case of the sweptback wing described above. It can 
be shown that, if tlie sweepback of the leading edge 
be allowed to approach zero (m — o), then Eq. (38) 
describes the contour of the surface of a rectangular 
wing of aspect ratio 5/2. 

Wings with Raked Tips. 
tour of a wing having a raked tip, it is necessary, as 
before, to define a uniformly loaded sector that de- 
scribes a uniform load over the region shown in Fig. 10, 
where u = uw form, <t < more See. Thefunc- 


In order to obtain the con- 


tion describing these conditions is 


e€e=— @& 


u , 
Us (€, €0, €;) = — "RP. i log (39) 
vs 


e- & 


If sectors of the form given in Eq. (39) are super- 
posed at the tip station in the same way as previously 


uy { 2 Vi + m? t 
u(t, ms, m) = - log |Ay — — , 
2a m i+ V14+f 


and the contour associated with Eq. (39) will be 


/ 
ux) 2V1+ m? 
Z(t, Ms, m) = {— 


(m — t) log | Ay — —— 
' 1+ V1 a 2 
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Uniformly loaded sector for m, < t m. Fic. 


Fic. 10a (top). < 
Uniform load in the e plane for e, < € < ©. 


10b (bottom). 


outlined for the streamwise tips, it is possible then to 
find the vertical velocity corresponding to the resultant 
streamwise component as 


_2 V1 + m,? lo | t 
ms . 1+ Vi 

l i. 1 
i = ; sinh! i) + constant (4)) 


‘ F i ° 
t 2V1+ m,? ; 
(t{ — m,) X 


4nV| m? m,* 
Vitm? Vitm? Ll. st. 1 
log = —A, 4+] + im — ~ sinh~!— + 
ll+Vil+? mA mA, A, Ao t 
1 ] /- ‘ V1 + m? V1 + m,? 
— V1l+e+ + + constant; (41) 
ms m m Ms \ 
Triangular Plan Form.—The triangular plan form _ sector, take the form 
may be achieved without much difficulty by canceling sto —s 
the lift behind a uniformly loaded triangular sector us(€, 0) = —— R.P. 2 log (42) 
us a. 


described in Eq. (19). The uniformly loaded sectors 
necessary for canceling the loads, behind the triangular 


where u = uw fore) <e€ < lorm <t< ~. InFig. lla, 
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Fic. lla (top). Uniformly loaded wing form <t< ». Fic. 11b 
(bottom). Plot of uniform load in e plane for e) < « < 1. 

















Eq. (42) is shown in the physical plane, and, in Fig. 
lib, it is given in the « plane. The vertical velocities 
and the contours associated with Eq. (42) may be ob- 
tained from Eqs. (40) and (41), respectively, if m, is al- 
lowed to approach infinity and A, to approach unity. 

The required superposition for the triangular wing 
may be characterized as 


u = U(x, ¥) — W(x — Co, ¥) + Us[x — Co, y — b/2)] + 
us’ [x — co, y + (b/2)] (43) 


where u;’ is the mirror image of u;. Its formula is 
given by 
Uu : e+ 1 
us) = — R.P.i log (44) 
T e+ € 
The vertical velocity component corresponding to Eq. 
(43) is given as 
WwW = W(x, y) — wi(x — Co, vy) + ws[x — co, y — (0/2)]+ 
Ws’ [x =— te > + (b, 2)| (45) 


and the contour will then be 


) + 25[x — co, y — (b/2)] + 
Zs'[x — co, y + (b/2)] (46) 


z= 2(x, y) — a(x — Co, y 


The superposition described in Eqs. (43), (45), and (46) 
is shown in Fig. 12. 
PRANDTL-GLAUERT TRANSFORMATION 


The solutions derived in the previous sections are 
those that satisfy Laplace’s equation [Eq. (1)], a well- 
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known linear differential equation of mathematica] 
physics. Solutions of this equation are of special in- 
terest to aerodynamicists, because it can be shown to 
be a special form of the linearized potential equation for 
compressible flow given as 


(1 — M*)Sez + Sw + Siz = 0 (47) 


Eq. (47) reduces to Laplace’s equation if the Mach 
Number of the flow is zero or if, in the case where the 
Mach Number is greater than zero and less than unity, 
the following transformation is applied to the coordinate 
system: 


eg! = —x/V1 — M?; y'’=y; 2’ =32 
in which case Eq. (47) becomes 
Serer + Sytgt + Serge = 0 (48) 


The solutions for S are now harmonic functions of 
the variables x’, y’, 2’. It is evident, therefore, that 
the solutions obtained in the preceding sections for 
M = 0 can be used to find solutions at Mach Number 
M. In order to find the vertical coordinates z of a 
cambered wing to support a specified uniform load at 
M, the procedure is to stretch the x coordinate in the 
equations for J = 0 by the factor 1 V1 — M*. If 
the. value of the slopes (i.e., w/V) for the wing in in- 
compressible flow (J = 0) are transferred unchanged 
to the corresponding points on the wing for Mach 
Number M, then the streamwise component of the 
velocity u (pressure) is increased by the factor 1 + 
V1 — M?. On the other hand, if it is desired to maintain 
a given u (i.e., pressure) at corresponding points of the 
two wings, then the slopes of the wing at Mach Number 
M are to be decreased by the factor V1 — M? of the 
values of the slopes of the MW = 0 wing. 

To illustrate, Eq. (17) given for 1 = 0 will be trans- 
formed to give the camber of a wing at Mach Number 
M to support a specified uniform load uo. 


(Continued on page 138) 
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Fic. 12. Superposition for obtaining a triangular wing of A.R. = 
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Calculation of Thermal Stresses in a 


Wedge-Shaped Wing’ 


GEORGE H. EISENHARDT?# anp WARREN M. ROHSENOW?} 
Massachusetts Institute of Technology 


SUMMARY 


A method for calculating thermal stresses is presented and 
applied to a solid symmetrical double-wedge steel wing flying 
at an altitude of 50,000 ft. The wing, after flying steadily at a 
Mach Number of 1.4 and having a uniform temperature dis- 
tribution of 531.6°F. absolute, is accelerated at a constant rate 
of one Mach Number per 30 sec. until a Mach Number of 6.0 is 
reached. Thermal stresses at the instant the wing reaches a 
Mach Number of 6.0 are calculated using the elastic properties 
of steel at room temperature. An interpolation method is pre- 


sented for estimating the thermal stress during the acceleration. 


NOMENCLATURE 


A = area, sq.in. 

C, Ci, Ce =,.arbitrary constants 

E = Young’s modulus Ibs. per sq.in. 

E’ = E/(1 — v?) 

Fi, = As(T — Tm)n 

ie = moment of inertia of wing cross-sectional area 


about the x-axis 
i, = moment of inertia of wing cross-sectional area 


about the y-axis 


K = —(aE + v'a’E’) 

M = Mach Number, wing velocity + sound velocity 
in main stream 

7 = actual temperature at a point 

Ti = initial temperature of the wing = 531.6°F. abs. 

Tn = mean temperature of the wing = (1/A \S aT aA 

T wes = temperature at a point at 7 = 6.0 immediately 
following the acceleration 

a = linear thermal coefficient of expansion, in. per in. 
“Fr. 

a’ = (1 + Via 

V2 = (Q? Ox?) ao (oO? Oy?) 


= Poisson’s ratio 
= normal stress, lbs. /in.? 


v 
o 
> = (T — 7;)/(T mec. 


7;), stress parameter 


= summation of terms from 0 to n 


0 
T = shear stress, lbs. /in.? 
0, = time from M = 1.4, sec. 
Omax. = total time from VW = 1.4 to 6.0, sec. 
6 = 0 O max. 
INTRODUCTION 


Fenn TRENDS IN AIRCRAFT and guided-missile 
flight indicate an ever increasing interest in high- 
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speed performance. At extremely high flight speeds, 
frictional effects in the boundary layer produce ex- 
tremely high temperatures and temperature gradients 
within the aircraft structure. In order to appraise 
the severity of the problem, an investigation of the re- 
sulting thermal stresses was undertaken. Although a 
specific example is presented, the method of attack used 
is applicable to a wide variety of such problems. The 
specific case treated here is a diamond airfoil for which 
the time-space temperature distribution is known.! 
The numerical procedure utilized, although cumber- 
some, readily yields interpretable results and can be 
applied to bodies with boundary curves that are ill- 
suited for analytic analysis. No attempt was made 
to consider plastic flow in the wing, although the 
stresses computed in some portions of the wing on an 
assumption of perfect elasticity are beyond the elastic 
limit of any suitable material. 


EXAMPLE TREATED 


The specific case considered was that of a solid sym- 
metrical double-wedge wing shown in Fig. 1 of the fol- 
lowing size: 


Chord length = 5.5 ft. 
5S per cent 


and thickness/chord ratio = 


The following assumptions are made concerning aero- 
dynamic, thermal, and structural conditions: 


(1) Acceleration is uniform and at a rate of one 
Mach Number per 30 sec. from M/ = 1.4 to 6.0. 


(2) The supersonic flow of air by the wing is as- 
sumed to be two-dimensional, independent of the span. 
Fig. 2 shows the flow pattern around the wing, with 
shock waves at the leading and trailing edges and an 
expansion wave at the mid-chord. 

(3) The wing is made of steel with constant values 
of the thermal conductivity, specific heat, density, 








z 
Fic. 1. Symmetrical double-wedge-shaped wing. 
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Fic. 2. Airfoil and flow pattern of wedge-shaped wing. 


modulus of elasticity, thermal expansion coefficient, 
and Poisson’s ratio. 

(4) Spanwise heat transfer is negligible. 

(5) Heat transfer by radiation to or from the wing 
surface is neglected. 

(6) The flight of the wing is at a fixed altitude of 
50,000 ft. 

(7) The flight of the wing begins with a uniform wing 
temperature of 531.6°F. absolute at a Mach Number 
of 1.4. 

(8) The root of the wing is rigidly built into the 
fuselage so that the root section is at all times plane. 
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Fic. 3. Notation for increment procedure. 
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Fic. 4. Definition of coordinates for interior points in the wing. 
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(9) The span of the wing is reasonably large with 
respect to the chord. 

Assumptions (8) and (9) permit computation of ther- 
mal stresses at the root of the wing which correspond 
to those at a large distance from the free end of a semi- 
infinite wing. 


TEMPERATURE DISTRIBUTION 


The resulting temperature distribution immediately 
after acceleration to a Mach Number of 6 is contained 
in reference 1 and is reproduced here in Chart A of Fig, 
5. The solution for the two-dimensional temperature 
distribution was accomplished by a stepwise procedure. 
Attempts to obtain an exact analytic solution failed. 
It was found that a good empirical approximation of 
the temperature distribution could be obtained by 
means of a product series that converged slowly; hence, 
an analytic solution for the thermal stress distribution 
became impractical. 


ANALYSIS 


Consider a section of the infinite wing with axes de- 
fined as shown in Fig. 1. Although the temperature dis- 
tribution exists only in numerical form, 
the function 7(x, y) as representing that distribution. 
It will be treated mathematically as a continuous func- 
tion with continuous first and derivatives. 
Since the problem under consideration is one of plane 
strain, the following quantities are zero, €,, Tr2, Ty2, and 
the equilibrium and Hookian equations become” 


let us consider 


second 


(0c,/O0x) + (07;,/0v) = O 
(Oc,/Ov) + (O77,/0x) = 0 
Oa,/0z = 0 1) 
er = (1/E’) (0, — v’'o,) + a’/T 
ey = (1/E’) (0, — v’orz) + a’T | 
0 = o, — v(o, + o,) + EaT 
Utilizing the condition of derivative continuity the 


following compatability equation is obtained: 


07(o, + «a,) 07(o, + «a,) _ 
ox? oy? 7 


V2(o, + o,) = —a’E'V?T (2) 


Combining this with the third of the Hookian equations 
results in the following equation for the stresses in the 
z direction: 

— v'a'E')\V?T = KV?T (3) 


V°c, = (—aE 


where 
K = —(aE + v’'a’E’) 


The solution of this equation has two parts; the first, 


or particular, solution is simply 
= KT (4) 


The second, or homogeneous, solution is one that satis- 
fies the equation 
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figure is concluded on page 118.) 
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Fic. 5 (Concluded). 


V7e, = 0 (5) 


As a first approximation, the following form was chosen 
as a solution of Eq. (5): 


og, = Co ao Cix + Cry (6) 


Thus we have, as a complete solution to the stress 


equation, 
0, = KT oS Co + Cix + Coy (7) 


Before determining the coefficients in the solution, one 
more step must be taken. We have, up to this point, 
taken the wing to be infinite. This includes the as- 
sumption that the wing sections remain plane and 
parallel. An examination of the curves of temperature 
distribution shows that there exists a certain amount of 
antisymmetry about both x and y axes. Under actual 
conditions similar to a semi-infinite wing, this will 
result in a curvature of the wing, thus reducing a cer- 
tain amount of the spanwise stress. If the infinite wing 
of the original assumption were to be cut, this curva- 
ture would result. In this case the change in the stress 
distribution would be numerically equal to the stress 
created by the application to the wing of a bending 
moment great enough to produce the same curvature. 
When the infinite beam is cut, there is an accompanying 
phenomenon that must be taken into consideration. 
Since there no longer exists a spanwise restraint on the 
now free end of the wing, the section will no longer be 
plane. However, as previously mentioned, the con- 
ditions at the end of the wing are of minor importance; 
therefore, we are interested only in the effect of the free 
end on the root of the wing. If the assumptions pre- 


viously stated are valid, the effect on the stresses at the 
root of the nonplane end will be negligible. 


Returning to the solution of the stress equation, the 
undetermined coefficients may now be evaluated by 
applying the necessary conditions of equilibrium for 
the free wing. These are: 


(1) The resultant spanwise force on any section 


must be zero. 
A a, dA = 0 (8) 
A 


(2) The resultant moment about the x axis must be 


/ yo, dA = 0 (9) 
A 


(3) The resultant moment about the y axis must be 


I xa,dA = 0 
A 


Applying condition No. (1) 


I ¢,dA =0=K / T dA + cy fa + 
A JA A 
Ci f xdA+C, I y dA 
A A 


Since the area is s7.n netric with respect to the x and the 


zero. 


zero. 


(10) 


y axes, we have 


fe dA = fo dA = 0 
A A 
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Hence, 


Co = -K| ( A) f r aA | 
A 


It is seen that the expression (1/A) f° 4T dA is simply 
the area mean temperature, which we shall define as 
Tm; hence, 


m) 


= K(T es T =) + Cx + Cry 
Applying condition No. (2), 


[ dA=0=K f y(T — T,) dA + 
JA 1 


C, J xy dA + Cy f y? dA 
JA A 


Since the x and y axes are also the principle axes of 
inertia, the expression f° 4 xy dA, being the product of 
inertia, is zero. Defining J°4 y*® dA, the moment of 
inertia of the area about the x axis, as /,, we have 


l 
C2 = —K [ WT am fie dA 
i. « A 


In a similar manner condition No. (3) yields, for Ci, 
| e 


CQ, = -—K “ P x(T — 


T,) dA 


and the final expression for the spanwise stress is 


l 
foxe- 
ae a ee 
ee = 
J x(T — Ta) dA | (11) 
| ee 


T,) dA — 
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leading and trailing edges during acceleration from M1 = 1.4 to 
6.0. 


There now remains the numerical evaluation of the 


above equation. 


METHOD OF SOLUTION 


In computing the temperature distribution, the cross- 
sectional area of the wing was subdivided into small 
areas, and the temperature of each of these subdivisions 
was assumed constant over the area. If this assump- 
tion is continued, the computation of stress is made 
relatively simple. Defining 7, as the temperature at 
any particular point and A, as the area corresponding 
to that point, Fig. 3, then the expression for the mean 


temperature becomes 


T, = (1/A) >> T,An (12) 
0 
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Fic 9. Calculations for thermal stress in wedge-shaped wing 
during acceleration from M = 1.4 to 6.0 at the rate of lg. 


Now defining F, = [(T — T»),A»], the final result 
is an expression for o, in the form 


1 } / l } 
a = K | (7 — Ewe) = L, Sala ei ¥ ya | 
0 0 


(13) 


NUMERICAL SOLUTION 
In order to obtain numerical results for the existing 
Since 
data were not available for a material that would serve 


stresses, values for /, v, and a must be chosen. 


satisfactorily at the elevated temperatures computed, 
the constants for steel at room temperature were 
chosen. The resulting stress values are therefore 
rather meaningless on an absolute basis, but they do 
serve as a means of studying the relative distribution 
of the existing stresses. Using the values E = 30 X 
10° Ibs. per sq.in., » = 0.30, and a = 6.5 X 10> in. 
per in.°F., we get 


E’ = 33 X 108; »v’ = 0.428; a’ = 8.45 xX 10-* 
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Hence, 


K = —(aE + v'a’E’) 


— 280 lbs. per in. per in. °F, 


The actual numerical calculations have been carried 
out on Fig. 5. Chart K shows the evaluation of the 
stresses introduced by cutting the infinite wing. As 
determined by the form of the homogeneous solution, 
Eq. (6), to the stress equation, these stresses vary line- 
arly across the wing. It will also be noted that their 
magnitude is about 12 per cent of the stresses given by 
the particular solution. - In view of the small influence 
of the homogeneous solution on the final stress dis- 
tribution, it is logical to assume that the solution 
chosen, although simple, provides adequate accuracy, 
Hence, Chart L, Fig. 5, represents the distribution of 
thermal stresses in the wing at the instant it has reached 
a Mach Number of 6.0. Values from this chart are 
shown graphically in Figs. 10 through 16 for various 
concentric diamonds in the interior of the wing. 

Only o, stresses are shown, because in preliminary 
calculations it was found that the o, and o, are much 
smaller than these. Fig. 16 is a schematic diagram 
showing the stresses at Mach Number 6.0. The dotted 
line represents the region of zero stress. The stresses 
inside the dotted line are tension stresses, and the 
stresses outside the dotted line are compressive stresses. 


STRESSES DURING ACCELERATION 


Inspection of Eq. (13) shows o, to vary nearly pro- 
portionally to (7 — T,,) at any point, since the last 
two terms decrease in magnitude when the first term 
does. Although the actual variation may not be 
linear, the assumption is fairly good, since the last two 
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Stress distribution on diamond with parameter — = 


Fig. 11. 
0.857. 


terms account for only about 10 per cent of the stress. 
As a good approximation, the stress during acceleration 
may then be calculated by the following equation: 


% = (T - Tul [(o2)a1—6-0 € g a T m)u—6-0] (14) 


Plots of temperature as a function of time were made 
for a series of points on the surface of the wing. It was 
found that the dispersion of all the points with excep- 
tion of the leading and trailing edges was less than 
5 per cent. Fig. 6 shows the distribution in the im- 
mediate vicinity of the leading edge. It is seen that 
the final chart is applicable at all points except the 
leading edge and trailing edge. Therefore, one curve 
would be sufficient for the temperature-time function 
on the surface of the wing. To make the curve useful, 
it was reduced to dimensionless parameters, which are 


9 6, actual time 
Ones time from M = 1.4 to 6.0 
DL = (T — 1)/(T naz. — Ts) 
where 
T = actual temperature at a point 
T; = initial temperature, 531.6°F. abs. 
T maz. = temperature at a point at M = 6.0 


Further investigation of the temperature distribu- 
tions showed that, for any diamond concentric with the 
surface, the time-temperature function was the same. 
The parameter £ was used to define a particular diamond 
as follows (Fig. 4): 

. 2d =—-A(|y| + |x| tan B) 

ee l 
The resulting family of curves, showing the variation 
of the parameter >> as a function of @ for various values 
of £, is shown on Fig. 7. 
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Fig. 8 is a plot showing the mean temperature of the 
wing cross section as a function of the time parameter. 
Chart A of Fig. 9 presents values of [o,/(T — T,,)] at 
a Mach Number of 6 immediately after the constant 
This function becomes infinite when 
(T — T,) is zero. The omitted values on Chart A, 
Fig. 12, are in this region. This region is of little in- 
terest, since the extremes of thermal stress occur near 
the center and at the outer edges. 

To obtain the thermal stress at some point in the 
wing at some time during the acceleration, determine 
> from Fig. 7; T, — 7; from Fig. 8; [¢2/(Tmaz. — 
T»)]| at M = 6 from Fig. 9, Chart A; and (Taz. — 
T;) from Fig. 9, Chart B. 

Then get T — 7, from 


. 1 . 
7 oa Fn _ (F sas. —_ 


acceleration. 


T)>> — (Tn — Ti) (15) 


and oa, from Eq. (14). 


SAMPLE CALCULATION 


— 14.143 in., 
For 


Suppose the stress at the point x = 
y = 1.1314 in., and @, = 103 sec. were desired. 
this point the time parameter would be 

0 = 04/O mar. = 103/138 = 0.75 
and the diamond parameter would be 
- 2(|y| + |x| tan 8) _ 2[1.1314 + 14.134(0.08)] _ 
l 5.380 
4.524 
- = 0.840 
5.380 


To find the stress, the following procedure is used: 
The parameters @ and é are taken to the Fig. 7, and the 
constant > = 0.45isfound. From Fig. 8, (7, — 7;) = 
152°F.; from Chart B, Fig. 9, (Tnaz. — Ti) = 456.7°F. 

From Eq. (15), 


T — Tm = (456.7) (0.45) — 152 = 54 
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From Chart A, Fig. 9, 
lo, CF ene. bat Tm) lar =6 = —106.1 


therefore, from Eq. (14), 
o, = (—106.1) (54) = 5,730 Ibs. per sq.in. 


At M = 6 the stress at this point is 12,600 lbs. per 
sq.in. as read from Chart L of Fig. 5. 
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DISCUSSION OF ASSUMPTIONS AND METHOD 


Under actual conditions the wing might be tapered 
from the fuselage to the tip, both in thickness and 
chord. It will be recalled that the temperature dis. 
tribution was computed on the basis of an infinite wing 
having the section of the root of the actual wing; there. 
fore, the temperatures as computed are probably valid 
only in the immediate vicinity of the root. Because of 
the continual decrease in cross-sectional area as the 
tapered wing tip is approached, it is a logical assump- 
tion that the heating of the wing will occur more and 
more uniformly, thus reducing the gradients. Since 
the stresses and temperature gradients are closely re- 
lated, it is seen that the maximum stresses will occur at 
the root of the wing. It was thus decided that a 
method of computation that would give the root 
stresses would be sufficient. These assumptions per- 
mitted the computations to be made for a semi-infinite 
wing with a section identical to that at the root of the 
wing. This assumption also serves to reduce the stress 
problem from a three-dimensional problem to one of 
plane strain. 

The problem of determining the stress distribution 
throughout the wing as a function of time was resolved 
into two sections for simplicity. First, the stresses 
corresponding to the temperature distribution at the 
completion of were computed. Then, 
rather than repeat the necessary computations for all 
time intervals from JJ = 1.4 to 6.0, a parametric 
method, based on the temperature distribution, was 
evolved. From Eq. (13) it is seen that the stresses are 
nearly proportional to the difference between the ac- 


acceleration 


tual temperature and the mean temperature at any 
given time, the digression from proportionality being 
caused by the nonlinear nature of the antisymmetry of 
the temperature distribution. However, since the re- 
laxation stresses are, in themselves, of small magni- 
tude, the variations in these stresses will be negligible, 
and the assumption of a proportional relation between 
the stresses and the temperature difference is a valid 
one. Unfortunately, the use of the parametric method 
described involves the distribution of four variables 
the two-dimensional location of a point in the wing, the 
value of the temperature difference at that point, and 
the time from \/ = 1.4. In order to reduce the com- 
putation of stresses to its simplest possible terms, the 
construction of a family of curves was necessary; how- 
ever, since four variables existed, it was necessary 
to determine a constant relationship between two of 
them. It was for these reasons that the seemingly 
complex procedure set forth was derived. 

At the present time there exists no information on 
the possible materials for the construction of a wing 
Therefore, the computa- 
In order 


such as has been analyzed. 
tion of actual stresses is out of the question. 
to obtain a concept of the relative magnitudes of the 
stresses involved, sample calculations were made using 
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the properties of steel at room temperatures. An entire 
calculation was carried out on this basis as a demon- 
stration of the method of application of the computa- 
tion procedure. The numerical values presented are 
meaningless and should be considered only as repre- 
senting the range and distribution of stress in the 
wing. The calculation procedure as described is appli- 
cable to any material and temperature distribution, pro- 
vided the stated assumptions may be made. 


DISCUSSION OF RESULTS 


Assuming that the stresses as found are valid for the 
conditions, an analysis of the stress results is in order. 
Of primary importance are the leading and trailing 
edges. It will be noticed that the values of stresses 
presented for these locations are beyond the range of 
credibility. This is a direct result of the assumptions 
of elasticity. Obviously, plastic flow would have 
taken place long before the stresses computed were 
reached. The existence of this flow is, in effect, a stress 
relief. When the condition of high temperature gra- 
dient has been removed, either by stabilization of speed 
at M = 6.0 or by return to a lower speed, the com- 
pressive flow that took place will reverse sign and be- 
come a residual tensile stress. The following cycles of 
stress will then be forced to overcome this tension before 
the compressive stresses can be produced. The com- 
bination of the high stress gradient and the thinness 
of the section at the leading and trailing edges would 
result in a buckling of the edges, since the ability of a 
thin section of that nature to sustain a high compressive 
load is limited. 

To validate the buckling of the edges, a prototype of 
the first 4 in. of the leading edge was made, and a com- 
pressive stress was induced by bending. The addi- 
tion of a heavy base moved the neutral axis of the 
wedge below the wing section so that the whole section 
was subject to a linearly increasing compressive 
stress. It was found that buckling occurred along the 
edge simultaneously with the commencement of plastic 
flow. Waves were formed having a length of a little 
over 1 in., and, as the loading was increased, the buck- 
ling waves increased in amplitude and extended fur- 
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of stress in the wing at VM = 6. 


ther from the edge of the section. Because of the 


symmetry of the buckling pattern, it was felt that its 
effect on the aerodynamic conditions would be small. 
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Buckling Load of a Stepped Column a 


FREDERIC M. HOBLIT* 
Lockheed Aircraft Corporation 


ABSTRACT of applying the resulting equations to obtain the critical 
A procedure is developed for computing the critical buckling load are given. 
load of a column having a finite number of steps with constant The present method applies not only to columns ) 
EI within each step. The procedure is applicable to columns having any combination of pinned, fixed, and elastj- 
with any end conditions and with loads applied within the span, cally restrained ends but also to columns having axial ¥, 


as well as at the ends. : sane 
‘ loads applied within the span, as well as at the ends. 


The column to be discussed is shown in Fig. 1; the 

















NOTATION ’ : : : ‘ 
Q;’s may act in either direction. 
The ith station is at the right-hand, or ‘“‘far,”’ end of the ith I 
segment. DERIVATION OF EQUATIONS 
L = length of the column 
\; = length of the ith segment Considering any point along the segment shown in 4 
E; = modulus of elasticity of the ith segment Fig. 2, isolating the portion to the left of this point, and 
fy ° = UREEOE 98 Senet af the a engenent : taking moments about the point, one obtains the fol- 
P; = axial compressive load carried within the ith segment ; ; ’ - ‘ 
(Fig. 2) lowing linear differential equation of the deflection 
rve: : 
Bi = V Pi/ Eid: awe anc 
x, y = coordinates shown in Fig. 2 (x is measured from the dty 
beginning of each segment) Edi = —Pilv — Yin) + Mis + Sx (1) V, 
9:1 = deflection at the ith station dx” ag 
yi’ = slope at the 7th station ‘ ‘ F ‘ ' , : 
‘M; © indie memeat ot Ge th ation The solution of this differential equation, with 
M; = é Ste < ¥ o,.° , Pa ss 
Q; = external axial load applied to the column at the ith boundary conditions y = y;_; and y’ = y,_1' at x = 0, 
station. (When load is applied only at the ends, is 
Q= P = P;) ] 
S = shear force in the column (Fig. 2) os ee Minn = Mi cos Bx 4 ent 
w = spring constant, moment per unit angle, for the end J hata P, P, ey « 
restraints or for a rotational spring joint in the js S/P,) Sx 
column yi-1 =~ f-*' sin By x + oe (2) 
s s ” i- \< 
0 = subscript denoting left, or ‘‘near,’’ end of the column B; P, 
n = number of segments and, hence, as a subscript denotes ; — tio 
the right-hand, or ‘‘far,’’ end of the column where 6; = VP, Eli. me 
INTRODUCTION 
IT REFERENCE 1, THOMSON DEVELOPS a method for a ' ae 
obtaining the critical buckling load of a column i 
; - ; : ee e} 
having a finite number of steps with different E/ in the - | 
various steps. Thomson's solution, however, does | os SS | 
. z | ; | 
not include the effect of the shear force that is present ° i-t i — STATIONS—n-I n 
ia . _-_—’ —_ eee’ 
in a single-span column whenever the two end moments i I SEGMENTS A 
are unequal and, hence, is applicable only when the con- ke elie RO AOR | 
genan . ’ fide: dapat ppice , Fic. 1. Single span column of » segments. eq 
ditions of end restraint are such as to make this shear ca: 
force zero. It is not applicable, for example, to the Su 
column fixed at one end and simply supported at the fu 
other or to a column fixed at both ends unless the f | x 
member is symmetrical. It is, of course, applicable to Vier y (a 
a column simply supported at both ends and to a canti- | P; | yi at 
lever column. S M 
In the present paper, the basic equations of reference Mit S i (4 
1 are rederived to include this shear force, and methods 
a h M; 
Received April 24, 1950. i 
* Stress Engineer. Fic. 2. Free body diagram of ith segment. 
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STATION i-lji 


Fic. 3. Rotational hinge joint in column. 


Setting x = A; gives y = yi: 


Mi1 Mir 
n= Ya + — p08 Bide + 
Yea’ — (5/Ps) sin BA; + oe (3) 
B; P, 


Differentiating Eq. (2) and again setting x = , gives 
1 = Meg sin Bet 
v4 P i ii 


( P >) iii $ (4) 
Vin ppm +5 


and 


9 


d*y 
M;, = Eid; — = M, 1 COS Bid; - 
1 


dx? 
5) ae 
43’ —- sin BA; (i 
Vii P, B, S1 oO 


If for any segment the axial load is zero, the differ- 
ential equation becomes 


EI (d*y/dx?) = Mis + Sx (6) 
Proceeding as before, one obtains the following equa- 


tions, which replace Eqs. (3), (4), and (5) for the seg- 
ment in question : 


_ + \ 4 Mid? ‘: SAi3 (3: ) 

a a a a 

, 4 Med, SA sins 

F = Vi zs : a 
eS di 

M, = Mia + Sr (5a) 


If within any segment the axial load is tension, the 
equations given do not apply, but equations for this 
case can be derived similarly to Eqs. (3), (4), and (5). 
Such equations will be found to contain hyperbolic 
functions. 

If a rotational spring joint is inserted in the column 
(a condition frequently approached in hydraulic actu- 
ators), as shown in Fig. 3, with a spring constant 
M/@ = w;, the following equations replace Eqs. (3), 
(4), and (5) for the joint: 


% = Vi (3b) 


yi = ver’ + (Mi1/a)) (4b) 
M, = Meir (5b) 


LOAD OF A STEPPED COLUMN 


—_ 
bo 
uo 


GENERAL PROCEDURE FOR COMPUTATIONS 


Eqs. (3), (4), and (5) give y, y’, and M at the far end 
of any segment in terms of their values at the near end. 

In using these equations to determine the critical load 
of a column, a trial-and-error procedure is followed. 
First, a value is chosen for the end load, Q; or if loads 
are applied within the span, as well as at the ends, a 
set of values is chosen for all the Q,’s. As successive 
trials are taken, all the Q,’s should ordinarily be varied 
in the same ratio. In the following discussion, all refer- 
ences to a value of Q apply equally to a set of values of 
the Q,’s. 

Values are then chosen for y, y’, and M at the near end 
of the column. The chosen values must, of course, be 
consistent with the end restraints (for a pinned end 
yo = Mo = 0, for a fixed end yo = ’ = O, and for 
elastic restraint yo = 0 and Mo/yo’ = wo). One of the 
three quantities yo, yo’, and Me is always arbitrary and 
may be given any convenient value. 

Using Eqs. (3), (4), and (5), one then proceeds seg- 
ment by segment along the column until y, y’, and M 
are obtained at the farend. Tabular forms can readily 
be devised to facilitate this numerical work. 

If the values obtained for y,, y,’, and M,, satisfy the 
end conditions at the far end, the column is in equilib- 
rium in its deflected position under the action of the 
chosen end load, Q, and, consequently, this value is the 
critical load. 

If the values obtained for y,, y,’, and MW, do not sat- 
isfy the end conditions, then the trial value of Q was not 
correct and additional trials must be made. However, 
for any trial, the values of y,, y,’, and M, which are 
obtained may be regarded as defining some new col- 
umn, with different end conditions, which does have 
the trial value of Q as its critical load. Usually, as 
discussed in the section on “Special Cases,”’ it is clear 
from the values of y,, y,’, and M,, obtained whether the 
new column is stiffer or less stiff than the given column 
and, hence, whether Q has been chosen too high or too 
low. 

Varying the value assigned to whichever of the three 
starting quantities yo, yo’, and My is arbitrary will, of 
course, have no effect upon the value obtained for 
critical load. If a column is in equilibrium at one 
small deflection, it will be in equilibrium at any other, 
and the magnitude of the arbitrary starting quantity 
simply determines the magnitude of the deflection. 

For each trial, all the constants in Eqs. (3), (4), and 
(5) are known except S, which requires special treat- 
ment depending upon the end conditions. Specific 
procedures for evaluating S and for determining from 
Vn» Yn’, and M, whether the trial value of Q is less than, 
equal to, or greater than the critical load are given in 
the section on “‘Special Cases.”’ 

It is always desirable to sketch the deflection curve 
for each trial value of end load, both to catch obvious 
errors and to assure that the critical load obtained is 
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the smallest one rather than a higher one associated 
with an unstable, multiple wave buckle form. If 
Yn» Yn’, and M, have been obtained by proceeding seg- 
ment by segment along the column as described above 
(rather than directly in terms of yo, yo’, and My as by a 
matrix solution), then all the y,’s, y;’’s, and M,’s are 
known. The deflection curve can therefore be drawn 
quickly and with good accuracy simply by spotting in 
the location of the curve (y;) at each station, drawing 
its tangent (slope y,’) at each station, and sketching a 
smooth curve between. 

The method presented herein is applicable in the plas- 
tic, as well as the elastic, range if E; is always taken as 
the tangent modulus corresponding to the axial load 
P,in the segment.’ 


SPECIAL CASES 


Fig. 4a, Span with Both Ends Pinned, End Loads 
Only.—S = 0. Since for this case Mi. = —Pyi-1, 
Eqs. (3) and (4) can be simplified to 


(3c) 
(4c) 


Yi = Vis cos BA; + (Vi-r’/Bs) sin BA, 


yi’ = —yi-18; sin BAY + Vir’ cos Bid; 


Eq. (5) is then not needed. 

Start with yo = 0; yo’ = arbitrary. For convenience, 
take yo’ = 8). 

The criterion for critical load is y, = 0. If y, >0,a 
longer column would be in equilibrium under the same 
load, so the chosen P is too low. If y, < 0, a shorter 
column would be in equilibrium, so the chosen P is too 
high. 

Fig. 4b, Span with Left End Fixed or Elastically Re- 
strained, Right End Pinned, End Loads Only.—Start 
with yo) = 0; Mo = arbitrary; yo’ = 0 if left end is 
fixed or yo’ = Mo/w» if left end is elastically restrained. 

S = —M)/L in Egs. (3), (4), and (5). 

The criterion for critical load is y, = 0 or M, = 0. 

If for convenience the arbitrary value of Mo is taken 


as My = P, then S = — P/Land Egg. (3), (4), and (5) 
become 
a 
MW = Hi P P cos BA; 
arid LAL)... dA, 
ws —_— sin B.A; — L (3d) 
= BO6 sin BA. + +4) BM = 
ye. = P i SIN PA; ¥i-1 L COS Pj; L 
(4d) 


M; 


M cos Biv; — (v0 + ) i r ( 1) 
4 = 0: os i e 
i—1l it i—1 rd Tonle i IC 


Fig. 4c, Span with Any End Restraints and Several 
Loads.—Start with yo = 0. When left end is fixed or 
has large elastic restraint, My = arbitrary, yo’ = 0 if 


AERONAUTICAL SCIENCES—FEBRUARY, 1951 
a Pp? | OSW,so O< uA <o 
We, p G J 
p~e__ “t= 
(a) | (b) (c) 
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al - 
(4) (e) (f) 





Fic. 4. Various cases of end restraint and loading. 


left end is fixed, and yo’ = Mo/w if left end is elasti- 
cally restrained. When left end is pinned or has small 
elastic restraint, yo’ = arbitrary, My = 0 if left end is 
pinned, and My = yo'wo if left end is elastically re. 
strained. For convenience, take arbitrary My = any 
Q or P, arbitrary yo’ = 6;. 

Carry S in literal form through the tabulated cal- 
culations. Various quantities in the equations will 
then appear in the tabulation in the form (a + )S), 
where a and 6 are the numbers tabulated. 


Use all of Eqs. (8), (4), and (5). Use Eqs. (3a), 
(4a), and (5a) for any segments having P; = 0. 


At the far end of the column, y,, y,’, and M,, are all 
functions of S. Obtain S from the condition y, = 0, 
and obtain y,,’ and M,,. 

The criterion for critical loading is — M,,/y,’ 
It should be noted that the value obtained for (— M,+ 
yn’) is the restraint that would be required to make the 
chosen value of Q the actual critical load. In other 
words, as various values of Q are tried, each trial results 
in finding some possible physical column for which the 
chosen value of Q is the eritical load. If for any trial 
the actual end restraint is lower than the required, the 
column clearly is unstable, and the trial value of Q is 
too high. Conversely, if the actual end restraint is 
greater than the required, the column does not buckle, 
and the trial Q is too low. 


= Wy. 


Negative values of end restraint mean simply that the 
end moment ./,, applied by the restraint acts to increase 
rather than decrease the end slope y,’ and could be 
provided physically by continuity with another col- 
umn less stable than the one in question. 

Fig. 4d, Symmetrical Span—S = 0 in Eggs. (3), 
(4), and (5). Start with the same values of yo, yo’; 
and M, as for Fig. 4c. Treat only half the member. 
The criterion for critical load is, then, y,’ = 0. 


Fig. 4e, Cantilever, End Loads Only.—S = 0. Start 
with yo = 0; Mo = arbitrary; yo’ = O if left end is fixed 
or Yo’ = Mo/w» if left end is elastically restrained. 

Since for this case 

M; = PO ne: _ Vi) a P[(M)/P) oe 1] 


(Continued on page 138) 
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On Some Special Problems in Linearized 
Axially Symmetric Flow 


BERTRAND DES CLERS* ann CHIEH-CHIEN CHANGT 
The Johns Hopkins Unwersity 


SUMMARY 


In view of the considerable amount of material published on 
the linearized treatment of axially symmetric flow, both on the 
subsonic and the supersonic régime, and because of the com- 
plexity of the results, it was thought advisable to work out an 
example that, by its very simplicity, would give a basic ‘‘feel- 
ing” for the behavior of bodies in axially symmetric flow. 

An infinitely long body with sinusoidal corrugations of the 
type shown in Fig. 1 was considered, and the first part of this 
work deals with the study of this body immersed in both purely 
subsonic and purely supersonic flow, whether in or out of wind 
tunnels, etc. The second part of the work considers the inter- 
action of concentric annular regions of subsonic and supersonic 
fow around the same body. The case of a subsonic region is in- 
vestigated in detail because of its resemblance to the supersonic 
boundary layer. Whenever deemed necessary, the results were 
compared or extended to the two-dimensional case by letting the 
radius of the body become infinite. 

Some interesting results are obtained, such as wind-tunnel 


corrections, deviations from the 1/\/1 — M? two-dimensional 
Prandtl-Glauert compressibility correction, and the fact that the 
pressure distribution of an axially symmetric body in supersonic 
flight is no longer proportional to the local slope of the body but 
lags behind it by an angle that is a function of free-stream Mach 
This lag has already been observed 
in particular problems worked out by the Fourier integral 
method, such as the minimum drag bedy problem of thin bodies 
of revolution. It is also found that the wave drag is finite at a 
Mach Number of unity, as was deduced by von Karman in 1935? 
for a similar body. The boundary-layer analogy also shows a 
few interesting results and, particularly, the effect of the thickness 
of the subsonic region on the wave drag and pressure distribution 


Number and body radius. 


of the body. 

It is shown briefly in an example how the whole procedure 
can be extended to a more general axially symmetric body and to 
any number of concentric annular flows, whether supersonic or 


subsonic. 
List OF SYMBOLS 

(x) = equation of the body 

r(x) = equation of the gaseous interface 

r = mean radius of the body 

1 = mean radius of the gaseous interface 

€ = amplitude of bumps 

L = wave length of bumps 

a = 2r/L 

h = distance to outer wall (radius of wind tunnel) 
M = U/a = free-stream Mach Number 

8B = Vi-M 

B = VM?—-1 

r’ rBi; fo’ = ToRi:- h’ = hp; 2’ = 2B, 
af = 1Be; ro" = ToB2; h” = hr; 7” = iPe 
i = 


7 — ry) = boundary-layer thickness 


Received April 28, 1950. 
* Graduate Student, Department of Aeronautics. 
t Associate Professor of Aeronautics. 


aU “ T,(aro’) Ko(at’) 
B = ~ 1 In(at’) + r a 
p2U, Ki (aro’) 
U2B B 1, (aro’) Ki (at’) 
Cc = Ii(at’) - dec - _ ] 
UBe K,(aro’) 
INTRODUCTION 


4 I {HE EQUATION OF MOTION of an isentropic, inviscid, 
irrotational fluid can be linearized if we assume 
that we have a parallel uniform mean flow upon which 


is superimposed a small disturbance velocity g’ such that 


iq’/U| < 1.0 and Ve = q’ = u,'iz + u,'i, where ¢ 
is the disturbance potential. 


In the case of axial symmetry, the differential equa- 
tion governing the motion is, in cylindrical coordi- 
nates, 


d%0 
— + +—-(1-M) =0 | (I) 
Ox 


This equation is linear and will allow us to form a 
general solution, with any particular solutions we can 
find, by simple superposition. If our boundary condi- 
tion is bounded and periodic, it can always be expressed 
as 


_ nas — - 


© 


| n(x) = > fa, sin (nx) + b, cos (nx)] | (2) 


a= 


and it will be sufficient to obtain a solution ¢ 5, OF Yen 
of Eq. (1) for the simple boundary condition a, sin (nx) 
or b, cos (nx), respectively, to build up the solution 
¢r, corresponding to the more complicated boundary 
of Eq. (2): 


ts) 


¢ro = p 2s [sn + Pen] 
n=0 

Therefore, our problem reduces to the finding of a 
solution of Eq. (1) for the simple boundary r, = 7% + 
a sin (nx) or % = % + 6 cos (nx). This has already 
been solved in two dimensions (see reference 5), and the 
results are well known. It is found that, in the case of a 
sinusoidal wall, say, the pressure distribution is also 
sinusoidal, being proportional to the slope—that is, 
90° out of phase—in the supersonic case and in phase 
with the wall in the subsonic case. We shall see that 
this no longer applies in three-dimensional flow. 


127 





128 JOURNAL OF THE AERONAUTICAL SCIENCES—FEBRUARY, 1951 














































WIND-TUNNEL WALL 
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Fic. 1. Picture of body investigated. 
LIMITATIONS OF THE THEORY written CI 
( 
In order to linearize the equations, we assume small (1/U) (O¢/dr),, = dr,/dx Wec 


perturbations, u,’/U, u,’/U < 1.0 and irrotation- 

ality. Moréover, in the process of solving the equa- and if we assume the surface of the body to have the 
tions, we introduce another limitation: Instead of equation 

applying the boundary conditions at the exact surface 

of the body, we apply them at 7) the mean radius.* We ry = % + € sin (ax) whic 
are therefore limited to €/rp5 < 1.0. 











the boundary condition reduces to Press 
ANALYSIS r 
es a W 
Boundary Conditions | (O¢/0r)r, = Uae cos (ax) (3) 
The condition that the flow must follow the surface cay it be 
of the body applies in all cases. This is approximately The other boundary condition will be furnished by inter 
* This is a standard procedure in thin aerofoil theory, where Whatever restrictions we impose on the rest of the flow Fo 
the slope of the surface is given on the chord line. field. 
Subsonic Case: M < 1.0 whic! 
By introducing the new variable r’ = 78, we can reduce Eq. (1) to the Laplace equation Ca 
i Ake et oan ae . 4 
| (0%~/dr’*) + (1/r’) (O¢/Or’) + (0%¢/dx?) = 0 (4) 
and if we assume the solution to be of the form g = R(r’)-X(x), we obtain the solution 
[ = ——— —_— “| The 
| ~ = [acos (kx) + b sin (kx)] [Alo(kr’) + BKo(kr’)] | (5) dime 


where J, and Ky are modified Bessel functions of order zero. 
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by 


low 
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Since we have changed variables, the boundary condition at the wall must be modified accordingly, and Eq. 


(3) becomes 


Bi(O¢/Or’), = Uae cos (ax) (6) 


Ririiaiescdiiesate ten eat _ 


and if we apply it at 7, this allows us to evaluate all the constants in Eq. (5) but one, and we get 


l le E ( ar’) + DKo( =] 
¢g = — cos (ax) : 
Bi Ti( aro’) — DK\( aro’) 


We now have to introduce another boundary condition to evaluate the constant D. We will consider the three 
following choices: 
Choice A, Solid Wall at r = h.—The fact that the flow must be parallel to the wall tells us: 
(0¢/Or), = 0 = (0¢/Or’), 
and if we introduce this condition into Eq. (7) for the potential function, we obtain D = J,(ah’)/K,(ah’) and, re- 


placing in Eq. (7), 


; Ue ; Io(ar’) K,i(ah’) + Ko(ar’)I\(ah’) 
g = cos (ax) ae ; : - (8) 

1 T(aro )Ki(ah") — Ki(aro \Ti(ah ) 
Choice B, Free Jet Boundary atr = h.—In this case, we must have ambient static pressure at the boundary, since 
itcannot support any force; therefore, (0¢/0x), = 0, which allows us to obtain D = —J)(ah’)/Ko(ah’) so that the 


potential function is 


(0) 


Ue eae —s Ko(ar’)Io( ah’) | 
cos (av : : : ; — (9) 
1 T\(aro’) Ko(ah ) + Ki(aro )Lo( ah ) 





Choice C, Free Flight—The only condition that we can impose here is that our solution remain finite at infinity. 


Wecan therefore obtain it by letting h ~ in Eq. (8), say, and we get 


Kol ar’) | ( 10) 
Ky ary’) 


yg = (— Ue/B;) cos (ax) | 


which vanishes when we let r > ©, and this satisfies our boundary condition. 


Pressure Correction Due to Wind-Tunnel Wall Interference 
We now have obtained the velocity potentials, and, since the pressure coefficient is defined as 
Cp = (—2/U) (0¢/dx) 
it becomes a simple matter to determine the pressure distribution at the surface of the body, in which we are most 


interested. 
For free flight, we get, from Eq. (10), 


C (—2ae/f;) sin (ax) Bed (11) 
>. = (—Zae sin (ax > 
PR F. é 1 K,(aro’) 


which we can compare with the wind-tunnel pressure distributions. 
Case A (Closed-Throat Wind Tunnel).—We obtain, from Eq. (8), 


C 2ae . ( | eet + Be 12) 
Neen sin x — . A eal 2 
Pw.t Bi wi I,(aro’)Ky(ah’) — Ki (aro’)Ii(ah’) 


The ratio of Eqs. (11) and (12) has been plotted in Fig. 2, for various values of ary’ and compared with the two- 


dimensional correction corresponding to ary’ = ©, which is given by 


(Coy +. Cry p )2-dim. = coth (ah’) 
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It can be seen that this ratio is always greater than unity; the reason the two-dimensional effect is greater than 
the three-dimensional one for large values of ah’ is because of an area effect. 
Case B (Open-Throat Wind Tunnel).—From Eq. (9), 


To( ary’) Kol ah’) ed Kol ary’) Io =| 
(13) 


2ae 
Seen © sin (ax 
Pw.T. By at on) Feterrean co Ki(aro’)Io(ah’) 


and the ratio of Eqs. (11) and (13) has also been plotted on Fig. 2, showing that the measured pressures will be 
lower than free flight. The corresponding two-dimensional expression (arp’ = ©) is 


(Cewr. ‘Cry » )2-dim. = tanh(ah’) 


It is seen, from the overall nature of Fig. 2, that the wind-tunnel interference becomes negligible for values of 
(ah’ — ary’) exceeding 2. This shows what size the wind tunnel should be in order to approximate free-flight 
conditions; this size is seen to increase for any particular model as the Mach Number approaches unity. 


Compressibility Effects (Free Flight) where Cp; is the pressure coefficient in the incom- 


Considering the wealth of information available on pressible case. In order to obtain Cp, for our axially 


the results of potential theory and the relative ease Symmetric body, we merely set the Mach Number 

with which we can compute the flow of an incompres- / = Oin Eq. (11), and we obtain 

sible, inviscid fluid about any particular body, we 

would like to know what happens to our calculated Cp | eee (15) 

pressure distribution when we put the same body in a Ce, Bi LKy( aro’) Ko( aro) ' 

compressible fluid. The pressure correction has been 

worked out for two-dimensional flow, known as the This equation has been plotted in Fig. 3 for different 

Prandtl-Glauert correction, values of 7)/Z and is compared with the two-dimensional 
expression of Eq. (14) and with the Géthert correction 








(Cp/Cp;)2-aim. = 1/8: (14) 
Cp log, 
saeeeial + one at enc . (5 =1+— Be Bs (16) 
1} Cp; G.c wees a log, (7%o i) 
| a | fs + } } } } } } } 4 4 
\ -_ : ‘ ; , 
=e ae originally obtained by considering the flow about thin 
= = il ellipsoids of revolution. 
| | It can be shown that Eq. (14) is merely a limiting case 
ir Hit idl <a aii: dias “ala ili’ im cian eae’ of Eq. (15) for ary’ > © and that Eq. (16) corresponds 
OS Ee ee oe oe to taking only the linear terms in the series expansions 
i] | | ma for the Bessel functions; this last expression is shown 
CLOSED not to agree too well, however, even when 7/Z is as 
THROAT 
pt a Sa A $0 
| | a | 3 ++ r 7 
' ee © a 
+ + + + + + + 
a 7" HG & 2 1Q7! 
NS - | 4 4 } } } 4 4 i: am . GC L -— 
Ye. ee | + . 
; = yee +(@) > ee 107! 
ars 02 ] ee ne | t] —t— 
— a es a (1) GG —=10 
Z erie. 7 ce tL : 
@{ t1'| 
OPEN THROAT 
10 
2 
oh-(08), ep * 2nVJi- M® nt ——— Fic. 3. Compressibility corrections—curve (2), two-dimen- 
sional—curves (1) and (6), Géthert correction—curves (3), (4), 


Fic. 2. Wind-tunnel corrections. (5), exact three-dimensional expressions. 
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small as 10~*, and to be very far off for a value of ™/L = 107! over all the practical range of Mach 
Numbers. 


Supersonic Case: M > 1.0 


By introducing a new variable, r” = 72, we can transform Eq. (1) into the wave equation 
(0?y/Or”?) + (1/r”) (O¢/Or”) = 0°y/Ox? (17) 
and, again assuming the variables to be separable, we obtain a solution of the form 


| go = [acos (kx) + b sin (kx)] [AJo(kr”) + BYo(kr”)] (18) 


and we must again modify our boundary condition at the body in terms of the new coordinate system, so Eq. (3) 


becomes 
| | 
| Bo(O¢/Or"),, = Uae cos (ax) | (19) 


which can be applied immediately to Eq. (18) in order to evaluate some of the constants; we get 


Ue Ji(ar”) + DY (ar") ; 
g=- cos (ax) = - > (20) 
Bo Ji(aro ) +. Dy 1(ary”) 


which contains only one arbitrary constant to be evaluated by introducing one more boundary condition. We 
will again consider three choices for that extra boundary condition. 
Choice A, Solid Wall atr = h.—We must have the flow parallel to the wall or 


(0¢/Or), = 0 = (O¢/Or”"), 


which gives us D = —J;(ah")/ Yi(ah"), where h” = h®.; if we replace into the potential function, we get 
| le Jy(ar”) Yi(ah") - Yo(ar”) J\(ah") | 9 
eee cos | ax) VAY t4 ” r ” ” (21 ) 
Bo Ji (ary )Vilah ) —_ Yi(aro \Si(ah ) 


Choice B, Free-Jet Boundary at r = h.—In this case, we must specify the pressure at the boundary to be ambient 
static pressure or 


(Oy/Ox), = 0 


which gives us D = — Jy(ah")/ Yo(ah") and the potential function becomes 
Ue Jy(ar”) Vo(ah") — Yo(ar")Jo(ah") pa 
situa Nest cos (ax WW ” , ” ” (22) 
2 J\(aro )Vo(ah ) — J 1( a7 )Jo(ah ) 


Choice C, Free Flight—The only further restriction that we can impose on the flow is to eliminate the disturb- 
ances from infinity—i.e., the incoming set of Mach waves. This is done by considering the Mach waves to become 
very nearly plane waves at a large distance from the body where the potential function will therefore be of the 
form 

. A J ” 
g = F(x — r”) + G(x +7") 
and making G=0. Using the asymptotic expansions of J, and Y,, we can evaluate the constant.* We get, after 
also applying the boundary condition at the body, 
<a Ue j ” ” , ” , ” ” r " 
p= I Ji(are”)Jo(ar”) + YVilaro”) Yo(ar”)] cos (ax) + [Si(aro”) Yo(ar”) — (23) 


Bo [J17( ary”) aa Y,7( ary”) | 
Vi(are”)Jo(ar”)] sin (ax) } 


*The same result can be obtained by taking ¢ = Re| Ae’ Hy (ar") |, which represents the outgoing family of waves,? while 


. » Hy) (ar”) represents the incoming family. 
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Free-Flight Pressure Distribution 


oxi wad saa odin with pennant to x and nee it by —2/U, we obtain 


— 2a 


~ Bal Si? 2( aro") + Y12( aro wy le aro "\Si( ary ”) + 


y = tan 


JOURNAL OF THE AERONAUTICAL SCIENCES- 


a 
lrary” [JoCaro” )Ji(arry” y+ Vo (aro ro") Viana”) 


FEBRUARY, 1951 


wr (ax — y) (24) 


4 
Yo( (aro ") Vi(aro”) |? +- en ra "2 f 





> 


vi (25) 


This shows us that the pressure distribution is again sinusoidal but out of phase by an angle y with the body sur. 


face. Thisimmediately tells us that we will have drag. 


If we compare Eq. (24) with the corresponding two- 
dimensional expression 


(Cp) " = 


2-dim. 


— (2ae/ Be) sin [ax — (r/2)] 


we see that their ratio and also the phase angle y 
contain the different parameters 7, L, and M only in 
the form of the product 


2r(ro, 'L) Be 


” 
a = 


Therefore, their value will be the same for equal values 
of aro”, irrespective of the contributions of the indi- 
vidual parameters. 

In Fig. 4, we have compared the two- and three- 
dimensional expressions, and it is seen that this ratio is 
close to unity from a value of arp” = 4.0 0n. In other 
words, as far as the ompotnte of the pressure coefficient 
is concerned, a body with ary” 2 4.0 is two-dimensional 
to all practical purposes. However, it can be seen from 
Fig. 5 that the phase —_ does not go to the two- 
dimensional value of 7/2 quite so fast. (The scale in 
terms of Mach Number is shown for a body of 7/L = 
1.0.) 

We can investigate what happens at extremely small 
values of ary” by taking only the linear terms in the 
series expansions of the Bessel functions. We find that, 
for ary” < 1.0, 
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Fic. 4. Supersonic pressure distribution. 


—4e 
Cp = 


roe? 


+5 + — [loge (aro”)}? , : sin (ax — yp) 


where 


y = tan! [1/2 log, (aro”) | 


and, similarly, if we investigate the case of ary” > 1.0, 
we find the two-dimensional expression for the Cp with 
y = 7/2 given above. The fact that the pressure dis- 
tribution lagged 90° behind the wall corrugations and 
had an amplitude proportional to ae (or the slope) in this 
case was already known. 


Free-Flight Wave Drag 


Taking a cylindrical control surface around the body 
and choosing it of radius 7 in order to use the ex- 
already obtained for the derivatives, we 


"L se Pa) 
D = — 2xnp J ( 4 ( °) dx 
0 Ox! ,,\OF / x 


which gives us 


pressions 
have 


D = Aap U*e?/ Bo? [J,?(aro”) + Y17(aro”) | 


and, if we introduce the drag coefficient, we get 





y 
aa a = Ee —_— 
= 
y 
, 
y 
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Fic. 5. Phase lag of supersonic pressure distribution 
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4 («/L)* (L/ro) (26) 


Cp = 9 W , 9 ” 7} 
? Be” [J1?( aro ) + Y,*(aro )] 


This equation has been plotted in Fig. 6a for a body of 
n/L = 1.0, «/L = 10~*, as well as the two-dimensional 


drag coefficient 


(Cp) = 41*(e/L)*/Be 


2-dim. 


for the same €/L of 10~°. 

One can see clearly that, although the two-dimen- 
sional drag coefficient becomes infinite at a Mach 
Number of unity, the three-dimensional one has a defi- 


nite value of 


(Co) yoyo = 40 *e/L)*(10/L) 


as had been obtained by von Karman in his 1935 paper.’ 
This finite value is explained by the fact that, although 
the Cp becomes infinite, the pressure distribution 
is exactly in phase with the wall as in the subsonic 
case. 

We have also plotted the ratio Co/ Cody aim. in Fig. 6b, 
which is again a function only of ary”. It can be shown 
to have a slope of r*/4 at the origin. 


Example 


Let us consider a boundary 


By + B,, cos (nx) 


r(x) = 
n=1 
where 
Be = (1/m) [sin 6 + (x — @) cos 6] 
B. = l |= (n — 1)@ sin (n+ aad 
"ne n— | n+ | 


which is a body of the type of Fig. 7 with parabolic 
bumps. 


(A) Subsonic Case.—Let us try to find a solution of 
the Laplace equation in r’ = r8, which will satisfy the 
above boundary conditions. We will take a solution of 


the form 


ie 3B [a, cos (kx) + & sin (kx)] X 
k 
[A,Jo(kr’) 4 B,K, (kr’) | 


and the boundary conditions are: (a) 


(3:1/U) (O¢/dr’),, = dr,/dx 


b 


flow must follow wall. We assume @ small and evalu- 
ateatr = 1.0. (b) Solution must remain finite at in- 


finity, from which we find 
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U <= Ko (nr’) 
>= B,| = > x) 
Y 8; > Bal sin (x 


from which we can find the pressure coefficient 


o = Ko(nB,) 
Cp=-— > nB,, | = cos (mx) 
1 n=l K,(")) 


This shows that we have merely introduced a weighing 
function [Ko(n8,)/K,(n6;)] into the equation of the 
boundary. 
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(B) Supersonic Case-—When M > 1.0, we have a solution of the form 
g = >> larJo(kr”) + be Yo(kr”)] cos (kx) + [ax Yo(kr”) — by Jo(kr”)] sin (kx) 
k 


where we have already introduced the fact that we only allow outgoing waves. 


Our other boundary condition js 


(B2/U) (O¢/Or"),, = dr,/dx 


flow must be parallel to the wall, which we evaluate at ry = 1.0, for small 0. 


_UéS Bn j 
” Bo Py J\?(nBe) + Y2(mBo) \ 


and we can find the pressure distribution at the body 


2 ge nB » 


C= -— — 
Be x Ji?(nB2) + Vi2(n 


3 4 [ V¥i(mB2) Yo(nBe) + Ji(nB2)Jo(nBe)] cos (nx) — 


We then get 


[Vi(mB2)Jo(nr”) — Ji(nBe) Yo(nr"”)] cos (nx) + 


[Vi(mBz2) Yo(nr”) + Ji(nB2)Jo(nr”)] sin (nx) 


9) | 
sin (mx) { 
TN Pe 


It is therefore a straightforward procedure to calculate the complete flow configuration for any such complicated 


boundary. 


MIXED FLow 


Let us consider the problem of having two concentric 
flows of different Mach Numbers, densities, velocities, 
etc., separated by a discontinuity surface or gaseous 
interface 7;(x) and with two arbitrary boundaries 7,(x) 
and 7,(x). Then, we still will have the differential 
equation 


dX 1d d? 
= + (1+ M) 2 =0 
x 


| Or ror 
| 


which will be satisfied by both yg, and yg: and for which 
we have the four following boundary conditions: 


(1/U,) (O¢./Or),, = dr,(x)/dx 
flow must be parallel to r,(x). 
—p,U,(0¢;/0x),, = —p2U2(0¢2/dx),, 
same pressure on both sides of the interface. 
(1/U;) (O¢:/0r),, = (1/U2) (O¢2/0r),, 
same flow inclination on both sides of the interface. 
(1/U2) (Og2/Or),, = drn(x)/dx 


flow must be parallel to 7,(x); theoretically, this allows 
us to solve for both g, and g, thus characterizing 
everything about the flow on both sides of the inter- 
face r;. It is clear that we can also solve this problem 
for any number of discontinuity surfaces. 

We will only solve a particular case here which was 
thought to be most important because of its interest as 


We know that the potential functions are of the kind 


| g1 = [ailo(ar’) + a2Ko(ar’)] cos (ax) + [aslo(ar’) + asKo(ar’)] sin (ax) 


a first approximation to the supersonic boundary-layer 
theory as represented by a step function velocity pro- 
file. 


ANALYSIS—FREE FLIGHT—MIXED FLOW 


We will limit the present investigation to the par- 
ticular case where the equation of the body is 


Yr, = ro + €sin (ax) 


as in the first section, where the region close to the 
body and up to 7; is subsonic at a Mach Number 
M, < 1.0, pressure P;, density p;, temperature 7, and 
for which the flow is supersonic and unbounded for 
r >r;, with flow functions M2 > 1.0, Pe, po, 72. 

The boundary conditions for this case are the follow- 
ing: 

(1) For ¢g, only, 

(1/U,) (O¢;:/Or),, = ae cos (ax) 


flow must follow wall. 
(2) For g, and gp, 


— l (0¢;/0x) 7 = pol ’2(Og2/0x) i 
same pressure on both sides of interface. 
(3) For g, and gs, 
(1/U,) (O¢:/dr); = (1/U2) (O¢ge/dr); 


same flow inclination on both sides. 
(4) For g only, we must only have the diverging 
family of waves—i.e., no disturbances coming from in- 


finity. 


(28) 
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and 


go = [bJo(ar”) + be Yo(ar")] cos (ax) 


AXIALLY 


SYMMETRIC FLOW 


+ [bi Yo(ar") — bJo(ar”)] sin (ax) (29) 


applying condition (4) to g2, as we had done in the first part of this work. After some labor, we get 


ous I he ] f U2p; z ; is oe 
a 9) : a —~—— Ky [C(Jo®? + Yo?) + B(Joi + Yo¥s)) — 
, BiKi(are”) (Bs, + CJo)? + (BY; + CYo)? Ube 1 | ) 0/1 1) 
U 
7 KolBU* + V2) + Ci + YoY)If 
where the J, and Y, are evaluated at ai” = (22/L)i®. and the K, at ai’ = (2x/L)i@, and similarly for the fol- 


lowing equations. 
Band C have the following meaning: 


piv, . 
=> To( 4 
. p2U» | sili K,(aro’) 


we also get 


sik — Use : aii : U2B, ; 2 
; BKy(aro’) pels U\B2— (at’ 
a = [a1 811;( ary’) - Ure] B,Ky(aro’) 
ag = a3[I;(aro’) K,i(aro’) | 


which determines ¢ entirely. 


} = l ie pil; U268; l 

», = ° ‘ 

83:Ky(aro’), pol» U Be (at’ 

te, —l ‘ie : Pl l f : l Bes : l 
, B,K,(aro’) pos UB. (at’ 


Now we can determine everything about both regions of 


Shape of the Gaseous Interface 
The local slope of the interface is given by 


(dr,/dx) = (1 


and differentiating ge, we get 


shee et) 
and 


U2pi 


EG - ee | 
Ui Be 


- Kil ary’) 


l 
2 (BM, + CIo)? + (BY: + CY)? 


We similarly obtain }, and dy. 


BJ, + CJo 

) (BI, + CJ)? + (BY + CY)? 
BY, + CYo 

) (Bs, + CJo)? + (BYi + CY)? 


the flow. 


U2) (Oge Or); 


dr; pil; (€ 1) | J? + ¥;* | . 
= is are : —- i cos (ax — wy) 
dx pol 2 Bi Ky(aro’) (BJ, a CJo)? a (BY, + CY,)? 
where 
—2C l ‘ 
yi = tan- | —- - ~ - i | (30) 
rat" B(J,? + Y;”) + C( Sos; + Yo Vi) 
and, if we integrate with respect to x above, we get 
r,=1+e¢ sin (ax — yy) 
where ¢; is the amplitude of the wiggles on the interface 
€; — aU, , ] : | J;? + ic if (31) 
€ p2U» at’ K,;(aro’) (BJ, -+- CJ)? ae (BY, + CY>)? 


which shows us that the interface is a corrugated pipe of mean radius 7, with bumps of an amplitude e; and out of 


phase with the wall bumps by the angle y. 


y, turns out to be negative, which means that the interface lies ahead 
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of the body by that angle. It is seen in Figs. 7 and 8* that, when we shrink the subsonic region to zero, the 
interface coincides with the body and that the amplitude of the wiggles on the interface decreases sharply as we 
increase the thickness of the subsonic region. We have also plotted on Figs. 7 and 8 the corresponding limiting 
case of 7)/L = © —e.g., the two dimensional case—for which we find 





€; 1 | 
= - (39 
(“) Ee aByT + (po? U2”, pil I?) - ( U2?B;? U7 Be”) sinh? = } 32) 
and (W,) 2-dim. = tan! [ — (po U2 ‘py l ,) +( U 2B, il Bo) tanh apy T] | (33) 





where 7 is defined as (i — ro), ‘“‘the boundary-layer” thickness. The amplitude e; is seen to decrease much more 
gradually, starting with zero slope at the origin, while the phase angle rises much more abruptly and later flattens 
out. 


Pressure Distribution at the Body 
Let us define the pressure coefficient at the body as 
Cp = (2/ U2) (Og; Ox) ro 
using U; as our reference velocity. This choice is consistent with our boundary-layer analogy. 


We obtain 


Cp { [a1 om U, aero Ko( aro’) |? + a3} ' . sin (ax = Yo) (34) 


i User Bi K1(aro’) 


where 


é a3 = 
¥2 = tan™! | : = | (39) 
a, — UjyaerKo( aro’) . 


and a, a; are the constants evaluated in the potential and have not been written out here in detail in order to save 
space. 

The two expressions above show us that the pressure distribution is not in phase with the body, which would be 
the case were it in an infinite subsonic flow field, and that we will therefore have drag. This was to be expected, 
and the drag should be just the one obtainable from Eq. (26) for a body of radius 7, bump amplitudes e;, and ina 
free-stream Mach Number M:. Thus, the corrugations of the interface send out waves that will show up as an 
overall loss in energy or drag. We have plotted the above results in Figs. 9 and 10 for the particular body con- 
sidered; as it should be, when we remove the interface further and further away from the body, the value of the 
Cp tends toward that given by Eq. (11) of the body in an infinite subsonic field (that is, 7)/Z = 10~!, «/Z = 107, 
M = ~/2/2) and represented by the dotted line in Fig. 9. Similarly, when we shrink the subsonic region to zero, 
the Cp tends to the value given by Eq. (24) for that same body in an infinite supersonic field of Mach Number? 
equals 7/2. 

When we examine the behavior of ¥2, we see that it tends toward the value given by Eq. (25) when the interface 
coincides with the body and vanishes when the subsonic field becomes large. 

We must remember that the pressure on the body is due (1) to a subsonic pressure, always in phase with the 
wall, and (2) to an out-of-phase pressure due to the influence of the interface. But the pressure must be the same 
on both sides of the interface; therefore, we can consider the additional pressure as the one on a body in purely 
supersonic flow and with corrugations, diameter, and position equivalent to that of the interface. With this pic- 
ture in mind, we see immediately that ly + yo must behave as our y of Eq. (25), since, when the interface moves 


* Figs. 7, etc., apply to a specific body of ro/L = 1071, e/L = 10-2 and to the case where the two flows were chosen with the 
same stagnation temperature, M, = /2/2 and M, = 1/2. The static pressures on both sides of the interface were taken 
equal (P; = P2), and the ratios of the flow functions were determined by using the relations: P = pRr, a? = yRr, 7ro/t = 
{1 + [(y — 1) /2\M?}, with y = 1.40. The same applies to the two-dimensional limiting case of ro/L = ©. The only case in 
which ¢/Z occurs is in Fig. 9. All of the others are valid for any «/L. 
+ Remembering that we must transfer our Cp in terms of go, 
Cp == (2/ l ly) (O¢; /OX) ry = (pol I, /p, U; ) [ = (2/ l 1s(Og2/OX) ro | 

with the help of boundary condition (2), since we must take into account the fact that the Cp given by Eq. (24) is Cp = —(2/U2)X 
(Og» OX) ry. 
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out to infinity, it moves 90° ahead of the body, but the pressure on its outside lags by 90° so that those two effects 
cancel out, and the additional pressure also gets in phase—i.e., 2 — 0. 


Wave Drag 
To calculate the drag, we again take a control surface at the body 


L " ” 9 mG 2a 
p= 2am fo (261) (261) ae = Seematie cost (ax) da) 
0 Ox To or To BLK i( ary’) 0 


and, if we introduce the drag coefficient Cp again, using free-stream conditions U2, p,, where Cp = D/p2U2*-mroL, 
we obtain 
of Ot (e/L)? ] I nee 
c, = 2 _— a —- " - ewe (36) 
p27 l P (4 L)* (Yo i 1B," Bo?K?( aro’) (BJ; + CJo)? + (BY, + CY,)? 

where, here again, the J, and Y, are evaluated at at”. 
In Fig. 12 we have compared this expression with the drag coefficient of the same body in an entirely super- 
sonic field, given by Eq. (26), and have also plotted the limiting case of ro/L, i/L becoming infinite but keeping 


(i —r)/L finite. We obtain for this case 


( Cp ) ws - (37) 
Cog p./ 2-dim ~ cosh? aBiT + (p2?U2?/p:2U42) « (U2? By2/ U2 62?) sinh? afiT 
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(Cp a? = [(€:/€)9 4im.]” 


This two-dimensional expression obviously also has zero slope at the origin. 


The reason for the three-dimensiona] 


drag not vanishing so fast as the two-dimensional one, even though the amplitude of the bumps on its interface jg 
smaller, is that the energy-radiating area increases as h, while it stays constant in the two-dimensional case. 

This comparison shows us that three-dimensional bodies are affected much more rapidly at first by boundary. 
layer growth but that the effect is gradual after that initial plunge. 





a 
L L 
Influence of boundary-layer thickness on wave drag. 


. ish 


Fic. 12. 
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Conical Flow Methods Applied to 
Uniformly Loaded Wings in 
Subsonic Flow 


(Continued from page 114) 


Since the x coordinate has been stretched by the 


amount 1/V1 — M? it is necessary to multiply the 
loading by V1 — M? so that it remains invariant, and 
the slopes must be diminished by the factor Vi - M?; 


therefore Eq. (17) takes the form 


i! on (1 ~ M?) Uo 


wV 


| 20.0" + log |x’? — 1 | + 


constant (49) 
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Flow, 


Buckling Load of a Stepped Column 
(Continued from page 126) 


Eqs. (3) and (4) can be simplified to 


Mo (> 
y= <i 
. P P 


Mo 
y i = S i .)a. sin Pir, + Vi 1" cos Bir, (4e) 


/ 


Vi . 
— Vi ) cos Bij; + * mi sin BA; 


i 


(3e) 


Eq. (5) is then not needed. 
The criterion for critical load is y, = Mo/P. 
If for convenience the arbitrary value of Mp is taken 
as My) = P, then M)/P = 1 in the above equations. 
Fig. 4f, Cantilever, Several Loads.—S = 0 in Eqs. 


(3), (4), and (5). Start with yo = 0; Mo = arbitrary; 
yo’ = 0 if left end is fixed or yo’ = Mo/w» if left end is 
elastically restrained. 

The criterion for critical load is 1, = 0, or Mo = 


Oi. 
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RIEF REPORTS of investigations in the aeronautical sciences 

and discussions of papers published in the JOURNAL are 

presented in this special department. Publication will be com- 

pleted approximately 8 weeks after receipt of the material. The 

Editorial Committee does not hold itself responsible for the opinions 

expressed by the correspondents. Contributions should not exceed 
800 words in length. 


Generalized Complex Forces in Flutter 
Calculations 


J.E. V. Nilsson and N. B. Langefors 
Stress Department, Saab Aircraft Company, Linképing, Swzden 
October 18, 1950 . 


B: THE INTRODUCTION OF generalized coordinates g; in flutter 

calculations, one encounters the problem of determining 
the generalized complex air force P; connected with each gener- 
alized coordinate (degree of freedom, mode of vibration). The 
usual way seems to be to calculate these generalized forces by 
forming the ‘‘virtual work done”’ 


bW = DOP; bq: = D> P2; bz; (1) 

Against this method the objection has been made! that such a 

“work done’’ cannot be formed with the complex air forces. This 

objection seems to be correct, since 6W in Eq. (1) is complex and 
cannot represent a work done in Lagrange’s sense. 

The introduction of generalized coordinates defines, however, 


acertain coordinate transformation 
z= Aq (2) 


from which, in the real case, the transformation rule of the forces 
isobtained by the knowledge of the fact that the force vector mul- 
tiplied by a deformation vector forms an invariant (the work 
done). In a similar way, the complex vectors form a complex 


invariant 


6 = >P;-8q:* (3) 


whete 6g;* is the conjugate vector of 5g;. That this is true is seen 
from the physical fact that the real and imaginary parts of 6W 
form the active and reactive energy. Ina similar manner to the 
real case, the transformation law of the complex fo1ce vector P is 
obtained from the invariant. It is 


P = A,*P, (4) 


where the subscript ¢ indicates that the matrix has been trans- 
posed. 

In flutter calculations, the generalized coordinates usually are 
introduced by a real transformation defining, for example, 
modes of vibration or the vibration amplitudes of a finite number 
of parts of the wing. Complex modes of vibration are then intro- 
duced by adding real forms multiplied by complex coefficients. 
In this case, when z and g are real, the matrix A is real, and, if 
P, is complex, one obtains 
SW = bWe + i8W, = > Pip 52; + 1d, Pojr 82; = 

>> Pirdgi +1 >> Pirdgi (5) 
From this expression it is concluded that, under these conditions, 


the usual method represented by Eq. (1) defines the same trans- 
formation for the real and imaginary parts of P. This is also the 


case in Eq. (4) when A is real. The two methods are thus 
identical, and the usual method therefore gives correct results 
when ¢ and qare real. 
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Vibrations of Rectangular Plates 


Harold Lurie 
The RAND Corporation 
November 6, 1950 


7 NATURAL FREQUENCY OF RECTANGULAR PLATES with 
either clamped or simply supported edges is readily available. 
However, there is little direct information for plates with mixed 
edge conditions. In view of the fact that the vibration and buck- 
ling problems for a plate are identical eigenvalue problems, it 
would seem profitable to utilize for the vibration problem the 
large amount of data that the N.A.C.A. has computed for the 
buckling of rectangular plates. In particular, if two opposite 
edges are simply supported, the analogy is extremely simple, no 
matter what the other two boundary conditions may be. It 
can be easily shown how the fundamental frequency is deter- 
mined as a function of the plate buckling factor k. 

Let the rectangular plate be simply supported along the two 
edges x = 0 and x = a and have various edge conditions along 
the other two edges y = O and y = b. Assuming a sine function 
in the x direction, the differential equation for f(y) for buckling 


is! 
d‘f 2m?x? d*f m'x4 Nz m*x? 
oo See em de Pome f=0 (1) 
dy‘ a? dy? a‘ D @ 
Similarly, the corresponding equation for free vibrations becomes? 


d‘f — 2m*x* d°f + mnt _ a) — om 
dy' a? dy? a‘ D/}- 


The general solution for both Eqs. (1) and (2) can be written as 
f(y) = Gcosh ay + CG sinh ay + C; cos By + CG sin By (3) 


where the corresponding values for a and 6 are: for buckling, 


mtx? [Nz mx? | mx? IN, mx? 


TY ; du, is a* a i a? (4) 


* Ve Vd «@ 


and, for free vibrations, 


lm?e? | mr? | reo? 

“Te *y ,<% a’ WD 
Considering the two edges y = 0 and y = b to have the same 
boundary conditions in the vibration case as in the buckling 
problem, the constants of integration in Eq. (3) will be identical 
in terms of a and 6 for the two problems. Hence, for equiva- 


(5) 


|p? 
D? 


lence, let 


pw?/D = (N,/D)(m*x?/a?) (6) 
The coefficient & for critical buckling is defined as 
N, = k(x?D/b?) (7) 
Substituting Eq. (7) into Eq. (6) yields 
w? = (m*r!D/u) [k/(ab)?) (8) 
Hence, for any boundary conditions along the edges y = 0 and 
y = b, the value of k can be obtained from Timoshenko! or the 


various N.A.C.A. reports,*4 and the natural frequency can be 


computed from Eq. (8). 
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The vibration problem differs in one respect from the buckling °° t/m= 0 ( RECTANGULAR) 








3 1.00 2 mome 
phenomenon—viz., that for the former the fundamental mode The ¢ 
always corresponds to m = 1. Hence, for the fundamental fre- the a 
quency, Eq. (8) becomes \to.8 2 a 

0.96 \ \ wk | value: 

, . 4 —— en — + —_— —__—_—___ o” 

w? = (r!D/p) (k/(ab)?] (9) cP bs 4. O(DELTA) P 0.169, 

. atte Ni | —T rial 

where & is the buckling constant form = 1. This may not neces- cr, ee er | —_ 
sarily correspond to the lowest value of k for a given a/b ratio, | ial | for 

0.92 - 4 accuré 





whereas in the case of buckling k and m are always chosen to give 
the lowest value of k. 
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. FIG.I- CENTER OF PRESSURE FOR CLIPPED DELTA fallac 
WING-BODY COMBINATION WITH NO AFTERBODY airpla 
airpla 
: : , ' P ’ nl 
On Supersonic Wing-Body Center of Pressure tions as compared to the exposed wing alone is uniformly less i 
predi 


than 7 per cent over the entire range of the parameters. 

G. K. Morikawa and T. F. Coleman 
Research and Development Laboratories, Hughes Aircraft Company, 
Culver City, Calif. 1 Morikawa, George K., Supersonic Wing-Body Lift, presented at the ny 
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October 93, 1950 Annual Summer Meeting, I.A.S., Los Angeles, July 12-14, 1950. (To be 2A 
published in the JouRNAL OF THE AERONAUTICAL SCIENCES.) Struct 


2 Spreiter, John R., Aerodynamic Properties of Slender Wing-Body Combi 1M 


— REPRESENTATION OF center of pressure for a general. sks at tibessic, Teensic ond Sutera Sucks, HALA. C2 ie 
family of wing-body configurations (circular cylinder body 1662, July, 1948 . 
with no afterbody and plan form made up of straight edges) in a 
stationary linearized supersonic flow may be considered simply 
for the limiting cases given by ‘“‘slender body theory” and “‘strip ° 
theory.”” The results of reference 1 imply that the limiting 
cases yield the bracketing limits of the center of pressure. Results Airplane Pitching Moment of Inertia Bas 
are presented here for the clipped-delta family of wing-body 
combinations ae ee are the diameter-span lnnes Bouton J, Le 
—_ d as o/s (where 0 of P = 1), the ratio of the tangent of the Chief Loads Engineer, McDonnell Aircraft Corporation, St. Louis . 
trailing edge angle (a, cf. Fig. 1) to the tangent of the Mach October 24, 1950 . 
angle, a; (where a; > 0), and the ratio of the tangent of the trailing vid 
edge angle ) to the tangent of the leading edge angle (wo , ° 
8 gle (w ; S aad dette < 8 é (ox ) HE PITCHING MOMENT OF INERTIA, J,, of an airplane is needed | 
a:/m (where 0 < ai/m < 1). The center of pressure C.P./C.P.o : ‘ : : ae A 
: . : 5 to evaluate horizontal tail loads in maneuvers with pitching 
is expressed as a ratio where C.P.o is the center of pressure ‘ - s ; 
; : accelerations. Since the computation of J, usually follows weight 
for the exposed wing alone referred to the root chord leading ‘ / ? és that, 
a estimates of the detail parts of the airplane, the value of J, is see 
edge. p ‘ : : supe 
es " = not available in the early stages of an airplane design when the , 
For a; = 0 (‘‘slender body theory,” cf. reference 2), : a ‘i ; of t 
tail loads are needed for design. Therefore, a reasonable estimate ni 
Ce a 1 ' 3 r42 ny 3 1) of J, must be made from the available preliminary data. ewe 
C.P.o 1—A 3 ~ Ae eK \ A correlation between radius of gyration in pitch, k,, and C 
: : é : DB 
gee airplane length, from which the moment of inertia may be ; 
which is independent of the parameter a;/m. , drag 
F (“strip th ” obtained t on ti computed, has been established for several classes of airplanes 3 
for a; = & (“‘s eory,”’ obtained by simple ‘grations), ; : 2 ‘ oe } 
: : y pe er The correlation formula is ky = Cla, where ky is the radius of we 
' 1 far gyration in feet, C is a constant depending on type of airplane, as 
CP. ae and /,4 is the overall length of the airplane. Table 1 notes the - 
—_— ok 7 | Xx values of C determined from data on almost 100 different aircraft. - 
C.P 6 1 fa: \? pury 
3\m - Fj 
2 TABLE 1 sider 
i. a: \* wt ~ z = “a 
1 /a\? (“:) Type Engines No. of Engines C -_ 
1—-{- +rA- , (1 — A*) + 2 In A] Turbojet 1 or 2 0.172 pe 
3 \m (1 — A)! é; Piston 1 0.176 that 
- > ) : 6 
Piston 2 0.163 effec 


Piston 4 or more 0.158 Thi: 


(2 
at 
lfa : (“:) 1 1 
padi Gliders 0.19 
1 — “ r -((1—A F 
[ Ata] Pee aacgy EM a 7 





2 of t! 
Eqs. (1) and (2) are plotted in Fig. 1 and indicate that the The standard deviation, o, between the theoretical value the 
deviation of the center of pressure for the wing-body ccmbina- determined with the constants above and the actual detailed pars 
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moment of inertia computations is approximately 20 per cent. 
The extremes of individual departures of the theoretical from 
the actual are plus 63 and minus 31 per cent, and these are 
yalues for extremes of configuration. A single value of C, 
9.165, could be used for all classes of airplanes with the standard 
deviation only increasing to 22-per cent. The same single value 
for C is also obtained in reference 2. Improvement in the 
accuracy of individual estimates may be obtained by considering 
the location and relative size of large masses, such as the engines, 
fuel tanks, and turrets, and by adjusting the radius of gyration 
up or down slightly when the large masses are at the extreme 
ends of the airplane or near the c.g. 

The method presented above is similar to the methods of 
Reference 3 presents a method that relates 


references 1 and 2. 
the pitching moment of inertia only to the weight as I,*/s 
0113W’. This latter method is of questionable validity and 
should be used with extreme caution in the design of airplanes, 
since the standard deviation for the 100 airplanes available is 
87 per cent and the individual extremes are plus 346 and minus 
68.7 per cent. 

The McDonnell XF-85 parasite fighter illustrates the essential 
fallacy in predicting the moment of inertia without considering 
airplane length. This airplane is an unconventionally short 
airplane for its weight, so the method of reference 3 estimates 
an], that is 316 per cent too high. The method presented above 
predicts an J, that is 26 per cent low. 
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Base Pressures at Supersonic Speeds 


J. Leith Potter 

Assistant Professor, Department 
University of Alabama 

November 16, 1950 


of Aeronautical Engineering, 


RECENT PAPER! PRESENTS an analysis of drag data for various 

two- and three-dimensional shapes. The author concludes 
that, depending upon body shape, the base drag coefficient at 
supersonic speeds assumes an approximately constant percentage 
of the theoretical limit, 2/y/?. After considering the data 
available to him, the author suggests that Coz = 0.7/M? for 
typical projectiles. Furthermore, he proposes to use a curve of 
Cos/Cop,, wr. aS a function of (APB/Q)svsonic tO estimate base 
drag coefficients in supersonic flow. The purpose of the present 
Writer is to expand the discussion of base pressures by reviewing 
some additional data? on simple bodies of revolution and circular 
discs. It will be shown that the processes suggested in reference 
1 may easily lead to excessive errors even when used only for 
purposes of rough estimation. 

First consider Fig. 1, where it will be noted that there is con- 
siderable effect on Cpg due to variations in body length and 
boat-tailing. No single function of Mach Number alone can be 
expected to satisfactorily express Cpg. In general, it would seem 
that a method for estimating Cpgz should take into account the 
effects of Mach Number, body geometry, and Reynolds Number. 
This is proved by reference 3 and several unpublished Naval 
Ordnance Laboratory reports? on the subject. The importance 
of these factors seems to have been recognized in reference 1, but 
the final recommendations as expressed in Figs. 13 and 14 and 
Paragraph 4 of the ‘“‘Results’”’ do not clearly indicate how they are 
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to be taken into account. It might be added here that all the 
available data are not plotted on Fig. 1 of the present comments, 
since the general trend is sufficiently illustrated by the points 
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shown. Also, the points for missiles having no boat-tails may be 
considered corrected for sting effects. Sting diameter for the 
boat-tail investigations was 1.27 cm., and the corresponding cor- 
rection has not been determined as yet. 

Reference to Fig. 2 will prove that it is risky to estimate base 
drag on the basis of (APB/Q)eisonice If (APB/g) 
—(.44 is taken‘ for a disc at high Reynolds Numbers, one may go 
to Fig. 14 of reference 1 and find Cpg/Cppg = ().7 at super- 


subsonic 


mar. 
sonic speeds. However, this value could easily be several 
hundred per cent in error, since the data of Fig. 2 in this dis- 


on s 
mar. rom 0.26 to 0.70 for 


cussion indicate a range of Cpg/Cppg 
Mach Numbers between 1.56 and 5.05. It is noted that 
(AP B/) suteonic Fig. 9” on Fig. 14 of 


reference 1, but a value of —0.5 is given for a simple disc on the 


—0.7 is labeled ‘‘disc 
referenced Fig.9. This is unexplained. 

One may conclude that there is a need for much more analytical 
and experimental investigation before base drag coefficients can 
be predicted with reasonable accuracy for many common aero- 
dynamic bodies. Typical data such as are presented in Fig. 1 
definitely establish the importance of boundary-layer flow in 
determining base pressures. 
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A Simple Graphical Solution of the Duhamel 
Integral 


E. V. Laitone and J. T. Ahlin 

Associate Professor, University of California at Berkeley, and Research 
Engineer, Santa Monica Plant, Douglas Aircraft Company, Inc., 
Respectively 

November 17, 1950 


SUMMARY 


A graphical solution of the Duhamel Integral is presented which will 
quickly determine any airplane motion caused by a given control surface 
deflection. This method also solves the inverse problem, determining the 
required control surface deflection for any desired motion. It also deter- 
mines the unit response to a step deflection of the control surface if the re 
sponse to any arbitrary control surface pulse-type deflection is known 


DISCUSSION 
F THE SO-CALLED INDICIAL RESPONSE to a step increment of the 
elevator deflection 6 is known, then the time variation of the 


response g to any arbitrary elevator deflection can be obtained 
from the Duhamel Integral as! 


te t di(r) 
g(t) = q(t — r) di(r) = q(t — 7) dr (1) 
0 0 dr 


In reference 1, Jones obtained a graphical solution of Eq. (1) 
by decomposing the time axis into equal increments. This ap- 
proach results in an extremely tedious and time-consuming 
graphical integration. Because of these difficulties, another 
numerical method has been proposed that does not use the 
Duhamel Integral; this is the A.A.F. Linear Increment Tabular 
Integration method of reference 2. This method depends on an 
approximation of the elevator motion by superposition of linear 
elements after obtaining the response to a linear elevator deflec- 
tion from the equations of motion. 


approximation, the approximation must be subtracted out and 


Thus, at the end of a linear 
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the succeeding one added in. A lengthy tabular integration is 
then necessary before the airplane motion is determined. More- 
over, this method is not readily applicable to recorded elevator 
motions that are highly irregular. 

The first form of the Duhamel Integral suggests a simpler 
graphical solution by a decomposition of the elevator deflection 
into equal increments whereby the elevator deflection 6 rather 
than time ¢ = 7 is taken as primary variable of integration. It is 
observed that Eq. (1) is approximately 


g(t) = x q(t — 7) Ab. (2 
As an example of this method, the pitching motion of an air- 
plane is obtained using a representative indicial response and 
elevator deflection time history. 
From the equations of motion for pitching, as given in refer- 
ence 2, the response g to a unit step in elevator angle 4, is 


a= = 1 — e~ (6/2) c+ (bc/Zw) + (¢/Zw)* 
i= 
c ay c — (b?/4) 
2 - — (b2/4) 
sin 4 ,./c — t+ tan! [ te ” \) 
UN 4 | (b/2) + (c/Zw) 
with 
Le (Uo + Zq) | 
b= — — + M, - - Mwy 
; Ss Lw ; i= Zu 
and 


c = [1/(1 — Zy)] [ZwMg — (Uo + Zg) Mw] 


The aircraft pitching velocity g response to a 1° elevator step 
(an instantaneous 1° elevator deflection from its trim position at 
time ¢ = 0) is plotted in Fig. Ia. 

If the plot of gq; in Fig. Ia is cut out and used as a template, 
then the graphical construction for an arbitrary elevator motion, 
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as in Fig. Ib for example, proceeds as a simple summation of the 
template ordinates for equal increments of Aé, as shown in Fig. 
Ie, 

The example demonstrates how the response can be obtained 
from an arbitrary elevator time history by summing the proper 
template ordinates at any desired time T to obtain g(t) at that 
selected time. To accomplish this, the elevator deflection is de- 
composed into equal increments (A6,), and the step function re- 
sponse g, for Aé,, is applied successively at the points h, to, ..., 
tk, where tz < T, each time adding the new ordinate to the one 
previously determined by placing the bottom of the template 
on the previous point and marking the intersection of the top of 
the template with the vertical line from 7. The points (tH, 
h,..., tk, ..., ty) correspond to the midpoints of each Aé, in- 
crement. If the elevator deflection is negative, the correspond- 
ing step function increment naturally should be subtracted from 
previously determined points. 

The usually difficult inverse problem of determining the ele- 
vator deflection 6(¢) required to produce a given motion g(t) can 
be solved readily by using the indicial response template and a 
simple graphical construction. The example in Fig. II illustrates 
this method. 

In this case, the airplane response is divided into equal time 
increments (dj, Ao, . .., Am), and vertical lines are drawn from each 
to intersect the curve g(t). The template is now moved hori- 
zontally along g(t) = O until its top contour passes through the 
intersection of g(t) and \;._ The origin of the template then gives 
the time (¢,) corresponding to the midpoint of the elevator step 
increment required to produce this response. At each of 2.2, 
Ms, .. . , Am, the intersection of the template profile is marked, 
thus indicating the airplane motion due to the step increment 
att = ¢,. 
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The template is moved vertically until the bottom rests on the 
point just marked on d». Then it is moved horizontally until the 
contour passes through the intersection of g(t) and A». The ori- 
gin of the template now gives the abscissa t = fz. Keeping the 
origin on the vertical line ¢ = f2, the template is moved vertically 
until the ordinate contributions due to the step at t = & have 
been added to each ordinate previously marked on dg, Ay, . . . , Am. 
This process is continued until the steps determined produce the 
final airplane motion. 

The amplitude of each equal step is determined already by the 
step function response. Hence, the cascade of steps shown in 
Fig. II is known. A smooth curve giving the required elevator 
time history can then be found by fairing through all the equal 
elevator step increments. 

The practical problem of finding the indicial step from any 
measured airplane motion due to an arbitrary elevator time 
history can be obtained by trial and error, using elements of this 
graphical method. After the elevator motion has been repre- 
sented by step approximations for equal Aé increments, the 
problem reduces to one of determining the template necessary 
to produce the given airplane motion. This is most easily ac- 
complished by attempting to reproduce the airplane motion with 
a trial template and correcting this template whenever necessary. 
A linear template should be used as the first approximation. 
Then the time abscissa y; is determined, beyond which the tem- 
plate fails to yield the correct airplane response. At some later 
time abscissa y2, any discrepancy between the result obtained by 
using this trial template and the true curve should be alge- 
braically subtracted from the corresponding trial template ordin- 
ate at t = y2. This new point should be connected with the 
ordinate at y; by a straight line, thus yielding a new trial tem- 
plate. A continuation of this iteration process determines the 
indicial response curve to any desired accuracy by a series of 
straight lines. 

The template method has proved simple, rapid, and accurate 
in practice. However, it is felt that the versatility of the method 
is its prime asset, because of the many applications and simpli- 
fications that appear. 
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—— HAS BEEN developed by means of which the pressure 
and density fields behind blasts passing over arbitrary sur- 
faces can be obtained in explicit analytic form. Assume that 
an originally plane shock front of given intensity (strong or 
weak) advances over the surface f(z) of a given thin airfoil or 
other flat structure. If the inclination of this surface with regard 
to the direction of propagation of the shock is relatively small, 
the advancing front will end perpendicularly on the surface at 
any instant, so that a curved shock front results, as well as a 
shock configuration that may be interpreted as Mach reflection. 
Because of the curvature of the shock, the flow behind (the front) 
is rotational and, for sufficiently strong blasts, is supersonic with 
regard to the airfoil. Such a case is treated here as a time-de- 
pendent perturbation problem, the undisturbed system being the 
advancing plane shock front and the disturbance being caused by 


the inclination of the surface f(z). 














8 
ISOPYCNICS a 
8 


DASHED LINES: ISOBARS a 


WA 


SOLID LINES: 


SHOCK FRONT: #=0.51 
AIRFOIL SURFACE* 2 = 0 


L.E. at —-& =-1.2I 
ct 




















Fig. a. 





Taking into account the rotationality of the flow, it turns out 
that the disturbance pressure 6p satisfies the wave equation. In 
addition, the usual boundary condition on the airfoil (or other) 
surface and the conservation laws applied across the disturbed 
shockfront must be fulfilled. Fixing (x, y, ¢) in the undisturbed 
flow behind {the shock) and subsequently introducing (%, 7, 2) 
obtained from the first system by means of a Lorentz transforma- 
tion, the following analytic problem may be formulated: 
Differential Equation: 

A(ép) — (1/c*?)(ép)iz = 0 

Boundary Conditions: 

(Surface) 

At 7 =0,% <0: (8p)y = Auf’ {al[cd + (¥/do)]} 
(Shock Front) 
Aty>0,% =0: D®(sp) =0 

where DD is a differential operator, of the second order, in # 
and / and where the parameters Ao, a, and \» are simple (rational) 
expressions in terms of the shock velocity U) relative to the airfoil 
and the velocity U and Mach Number of the shock relative 
to the undisturbed flow behind. 

It is now an essential point of this theory that, using some ‘“‘lost 
boundary conditions,” (3) can be replaced by the condition: 


(Reflected Surface) 
3 
Aty=0,%>0: (dp)y = D> Aif’{alcl — (%/a)]} (4) 
i=1 


where again the A; and }; are simple expressions in Uy, U, and 
M. 

This step makes it possible to solve the boundary value problem 
in terms of a Possio solution: 


5p(x, §, 2) = 
( 
4 if fe — (5p),(£, 0, 7) (5) 
7 Vi — 7)? — (1/c?)[((# — &)? + 57] 


After, substituting for (6p), from Eqs. (2) and (4), the integration 
is carried out over all values of £ and r for which the denominator 
is real. Expression (5) for 6p is valid for any surface f(x). This 
result is not restricted to conical fields as were the results of pre- 
vious investigations along similar lines. 

In order to compare with shock tube experiments, it is neces- 
sary to compute also the density field. Note that the ‘isopyenies 
that appear in any interferogram will differ from the isobars in 
rotational flow. It turns out that here the rotational flow is 
confined to a triangular domain 70C, immediately behind the 


shock front. One obtains 


1 
bp = — {5p(x, 9,2) — 2-5p(0, 9, 2 — (%/cM)]} (6) 
- 


For a diamond profile, Eq. (5) can be evaluated explicitly. 
Using physical coordinates, one obtains simply 


ct + Aox 
V (rorct + x)? + (1 — Av®)y? 
3 : ct — Ux 
> K;-cos™! 
i=1 


/ : 

V (Act — x)? + (1 — A?) 
and a similar expression for 5p. If the shock is behind the maxi- 
mum thickness point, three more cos~! terms have to be added 
in Eq. (7). Again, 2 and K; are simple expressions in Uo, U, and 


M. 


In Fig. 1, the curves of constant pressure and density are plot- 
ted for (Uy) — U)/c = 1.21, in accordance with Eq. 7 and the 
similar expression for dp. It is interesting to note how strongly 
and suddenly each isopycnic splits off from the corresponding 
isobar after traversing the demarkation line between irrotational 
and rotational flow. The point 0 is a singular point, in which a 
large number of isopycnics ends. Within a thin strip along 07, 
there is a steep density gradient, while there is no appreciable 
In other words, a slipstream is obtained. 


6p = Kocos™! 


(7) 


pressure gradient. 
Such result must be expected with every Mach reflection and can 
be clearly seen in W. Bleakney’s interferograms. 
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